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. Résumé. — On décrit deux expériences de diffusion des neutrons par des substances magné- 
tiques. Dans la première, on étudie au moyen d’un spectromètre àlcristal la structure antiferro- 
magnétique de FeCl, et l’action d’un champ magnétique extérieur sur la disposition des moments 


magnétiques. : 


La seconde, effectuée avec un spectromètre à temps de vol, a pour objet la diffusion critique au 
voisinage du point de Curie du fer. Elle permet de déterminer les deux paramètres caractérisant 
les fluctuations d’aimantation, la partie des corrélations «et leur constante de diffusion. 


Abstract. — This paper describes two series of experiments on neutron scattering by magnetic 


media. 


Series À gives the antiferromagnetic structure of FeCl,, and the effect of a magnetic 
field on the spin orientations as measured with a crystal spectrometer. 


In series B, the critical 


magnetic scattering of iron near the Curie Point has been investigated with a time of flight spectro- 


meter. 
which is related to their decay in time. 


4. Les neutrons qui sortent d’un canal de pile 
sont en équilibre thermique avec la matière dans 
laquelle ils se sont ralentis. Leur répartition spec- 
trale est donc maxWellienne et la longueur d’onde 
correspondant au maximum dulux est de 1,8 À. 
Cette longueur d’onde, du même ordre de grandeur 
que celle des rayons X utilisés pour l’étude de 
structures cristallines, permet de faire des études de 
 cristallographie en vue de déterminer la position 
des atomes dans la maille cristalline. De plus les 
neutrons possèdent un moment magnétique égal 
à 1,91 magnéton de Bohr nucléaire, qui peut 

interagir avec les moments magnétiques atomiques. 

Par étude de la diffraction des neutrons, nous 

pourrons donc localiser les porteurs de moments 
. magnétiques dans la maille cristalline. C’est là un 

des avantages principaux des neutrons sur les 
rayons X et qui justifie leur emploi, malgré les 
difficultés rencontrées (rareté des neutrons — inten- 
sité des flux relativement faible — impossibilité 
_ d’une détection sur film). 
Dans les conditions usuelles, c’est-à-dire pour des 
neutrons non polarisés et des noyaux non orientés, 
l'intensité neutronique diffusée est la somme de 
l'intensité diffusée par les noyaux et de l'intensité 


As a result one gets both therange of (he spin correlations and the scattering constant 


provenant de la diffusion magnétique, sans qu’il y 
ait entre elles d’interférences : 


TANT (1) 


L’intensité «nucléaire » diffusée par un éristal a 
pour expression 


IX = AÏEb; exp (:KR;)|°. (2) 

3 
La sommation est étendue à tous les atomes, 
repérés par le vecteur R,. À est une constante et K 


le vecteur de diffusion 


K—on(®+) = Xe (3) 


où A,est la longueur d’onde des neutrons incidents 


dont la direction de propagation est repérée par le 
vecteur unitaire S$,; À et s se rapportent aux 
neutrons diffusés. 

b; est l'amplitude de diffusion cohérente du 
noyau j. C’est une constante, ayant les dimensions 
d’une longueur, telle que la section efficace de 
diffusion cohérente soit 


* Gj-— LT b;,. (4) 
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L’intensité « magnétique » diffusée a une expres- 
sion analogue à (2), en remplaçant les amplitudes 
de diffusion nucléaires par les amplitudes de dif- 
fusion « magnétique » bf” qui dépendent : 


a) de la longueur du vecteur de diffusion K, ce 
qui traduit les dimensions finies des couches élec- 
troniques responsables de l’existence des moments 
magnétiques atomiques, d’où un facteur de 
forme f;(X) ; 


b) de la direction du moment magnétique M; de 
l’atome 7 par rapport au vecteur de diffusion K. 
Pour un neutron de spin s, la longueur de diffusion 
magnétique de l’atome j est 


= bo fi(E) [(eM;) e — M;].s (5) 
; a RCE 
où best une constante égale à Eu (y = rapport 


gyromagnétique du neutron). Dans le cas où tous 
les moments magnétiques sont parallèles ou anti- 
parallèles, l'intensité magnétique est proportion- 
nelle au carré du sinus de l’angle que font les 
moments magnétiques avec le vecteur de diffusion, 
comme on le voit facilement sur la formule (5). 
L'étude de l’intensité « magnétique » pourra donc 
nous permettre de déterminer la position et la 
direction des moments magnétiques atomiques. 
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2. Techniques expérimentales. — Les techniques 
utilisées dans la physique des neutrons n’étant pas 
encore très classiques en France, nous nous per- 
mettrons de rappeler très brièvement les dispositifs 
expérimentaux utilisés. Le but à atteindre est : 
1) Isoler, dans tout le spectre des neutrons sortant 
d’un canal de pile, un faisceau de neutrons mono- 

 cinétiques  (monochromatiques) de longueur 

d’onde À, et d'énergie Z,. 2) Etudier : a) la répar- 
tition des neutrons diffusés en fonction de l’angle 
de diffusion soit n(0) ce nombre ; b) le spectre 
énergétique des neutrons diffusés dans une direc- 
tion donnée, (0, E). 


À. SPECTROMÈTRE A CRISTAL. — Le plus simple, 
pour obtenir des neutrons monocinétiques est 
d'utiliser un monochromateur à cristal, analogue, 
à l’échelle près, à ceux que l’on emploie dans le cas 
des rayons X. Le faisceau de neutrons, collimaté 
dans le canal de la pile, tombe sur un monocristal 
dont une famille de plans réticulaires distants de d 
fait un angle 0 avec l’axe du faisceau. Les neutrons 
réfléchis sélectivement ont une longueur d’onde 
à = 2d sin 0. Il suffit d'éliminer tous les autres 
neutrons, ce qu’on réalise au moyen d’une lourde 
protection percée d’un canal. Dans les expériences 
de diffraction cristalline, les neutrons diffractés ont 
en majeure partie la même énergie que les neu- 
trons incidents (diffusion élastique) et pour les 
étudier, il suffit d'explorer les différents angles de 
diffusion au moyen d’un détecteur (généralement 
un compteur à BF; enrichi en isotope !B). Les 
figures 1 et 2 montrent un plan schématique et une 
vue d'ensemble du spectromètre utilisé à Saclay [1]. 

. Si on veut déterminer la répartition énergétique 
des neutrons diffusés par l’échantillon dans une 
direction donnée, il suffit d'analyser le flux diffusé 
au moyen d’un second spectromètre à cristal, 
comme le montre la figure 3. Maïs un spectromètre 
à cristal est un appareil dont la luminosité est très 
faible, et ce dispositif ne peut être utilisé que si le 
flux de neutrons sortant de la pile est suffisamment 
élevé, comme c’est le cas pour la pile canadienne de 
Chalk-River [2]. 


B. SÉLECTEUR A TEMPS DE VOL. — Lorsqu'on ne 
dispose que d’une pile de flux moyen, il est néces- 
_saire de recourir à un dispositif expérimental dif- 
férent, le spectromètre à temps de vol, dont le 
principe est assez simple, mais qui nécessite un 
appareillage mécanique et électronique très com- 
plexe. | è 
Considérons un neutron diffusé à l'instant 4 — 0, 
et dont la vitesse est v. Il atteindra le détecteur, 
placé dans une direction donnée, à une distance ZL 
de l’échantillon à un instant & — L/]v. Le faisceau 
de neutrons est haché par une roue à bras haute L, 
percée d’une fente radiale tournant à une vitesse 
angulaire ©, de lordre de 6 000 tours/minute. 


ÉTUDES DE.MAGNÉTISME RÉALISÉES AVEC DES NEUTRONS 


619 


L’instant origine est déterminé par le passage de la 
fente devant le faisceau. Ce passage déclenche un 
ensemble denumérateurs retardés. Ces numérateurs 
sont au nombre de V (NW — 100 dans le spectromètre 


Crisle/ 


monochromaleur 


Diffusaur 


Cristo/ 


ono/yseur 


Daltecleur 


TGS. 


£chanillon 


Complage 


Ex AAA} 


A 


Sélecteur en 7 


temps _ 7 x 


Fic. 4. 


utilisé à Saclay) et le numérateur n totalise les 
neutrons qui ont atteint le détecteur à un temps f 
compris entre (n —1)r et nT; c’est-à-dire les 
neutrons dont la vitesse est comprise ertre Lnr 
et L(n — {)r. x est une constante de temps de 
l’ordre de 10 us. Grâce à ce dispositif, le spectre des 
neutrons diffusés est directement enregistré dans 
les N numérateurs. 

Pour monochromatiser les neutrons, on utilise 
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une roue L:,, analogue à la roue L,, placée avant L, 
sur le trajet des neutrons sortant de la pile, la dis- 
tance séparant les deux roues étant /. La fente 
de L,, est décalée d’un angle « par rapport à celle 
de L.,. Les neutrons de vitesse v, = lw/a passent 
à travers les fentes des deux roues ; les neutrons 
plus rapides tombent en avant de la fente de L:, 
les neutrons plus lents tombent en arrière. La 
figure 5 représente la forme de la raie monochro- 
matisée ainsi obtenue. Sa largeur est due à la 
largeur des fentes ét au manque de parallélisme du 
faisceau de neutrons collimaté à l’intérieur du 
canal de la pile. 

Pour augmenter la luminosité de l'appareil, 


35 


45 5 


Fre. 5. 


chaque roue porte plusieurs fentes. Une vue 
d’ensemble du sélecteur mécanique installé devant 
la pile de Saclay a été publiée dans ce journal [3]. 


3. Structure antiferromagnétique du chlorure 
ferreux anhydre. — Nous avons entrepris des 
expériences sur le chlorure ferreux anhydre en vue 
de mettre en évidence l’action du champ magné- 
tique sur la disposition des moments magnétiques. 
En effet, les mesures magnétiques ont montré 
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que FeCl, était à basse température un antiferro- 
magnétique d’un type particulier, appelé méta- 
magnétique : son couplage antiferromagnétique est 
assez faible pour que, sous l’action d’un champ 
magnétique modéré, ses moments magnétiques 
soient parallèles. 


A. STRUCTURE MAGNÉTIQUE DE FeCl,. — FeCl, 
cristallise dans le système rhomboédrique. Sa struc- 


oui? 
101102 107 dé 105 


te AN RES 


Spectres de diffraction de neutrons de FeCl, 


Miciut 


ture est représentée sur la figure 6 avec la maille 
hexagonale correspondante. C’est un empilement 
de feuillets d’atomes de fer, séparés par deux 
feuillets d’atomes de chlorure, chaque maille cristal- 
lographique contenant trois feuillets de fer. La 
figure 7 montre les spectres de diffraction obtenus 
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à la température ordinaire et à la température de 
Vhélium liquide. On voit qu’il apparaît à basse 


_ température quatre raies supplémentaires qui tra- 


duisent l'existence d’une surstructure magnétique. 
On interprète ces raies en admettant que ‘les 
atomes de fer d’un même feuillet ont leurs moments 
parallèles, les moments de deux feuillets voisins 


- étant antiparallèles. Il en résulte que la maille 


«magnétique » a un paramètre c double celui de la 
maille cristallographique. Dans cette nouvelle 
maille, les raies de surstructure peuvent s’indexer 
(1, 0, 1) où ! = 1, 5, 7, 11. En outre l’absence de la 
raie (0, 0, 3) impose aux moments atomiques une 
direction parallèle à l’axe c. Ces résultats sont en 
accord avec ceux obtenus par Wilkinson et 
Cable [5]. La figure 8 montre l’accord entre les 
intensités mesurées et calculées à partir de ce 
modèle, ces intensités étant normalisées à leur 
somme. 


CL 


@ / théorique 


CHE expérimentale 
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B. ÉTUDE DE L'ACTION D'UN CHAMP MAGNÉ- 
TIQUE. — La direction du champ magnétique reste 
parallèle au vecteur de diffusion. Pour tous les 
cristallites qui participent à une réflexion sélec- 
tive (k, 4, l), le vecteur de diffusion K coïncide 
avec le nœud (4, 4, l) du réseau réciproque. Tous ces 
cristallites ont donc même disposition cristal- 
lographique par rapport au champ magnétique. Si, 
de plus le cristal est uniaxe et si la direction d’anti- 
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ferromagnétisme À est l’axe c, pour tous ces cris- 
tallites À fait le même angle © avec le champ 
magnétique. Dans un cristal uniaxe on peut, au 
moins en première approximation, admettre que 
lorsqu'on applique un champ magnétique H les 
vecteurs aimantation des deux sous-réseaux M, 
et M, primitivement antiparallèles dans la direc- 
tion de À restent dans le plan formé par la direction 
d’antiferromagnétisme et le champ magnétique H. 
Dans ces conditions les vecteurs M, et M, des 
divers cristallites contribuant à une même raie 
(2, k, l) sont parallèles entre eux. Donc l’étude de la 
variation de l'intensité d’une raie, en fonction du 
champ est équivalente, aux interactions magné- 
tiques entre cristallites près, à une étude sur mono- 
cristal. Les mesures magnétiques effectuées sur une 
poudre donnent au contraire des résultats qui sont 
une moyenne plus difficilement interprétable. 

Dans un cristal métamagnétique, les aiman- 
tations M, et M, des deux sous-réseaux seront 
parallèles lorsque le champ magnétique sera supé- 
rieur à un certain champ limite Z;(+). L’intensité 
d’une raie magnétique est alors nulle. Mais la 
variation d'intensité en fonction du champ diffère 
suivant la valeur de l’angle ©. Cette intensité est en 
effet proportionnelle à 

(sin &«, — sin œ,)? 


où «, et «, désignent les angles que font, avec la 
direction du champ les moments M, et M, des 
deux sous-réseaux (cf. fig. 9). Or la variation des 


à S 


H1G:-9: 


angles &, et «, en fonction du champ dépend de 
l’angle ©. Une étude théorique élémentaire montre 
que le champ limite Æ5(œ) est d'autant plus élevé 
que l’angle @ est plus voisin de 900. Comme le 
champ magnétique dont nous disposions ne pouvait 
dépasser 15 000 Oe, nous avions intérêt à étudier 
des raies magnétiques correspondant à des angles @ 
petits. Mais l’intensité d’une telle raie, propor- 
tionnelle à sin? +, est faible. Nous avons adopté 
une solution intermédiaire, et nous avons étudié la 
variation de l'intensité de la raie (1, 0, 7) en fonc- 
tion du champ, pour laquelle l’angle ® est égal 
à 570, La figure 10 représente le résultat de cette 
étude. La courbe en tirets représente la variation 
calculée en supposant, comme l’a fait L. Néel, que 
l’anisotropie peut être représentée au moyen d’une 
seule constante Æ et que X/n M? — 0,25, où n est 
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la constante de champ moléculaire et 7 l’aiman- 
tation. Le champ critique a été ajusté pour obtenir 
e meilleur accord entre la courbe et les points 
expérimentaux. On trouve A, = 8 500 Oe, ce qui 


J50D0 


10000 
Fi. 10. 


5000 


est en assez bon accord, vu la précision des mesures, 
avec les résultats des mesures magnétiques de 
Starr Bitter et Kaufmann qui donnent 10 000 Oe 
comme ordre de grandeur de ce champ. 


4. Diffusion critique au voisinage du point de 
Curie. — Dans la structure magnétique que nous 
venons d'étudier, l’arrangement des moments 
magnétiques était bien ordonné. Nous allons main- 
tenant envisager le cas d’une structure magnétique 
partiellement  désordonnée, c’est-à-dire dans 
laquelle il y a de fortes fluctuations dans l’orien- 
tation des moments magnétiques. De telles fluc- 
tuations d’aimantation existent dans un ferro- 
magnétique lorsque la température est voisine de la 
température de Curie [7]. Ceci se traduit expéri- 
mentalement par une diffusion magnétique anor- 
male qui n’est pas sans analogie avec la diffusion 
de la lumière par les fluides au voisinage du point 
critique. La figure 11 représente la variation de 
l'intensité diffusée à 39 en fonction de la tempé- 
rature dans le cas du fer, au voisinage du point de 
Curie. Nous avons choisi d'étudier le fer parce que 
lintensité diffusée, proportionnelle au carré du 
moment magnétique, est beaucoup plus forte que 
dans le cas du nickel. Le but des expériences est de 
déterminer les paramètres au moyen desquels on 
peut donner une description des fluctuations. Ces 
fluctuations peuvent se décrire en termes de corré- 
lations entre les moments d’atomes situés à la 
distance R. Soit 


Y{R < M(0).M(R) > moy. (6) 


1 
Dear 
la fonction de corrélation correspondante. A très 
basse température, tous les moments magnétiques 
sont parallèles, et la fonction de corrélation est 
égale à 1. A très haute température, dans le domaine 
paramagnétique, les moments magnétiques sont 
indépendants, et la fonction de corrélation est 
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nulle, quel que soit À. À une température inter- 
médiaire, c'est une fonction de À présentant un 
profil en «toile de tente » tel que le représente la 
figure 12. On peut la décrire grossièrement par la 
portée de corrélations 1/*x. Si la distance de deux 


Intansitæ 
au pic 


Pre 42: 


atomes est supérieure à 1 /x, leurs moments magné- 
tiques sont indépendants, donc y — 0. Si leur dis- 
tance est inférieure à 1/x, la corrélation est diffé- 
rente de zéro. Notons que pour À = 0, y = 1, par 
définition de y(R). On peut montrer qu’asympto- 
tiquement la fonction a pour expression 


e—xR 
F5 


i 


Y(R) = C (7) 


No7 


Caleulons le carré de la fluctuation d’aiman- 
tation. On a 


< ES M(R)? > = ZE < M(R) MR) > 
R R R’ 
= NE <M(R).M(R + R') > = N.M°>7y(R). 
R’ R 


Donc l'aire délimitée par la courbe y(A) est pro- 
. portionnelle au carré de la fluctuation d’aiman- 
tation. Cette fluctuation devenant très grande au 
voisinage du point de Curie, l’aire délimitée par y(Æ) 
devient très grande. Comme (0) = 1, on en déduit 


x (2° #2) 


10.10 


Tan °C 


820 


790 800 &0 
Fie49 


770 780 


que les ailes deviennent très grandes, donc que x 
devient très petit. Un calcul plus précis montre que 


ÉTÉ 


7) me Fe 


où a est la distance entre moments magnétiques. 
L’intensité magnétique diffusée dans une direction 
repérée par le vecteur de diffusion K est la trans- 
formée de Fourier de y(R) : 


I(K) = Z eiK.Ry(R). 
R 


À basse température, y(R) — 1. L’intensité est 
donc localisée sur des pics de réflexion sélectives 
_(K doit être égal à un vecteur du réseau réci- 
proque t). À très haute température, y — 0,iln’ya 
pas de diffusion magnétique. À une température 
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intermédiaire, on trouve en utilisant la forme 
asymptotique (7) : 
TR à 


R=j+ix : 


Nous avons étudié la diffusion en avant (tr — 0), 
le faisceau diffusé faisant un angle de 39 avec le 
faisceau incident. De la mesure de l’intensité dif- 
fusée, on peut déduire facilement la valeur de x. 
La figure 13 représente la variation de x avec la 
température. Ce résultat est en bon accord avec 
lexpression théorique de x donnée plus haut. 

L'expérience ayant été faite avec un spectro- 
mètre à temps de vol, on peut également déter- 
miner la dynamique des fluctuations. En effet, de 
même que Z(K) est la transformée de Fourier 
de y(r), la section efficace différentielle est la trans- 
formée de Fourier de la fonction de corrélation 
spatiotemporelle : ne 

TK 0) = fat E Y(R, 0) cilotK.R) (9) 
où Lo est la variation d'énergie du neutron el 
MR) n< M(0, 0).M(R, 1) >. (10) 

Les fluctuations sont un phénomène irréversible, 
eb il est raisonnable d'admettre que leur compor- 
tement est régi par une équation du type équation 
de diffusion : 

dY(R, t) 


——s, — AÀ.ARY(R, i) 


Si (11) 


À est un paramètre phénoménologique qui décrit le 
comportement dans le temps des fluctuations. 
D'une façon plus précise Ax? est l’inverse d’un 
temps de relaxation + qui est relié à la variation 
moyenne d'énergie par À © + fr —hAx?. 
L'équation (11) permet de déterminer la trans- 
formée de Fourier de Y(R, t), donc l'intensité dif- 
fusée. En pratique, la variation d’énergie est beau- 
coup plus petite que la largeur de la raie incidente. 
Il faut donc calculer la forme de la raie diffusée 
pour différentes valeurs de À, ce qui permet de 
déterminer la valeur de À qui cadre le mieux avec 
le résultat expérimental. Les premières mesures 
de A ainsi faites n’ont pu nous donner qu’un ordre 
de grandeur. Cependant les résultats semblent en 
meilleur accord avec les calculs faits en utilisant un 
modèle de Heisenberg pour le fer plutôt qu'avec 
ceux qui utilisent un modèle de bandes. 


Manuscrit reçu le 8 février 1958. 
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CALCUL DES VALEURS PROPRES DE L'OPÉRATEUR DE COLLISION ÉLASTIQUE 
D'UN GAZ DE LORENTZ. ÉTUDE DE QUELQUES CAS USUELS 


Par A. ARSAC, L. BASQUIN, J. F. DENISSE, J. L. DELCROIX et J. SALMON, 


Commissariat à l'Énergie Atomique. 


Résumé. — Cet article contient le calcul des premières valeurs propres de l’opérateur de collision 
élastique d’un gaz de Lorentz pour différents types d’interaction. 
On développe, en particulier, les procédés numériques qui ont permis d’obtenir ces valeurs 


propres. : 


Abstract. — This study includes the calculation of the first eigenvalues of (he elastic collision 
operator in a Lorentzian gas for different types of interaction. 
Particularly, they expand the numerical processes which permitted these eigenvalues to be 


obtained. 


I. Introduction. — Le calcul des valeurs propres 
de l’opérateur de collision élastique électron-molé- 
cule revêt une grande importance en théorie ciné- 
tique des plasmas. Dans lapproximation du gaz de 
Lorentz, Bayet, Delcroix et Denisse [1] ont en effet 
montré que la détermination de la fonction de 
distribution des vitesses des électrons était étroi- 
tement liée à la connaissance de ces valeurs propres. 
Celles-ci se divisent en deux classes. D’une part les 
valeurs propres v, relatives aux anisotropies ; leur 
calcul est relativement aisé. D’autre part les valeurs 
propres À, relatives aux écarts non-maxwelliens 
caractéristiques pour la partie isotrope de la fonc- 
tion de distribution. Leur recherche conduit à des 
calculs numériques assez complexes. Nous avons 
effectué ceux-ci pour quelques types d'interaction. 


Nous avons choisi le cas d’un potentiel répulsif 


en 1 /r#, celui du potentiel répulsif en 1 /r° ainsi que 
le cas des molécules considérées comme des sphères 
rigides. (Potentiel en 1 /r®). 

Seul ce dernier exemple correspond à une inter- 
action ayant un sens physique. Il était alors extrè- 
mement intéressant de comparer les résultats 
obtenus à ceux provenant de l’étude d’un poten- 
tiel 1/r4 pour lequel il s’introduit des simplifi- 
cations mathématiques considérables conduisant 
à des expressions très simples des valeurs propres. 
Enfin le potentiel en 1 /r5 assurait une manière de 
continuité entre ces deux cas extrêmes. 

Désignant ces valeurs propres par les symboles 
— 9, — Ày, ... — À, nous rappelons qu’elles ne 
dépendent que du rapport m/M des masses de 
l’électron et de la molécule et de la fonction-repré- 
sentant la variation de la fréquence de collision 
avec la vitesse. Le type d'interaction choisi déter- 
mine complètement cette fonction. Désignant par T 
la température des molécules, par X la constante 
de Boltzmann, par v la vitesse des électrons et 
par v, la fréquence de collision, nous poserons : 


Do = [2:AT Im]? (1) 


Le DO, (2) 
V1 = Vof(u). (3) 


L’équation (3) signifie que v, est mis sous la 
forme du produit d’une fréquence de collision 
constante v, par une fonction sans dimensions {(u) 
de la vitesse réduite. 

Nous rappellerons maintenant que les valeurs 
propres sont données par une équation de Sturm- 
Liouville dont l'expression s’écrit : 


£ [out u (y + | a 


et à laquelle sont associées des conditions aux 
limites homogènes pour les valeurs nulle et infinie 
de u. 


ut y = 0 (4) 


Posant : 
y = e—w plu) (5) 
m 
A = UT (6) 
; u = 24 V9 (7) 
il vient : 


a [tu u? 24e El + &au?e—w 6 — 0 (8) 
a est maintenant une quantité sans dimensions ; 
d'autre part, type d'interaction non seulement 
détermine f(u) mais, de plus, impose des conditions 
aux limites homogènes à l'équation (8), conditions 
homogènes à partir desquelles, on peut calculer 
numériquement les valeurs propres &5, &4, ... & 
de (8). 


IT. Cas d’un potentiel d’interaction en 1/r*, — 
La fréquence de collision est alors indépendante de 
la vitesse et on peut écrire : 


f(u) = 1 S (9) 


des sserene tn 


N9:7 OPÉRATEUR DE COLLISION ÉLASTIQUE D'UN GAZ DE LORENTZ 62 


On montre [1], de plus, que les conditions aux 
limites s’écrivent pour un tel potentiel 


WW do/du —:0 (11) 
u — co p—> ui nie ip —+ 0. (12) 


L’équation (10) s'intègre analytiquement et pos- 
sède les valeurs propres et fonctions propres : 


ag = 0 Po — 2 (13) 
5 3 
CN À Pr 10 (u®—:) (14) 
4 r,,2 15 F 
0) LE ont Len DUT (15) 
HA CR qi = H si + 1(u) (16) 


u 


Hi x(u) représentant le polynôme de Hermite, de 
degré 21 + 1. Sur les figures 1 et 1 bis se trouvent 
à un facteur constant près, les courbes représen- 
tatives des fonctions propres @(u) et y(u) correspon- 
dant à ces 4 valeurs propres (!). 


IIT. Cas des molécules considérées comme des 
sphères rigides. — Le libre parcours moyen d’un 
électron est alors constant et la fréquence de col- 
lision devient proportionnelle à la vitesse. On peut 
donc écrire : 


= it (17) 

Tu = (18) 

L'équation de Sturm-Liouville prend la forme : 
& [ue eu SE] + Aau? eu" o — 10; (19) 


D'autre part, un tel type d'interaction conduit [1] 
aux conditions aux limites : 


d ï ! 
Hi 0 en Ku = 0 (20) 

do < 
Hé op + 0, (21) 


K désignant une constante qui sera déterminée plus 
loin. 

Il n’est plus possible de traiter le problème de 
manière analytique dans son ensemble. Aussi, 
avons-nous mis au point une méthode qui utilise 
à la fois un développement analytique près de 
l’origine, un procédé classique d'intégration numé- 
rique et une méthode variationnelle. 

Ecrivons (19) sous la forme : 


u(d? p/du?) + (3 — 2u?) (dofdu) + ka p — 0 (22) 
et, cherchons un développement du type 
4 e = u’[b, + bu... + bu? + .,.] (23) 


(1) Le facteur constant a été choisi, de manière à rendre 
e égale à l’unité pour u égal à zéro. 


[1 


Ÿ(u) 


Fic. 1. — Potentiel répulsif en 1/4. 
Fonctions propres œ(u). 


Yu) 


0,20 


0,1 


Fic. 4 bis. — Potentiel répulsif en 1 /r4. 
Fonctions propres y(u) = e—w* o(u), 
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valable au voisinage de l’origine. L’équation carac- 
téristique a pour racine : 


ru 1) r = — 12, (24) 


Le développement correspondant à r ——72 ne 
peut convenir car il conduirait à une valeur infinie 
de + à l’origine, ce qui est physiquement impos- 
sible. 

Il subsiste le développement relatif à r = 0. En 
le portant dans l’équation (22) on peut lPexpliciter 
sous la forme : Ë 


@ = 1 — (4a]3) u + (2a°]3) u? — (8a/45) (1 + a?) u° 
+ (a2/135) (19 + a)ut. (25) 


On peut lui associer la condition aux limites 
homogène à l’origine. 
(dedu) + (ka]3) ® = 0 


et la constante X de l’équation (20) est liée à la 
valeur propres cherchée. 

Quant à la condition (21), elle résulte du fait que 
la fréquence de collision croissant indéfiniment 
avec u, les électrons tendent à se mettre en équi- 
libre de vitesse avec les molécules. y doit donc 
tendre vers 


u— 0 (26) 


et 
et l’on a bien 


u — co o—1 dofdu — 0. (27) 

Pour obtenir une valeur propre, nous choisirons 
une valeur arbitraire de a et nous intégrerons par 
un procédé numérique l’équation (22). 

Partant de l’origine avec la condition (26), nous 
intégrerons pas à pas et nous examinerons le com- 
portement de do/du lorsque # tend vers l'infini. 
La valeur propre cherchée est la valeur de a pour 
laquelle la condition (27) est satisfaite. 

En pratique, la présence d’un pôle à l’origine 
dans l’équation (22) soulève une première difficulté 
pour l'intégration numérique ; une autre provient 
du fait qu’on ne peut aller jusqu’à l'infini. 

Nous obtiendrons, cependant, les valeurs propres 
en utilisant l’artifice suivant. 

Dans un intervalle faible 0 —u, (uy — 0,01), 
nous utiliserons le développement (25) et obtien- 
drons les quantités w, et «6 

u — Uy Po — Po, 4) 

A partir de ,, nous intègrerons pas à pas à 
l’aide de la méthode de Runge-Kutta (intervalle 
élémentaire 0,005) et cela jusqu’à la valeur maxi- 
mum um — 9. En effet, la fonction propre initia- 
lement cherchée est y avec : 


y = et p (5) 
et, au delà de 5, y est pratiquement égal à 0. Il est 


donc sans intérêt de poursuivre l'intégration. 
En procédant de cette manière, nous constatons 
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qu'effectivement, il existe des plages de valeurs 
de a pour lesquelles d/du devient positivement 
très grand au voisinage de 5 et des plages de valeurs 
pour lesquelles do/du devient négativement très 
grand. Il‘existe donc une suite de coupures qui 
définissent les valeurs propres &s, 4 ... a, (voir 


fig. 2). 


so P (M) 


F1G. 2. — Équation de Sturm-Liouville. 
Recherches de la 1'e valeur propre. 
Cas des sphères rigides f{(u) = u. 
Variation de ® en fonction de w pour 3 valeurs de 4. 
IL. a — 1,1685163. 


IL 0="1,16851297% 
ELA UT6S: 


g (0) =1 


Les calculs sont effectués à l’aide d’une machine 
électronique Gamma A. E. T. Bull. 


Première valeur propre : 
5 10, 


Deuxième valeur propre : 


(22 u 1e 
du 
1,168 5 — 4h 923 
1,168515 5 + 0,606 
1,168516 5 + 116 
Troisième valéur propre : 
à u ce 
du 
2,9191 a + ‘6 215 
2,51955 5 2 0,198 
2,5196 ÿ — 776 
Quatrième valeur propre : 
de 
(12 u TA 
4,089 5 — 896 
4,0891 5 — 485 
2,0893 5 — 0,342 


| 


bg, ben ue 2/87: 


nn. «2 
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F1c. 3. — Cas des sphères rigides. 
Fonctions propres œ(u). 
— 1,168515. 
(II) @ — 2,51955. 
(III) az = 4,08929. 


Au voisinage de u — 0 on à la fonction 


Courbes I’, 11’, III’. 


Sur la figure 3 on a représenté les fonctions 
propres correspondant à ces quatre valeurs propres 
ainsi que leur développement en série autour de 
l’origine (25). On a également fait figurer les fonc- 
tions propres y(u) (fig. 3 bis). 


IV. Contrôle par la méthode variationnelle. — 
D’après les propriétés bien connues de l’équation 
Sturm-Liouville (4), À est extremum du quotient 


d’intégrales. 
co o 
f vilu) u? (5) et lu 
m Jo du a 
2m œ Fe Fa) 
JL u? p?e—w du 


SX — 0. (30) 
 Rappelant que: 
Ai=umIM (31) 


Ge SERA PTE (32) 


1u) 


0,05 


Fic. 8 bis. — Cas des sphères rigides, 
Fonctions propres y — e—# œ(u) 
1) ay = 1,168515. 
2) a = 2,51955. 
3) az = 4,0891. 


il vient également. : 


ee  (d?\* 
[ f{u) u (5) du 


d' — RS (33) 
ES > 3 3 
l Au? pe du 
û 


OURS (34) 


Il est donc extrêmement intéressant de calculer 
l'expression (33) au fur et à mesure de l'intégration 
numérique de l’équation différentielle. 

On choisira done toujours une valeur arbitraire 
de a mais on calculera simultanément œ(u) et 
(do/du) par l'intégration numérique et l’expres- 
sion (33). Celle-ci aura pour valeur 4,. Les valeurs 
propres seront les valeurs de a pour lesquelles a se 
confond bien avec a,. 

Bien entendu, le caractère approché des 
méthodes de calcul ne permet pas d’avoir un accord 
exact, mais il est intéressant de constater qu’en 
utilisant une méthode de Simpson pour (33) la 
différence relative entre & et a, est de l’ordre 


de 2/1 000 à 4/1 000. 


On obtient finalement le tableau suivant : 


Première valeur propre : 


ao = ÙÔ 


>, 
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Deuxième valeur propre : 


Valeurs Intégration Méthode Écart 
proposées numérique variationnelle relatif , 
“) ge: dy — 4 
a (5 (25) y | = | 
1,168 — 44 923 8,9697 
1,168515 0,60 1,170861 2/1 000 
1,168516 + 116 1,170917 
Troisième valeur propre : 
Valeurs Intégration Méthode Écart 
proposées numérique variationnelle relatif 
d @ dy — 4 
a (5) (u —= 5) dy | | 
2,5191 — 6 215 3,2820 
2,51955 7.019408 9,5280 3/1 000 
2,91 96 — 776 2,5399 
Quatrième valeur propre : 
Valeurs Intégration Méthode Écart 
proposées numérique variationnelle relatif 
dp ! y — à | 
(5) (u—=°5) dy | See 
4,0893 +. 16,98 4,10669 
4,0891 — 485,5 k,1198 4/1 000 
k,08929 — 0,34 4,10667 


Finalement, nous proposerons les valeurs : 


GÉE=A dy —1,17 ln 02502 d> — 4,09. 
Il est, bien entendu, possible d'obtenir une préci- 
sion supérieure à 3/1 000 mais les temps de calcul 


deviennent prohibitifs. 


V. Cas d’un potentiel d’interaction en 1 /r°. — 
Il est bien connu qu’une loi de force en 1 /r° conduit 
à une fréquence de collision variant selon une loi 
en "avec 

s — 5 


SE 


s —1 


Dans le cas présent, on aura, par conséquent : 


SARA AU Tati (34) 


et l'équation de Sturm-Liouville prendra la forme : 
1 14 
+ [uv re = + 4au?e—# © = 0, (35) 


Nous chercherons encore un développement 
limité au voisinage de l’origine. Pour ce, nous 
poserons : 

us (36) 
d’où 


10 æ) + 100 ax? @ = 0 (37) 
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et de là, les expressions : 


RAS than DT et 109 a° a7 | 
pére 1 10 a AO ES ET 
DD 2 EM OS 
œ{u) = 1 10 15% 
a? 106 4 
104 186 22107 CORP ATEN 
107 58 320 * ST 063 120 4 0e ES) 


F1G. 4. — Potentiel répulsif en 1 fr5. 
Fonctions propres o(u). 
I. a; = 1,066802. 
IL. a, = 2,21474. 
Ta = 3 41525 


Au voisinage de u — 0 on a la fonction : 


(12 
o(u) = 4 = 102 2 yo 
œ(u) 10 1% 
a? a? 
AO 00/6 AO 
F 1055 320 * 10 51 963 120 * 
Courbes I’; II’, III’. 


Nous adopterons, par conséquent, les conditions 
aux limites homogènes 


u'= 0 — = ,0 


(39) 


(40) 
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Fic. 4 bis, — Potentiel répulsif en. 1 /r5. 


Fonctions propres y) 


1) a; — 1,066802. 
2) dy = 2,218474. 
3) ay = 3,41545. 


eu, o(u). 


Nous utiliserons les mêmes procédés qu'aux para- 
graphes IIT et IVet obtiendrons le tableau suivant: 


a 


4re valeur 0 


[roses 


2€ valeur 1,066802 


1,066808 


3e valeur 2,21474 


2,21501 
ER 


EE 


4e valeur 3,41545 


3,43387 


d 
—. (u = 5) 
29 607 
+ 18,3 
+ 1 326 
HA 518. 
“ 2,58 
— 6 861 
— 11 389 
+ 428,78 
+ 80 274 


1,1492 
1,065% 
1,0657 
2,21412 
2,20946 
2,22043 


3,43914 
3,40398 


5,13953 


Aa 
a 


1/4 000 


2/1 000 


3/1 000 
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Finalement, nous proposerons les valeurs 


ag = 0 ai = 1,07 day — 2,21 ds = 3,42 
Sur la figure 4 on a représenté les fonctions 
propres correspondant à ces quatre valeurs propres 
ainsi que leurs développement en série au voisinage 
de l’origine selon la formule (38). On a également 


fait figurer les fonctions propres y(u) (fig. 4 bis). 


VI. Comparaison des résultats. — Le tableau 
suivant résume l’ensemble des résultats obtenus 


TYPE FRÉQUENCE VALEURS PROPRES 
DE COLLISION 

s n a 
0 

5 0 ’ É 
à 
3 

( 1/5 0 
4,07 
921 
3,42 

co 1 


pour des potentiels d'interaction en 1 /ravecs >5, 
c’est-à-dire, pour des fréquences de collision variant 
selon une loi en x" avec 0 £ n <1. 

Il est extrêmement intéressant de constater que 
les valeurs propres varient somme toute assez peu 
avec le type d'interaction, tout au moins pour les 
premières valeurs propres, car il semble que l’écart 
s’accentue peu à peu. 

Il serait enfin extrêmement intéressant de 
reprendre ces calculs pour une loi de fréquence de : 
collision correspondant à une forme d’interaction 
plus réaliste, liée à des considérations physiques 
précises, la fonction v,(u) ne se présenterait plus 
sous une forme analytique, mais il est sans doute 


possible de généraliser la méthode à de tels cas. 


Nous remercions vivement M. Cheneau pour la 
programmation et l’exécution des calculs. 


Manuscrit reçu le 12 février 1958. 
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ÉTUDE DE LA STRUCTURE DÉTAILLÉE 
DES COURBES DE DIFFUSION DES ONDES ÉLECTROMAGNÉTIQUES 
PAR LES SPHÈRES DIÉLECTRIQUES 


Par JEean MEVEL, 


Faculté des Sciences de Rennes. 


Résumé. — Des déterminations théoriques récentes de la courbe de diffusion. de la lumière par 
les sphères diélectriques, ont montré qu’elle présente, superposée aux ondulations principales, 
une série d’indentations secondaires dont l’origine n’a pas été élucidée. Dans ce qui suit, nous 
établissons que cette structure est due aux résonances des termes multipolaires constituant le 
rayonnement diffusé et nous interprétons chaque maximum secondaire. Nous avons dû pour cela 
formuler un théorème nouveau permettant de déterminer aisément l’intensité d’excitation des 
modes, nous en déduisons les positions et les intensités des maxima de la structure détaillée. Enfin, 
nous montrons que la courbe de diffusion actuellement connue est imparfaite et nous indiquons 
de quelle manière on doit la modifier. 


Abstract. — Recent theoretical computations of the light extinction curve by dielectric spheres 
have shown thatit has, superposed on the major oscillations, very rapid fluctuations of small ampli- 
tude : this feature has not been explained. In this work, we show that these fluctuations result 
from resonances of multipole terms which form the scattered radiation and we elucidate the origin of 
each of these secondary maxima. We demonstrate a new theorem which gives the amplitude of 
the multipole modes, and from this, we compute the positions and intensities of the detailed struc- 
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ture. 
what manner it must be modified. 


Introduction. — Soit un faisceau lumineux plan, 
monochromatique, de longueur d’onde À, qui se 
propage dans la direction OZ (fig. 1) et qui est 


Pre te 


polarisé rectilignement, de sorte que le champ élec- 
trique incident soit parallèle à OX. Plaçons au 
point O une sphère de rayon a et d’indice x quel- 
conques : elle perturbe et disperse selon un pro- 
cessus complexe (diffraction, réflexion, réfraction) 
une portion du faisceau incident qui se trouve 
ainsi affaibli. L’étude qui suit se rapporte à la 
détermination de cet affaiblissement ; on ne consi- 
dère pas la répartition de la lumière diffusée et on 
néglige tout eflet quantique. 

L'action de la sphère sur le faisceau incident est 
caractérisée par le « facteur de diffusion » ÆX(«, n), 
c’est le rapport du flux lumineux total diffusé par 
la sphère dans toutes les directions, au flux incident 
sur sa section droite. Cette grandeur présentant un 
grand intérêt pour l’optique des milieux troubles, 
de nombreux auteurs se sont attachés à la déter- 
niner théoriquement. Mie, le premier, en donna 


Finally we show that the present known extinction curve is not correct and we indicate in 


l'expression mathématique rigoureuse, elle est très 
complexe et le calcul numérique en est long et 
délicat. Lorsqu'on ne dispose que de machines à 
calculer classiques, les formules sont presque inuti- 
lisables pour des valeurs de &« — 2x a/} dépassant 
quelques unités. Des méthodes d’approximation 
dues à de nombreux chercheurs ont permis de 
compléter les renseignements fournis par la théorie 
de Mie et de déterminer la forme de la courbe de 
diffusion (fig. 2). 


K\' 


0 5 0 1$ 20 25 30 


FiG. 2, — Courbe de diffusion (valeurs moyennes). 


Une série de calculs entrepris par Penndorf [1}, 
selon la théorie de Mie, à l’aide de calculatrices 
électroniques, montra que la structure indiquée 
ci-dessus n’était qu’approximative. La courbe 
obtenue, publiée pour la première fois en 1953 [2] 
présente une série d’indentations que les calculs 
précédents n’avaient pas montrée. Les caracté- 
ristiques essentielles qui résultent de l'examen de 
cette courbe sont les suivantes (fig 4) : 

1) La valeur asymptotique X — 2 : l'énergie 
interceptée par une grosse sphère est le double de 
celle-qu’elle occulte géométriquement. 
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2) L'existence de grandes fluctuations dont les 
maxima et les minima se placent approxima- 
tivement à des valeurs fixes de 2x (n — 1). 

3) Le fait que la courbe est située en moyenne 
plus haut que son asymptote. 

. 4) L'existence d’une structure détaillée. 

Les trois premiers points ont été examinés par 
Van de Hulst [3] qui leur a trouvé des explications 
satisfaisantes, nous n’y reviendrons pas et nous 
nous attacherons ici à déterminer l’origine de la 
structure détaillée (« Ripple » de Van de Hulst). 


Rappel des résultats essentiels de la théorie de 


Mie. — Avec les seules hypothèses de l’existence 


du champ et de la vérification des équations de 
Maxwell, Mie [4] a pu établir les expressions rigou- 


. reuses du problème de diffraction des ondes électro- 


magnétiques par les sphères. Nous renvoyons pour 


_ le détail de la théorie au travail original de Mie et 


/ 


aux exposés plus récents qui en ont été 
publiés [3, 5]. 
Les résultats généraux sont les suivants : l’onde 


lumineuse incidente excite dans la sphère diffrin- 
gente des vibrations appartenant à une infinité de 
modes multipolaires électriques et magnétiques, la 
somme de leurs rayonnements constitue le champ 
diffracté. L’amplitude et la phase avec laquelle ces 
modes sont excités sont données par l'expression 
des coefficients de diffraction : 


Sm(8) S(oc) ee nS 5m (ct) Sm(B) (1) 
Sm(8) vi nes n D(æ) 5 m(B) 
(x) 

(œ) S 


ana, n) = (—1)n+1? 
. Sync) Ta Sn œ Sm(B (8) 
CALE m(B)" 


Dans ces formules, "= est le rang du multipôle, 
à est la longueur d’onde de la lumière incidente 
dans le milieu extérieur, a est le rayon de la sphère, 

— 27 a/à est le paramètre fondamental définis- 
sant cette dernière, 8 — nx, n est l’indice de réfrac- 
tion de la matière constituant la sphère par rapport 
au milieu extérieur. C’est un nombre réel ou com- 
plexe selon que cette matière est transparente ou 
absorbante. Il est imaginaire pur si la sphère est 
parfaitement réfléchissante ; le milieu extérieur est 
toujours supposé transparent. 

La variable z représentant « ou B, S,(2) et Ch(z) 
sont les fonctions de Riccatti-Bessel définies par : 


. Z Z À 
Sy (z) ee VÉAC) Cmz) = nuits (3) 


OÙ Jyyrp et J_» 17 sont les fonctions de Bessel 
d'ordre m + 1/2. Les dérivées #,(:) et C,(2) 
doivent être prises par rapport à z. On a posé de 


plus : 
Sy (z) Sn 1Cm(z). (4) 


Toutes les grandeurs physiques relatives à la 
diffraction par les sphères sont données par des 
développements en série, où apparaissent les coef- 


bylo, b) RE ie ar és 


(2) 


mto) — Dino 


Ph (z) = 
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ficients de diffraction. Nous renvoyons aux 
ouvrages cités en référence pour leur formulation 
et nous nous bornons à indiquer que l’expression 
du coefficient de diffusion est donnée par : 


[|@ |? + VAE (: 


OX 


Kio à] = 2 $ (0m +4) 
1 


Nouvelle méthode de calcul des amplitudes 
d’excitation des modes multipolaires. — Nous 
allons établir une propriété des coefficients de dif- 
fraction qui nous permettra de simplifier gran- 
dement les calculs. (Considérons les expres- 
sions (1, 2) des coefficients de diffraction et posons : 


À = ShlB) Snlx) — nSmlx) Sm(8) 
TI Sn (B) Co) Fes nCm(x) Sn (8) 


(6) 

l'expression (1) se simplifie en : 

a 

A + iB° 
En séparant les parties réelles et imaginaires, on 

obtient compte tenu de (4) : 


Am(c, n) = CRIE 


(2) 


Gm(&, n) = R[an(æ, n)] 
ë AB La + 
+ ülan(æ, n)] = EE [x + Rp? GE 72 al: (8) 
Remarquons d’autre part que : 
A? B? 
lance, n)|? TS (A? + 1) 
A A? 
(A? ne B?)? ve A2 ee B? 3 H Tam (ot, n)]|. (9) 


La forme analogue du terme b, montre que le 
même résultat est valable également pour ce terme. 
Nous pouvons donc énoncer le théorème suivant : 

« L’'amplitude d’excitation d’un mode multipolaire 
est égale à la valeur absolue de la racine carrée de la 
partie imaginaire du coefficient de diffraction corres- 
pondant. » 

Comme conséquence, l'expression du coefficient 
de diffusion devient : 

9 CO 
=D (2m 1 

pr 


K{«, n) - 1) an(x, n)]| + 11[bulx, n)]l 1. 


(10) 


Les diverses tables donnant directement les 
parties réelles et imaginaires des coeflicients de 
diffraction, on peut juger de l’avantage qu’apporte 
ce nouveau résultat pour le calcul numérique des 
phénomènes. En effet, le théorème que nous venons 
de démontrer permet de déduire à simple lecture 
l'intensité d’excitation des modes, alors que le 
calcul classique aurait exigé la détermination des 
carrés de chaque terme, puis lextraction de la 
racine de leur somme. L'opération eut été longue, 
l'étude d’un seul mode nécessitant le calcul de 
plusieurs centaines de points. 
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Remarque. — Rappelons une propriété intéres- 
sante obtenue par Van de Hulst à l’aide de raison- 
nements analogues. Il établit que le lieu des affixes 
des points représentatifs des coefficients de dif- 


4m (M pair) 
(MX impair) 


bm (M pair) 


Fic. 3. — Lieu des affixes des coefficients de diffraction. 


fraction est un cercle. Les notations employées ici 
étant légèrement différentes, ce résultat se trans- 
forme comme suit : les affixes des coefficients de 
diffraction se répartissent sur deux cercles de 
centres + 0,5 à, de diamètre 1 (fig. 3). Les termes 4, 
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sont situés sur le cercle (1) si » est pair, sur le 
cercle (2) si m est impair. Pour les termes b,, le 
résultat est opposé: cercle (1) pour m impair, 
cercle (2) pour m pair. 


Étude de l’excitation des modes multipolaires 
électriques et magnétiques. — a) TECHNIQUE DES 
CALCULS. — À l’aide de la propriété que nous avons 
établie au paragraphe précédent, nous avons déter- 
miné, dans le cas de z — 1,44, pour chaque valeur 
des coefficients calculés par Penndorf, les inten- 
sités des 34 premiérs termes multipolaires élec- 
triques et magnétiques. Les calculs sont effectués 
par incréments de 0,1 «. Rappelons que seule la 
valeur locale de l'intensité d’excitation est déter- 
minée, mais non le sens de sa variation avec «. 


b) RÉSULTATS DE L'EXAMEN DES COURBES 
D’EXCITATION DES MODES. — Sur les figures, nous 
portons en abscisses, le paramètre «, en ordonnées, 
l'intensité d’excitation des modes. Pour la clarté 
du dessin, nous ne représentons (traits fins) que 
les 10 premiers modes électriques 4, (fig. 4) et 
magnétiques b, (fig. 5). Nous superposons à chacun 
des réseaux la courbe de diffusion X(«, 1,44) (gros 
trait). 


FiG. 4. — Intensités des premiers modes multipolaires électriques (traits fins) et courbe de diffusion Æ{x, n) (gros trait). 
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Fic. 5. — Intensités des premiers modes multipolaires magnétiques (traits fins) et courbe de diffusion X{(«, n) (gros trait). 
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On constate que chaque courbe est constituée 
d’une série de maxima généralement flous ayant 
pour borne supérieure 1, séparés par des minima 
nuls. À ces oscillations amples se superposent des 
ondulations mineures qui nous paraissent difficiles 
à interpréter. Nous ne le tenterons pas, car ceci ne 
présente aucun intérêt pour la suite. Les caractères 
essentiels qui se dégagent de l'examen de ces gra- 
phiques sont les suivants : 

1) Le comportement des modes électriques et 
magnétiques de même rang mn est semblable. 

2) La largeur du premier maximum  décroît 
régulièrement quand m croît et finit, pour les 
grandes valeurs de m, par devenir très inférieure 
à 0,1 «. Il disparait donc de la tabulation, comme 
le montrent les valeurs suivantes relatives au 30° 
mode magnétique 


a = 24,1 Ib») — 0,0065 
24,2 0,0136 
24,3 0,0425 
24,4 0,0292 
24,5 0,0145 


Le premier maximum, égal à {, est situé entre 
24,3 et 24,4, il a échappé à l’analyse. 

3) Au delà de m — 10, le second maximum 
s’affine à son tour. A titre indicatif, nous donnons 
(fig. 6) la courbe d’excitation d’un mode supérieur : 
30€ magnétique. 

4) Les premiers maxima des modes magnétiques 
sont plus fins que les maxima électriques de même 
rang. 

o) Les cinq premiers maxima électriques et 
magnétiques, assez flous, se groupent au voisi- 
nage de &« — 4,5 : il en résulte la première ondu- 
lation principale de la courbe de diffusion. 


6) Pour x > 6 (fig. 5), les indentations secon- 
daires de la courbe de diffusion correspondent aux 
premiers maxima des modes magnétiques des divers 
ordres. 


21 0 25 «em 210 200 À Ga 


F1G. 6. — Intensité du 30° mode magnétique. 


Interprétation des courbes de diffusion. — Pour 
préciser le résultat que nous venons de signaler, 
nous avons tracé, sur un même graphique avec 
les mêmes coordonnées que dans le paragraphe 
précédent, les premiers maxima des modes magné- 
tiques (tirets) et électriques (pointillés), ainsi que, 
superposée à ces réseaux, la courbe de diffusion 
(trait plein). Pour faciliter l'exposé, nous divisons 
en trois régions l'intervalle de variation de «. 

a) 1Te RÉGION à < 18 (fig. 7). — 1) Pour «x > 6, 
à chaque premier maximum des modes magné- 
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Fic. 7. — Section efficace de diffusion (—) et 1°7 maxima des modes multipolaires cos 
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tiques successifs correspond une indentation de 
K(ax, n). 

2) Pour «x > 10, une nouvelle série d’inden- 
tations se rapportant au premier maximum des 
modes électriques des divers ordres, vient s’inter- 
caler entre les précédentes. On observe que l’action 
d’un mode sur la courbe de diffusion est d’autant 
plus marquée que le maximum d’excitation corres- 
pondant est plus aigu : c’est pourquoi les compo- 
santes électriques n’agissent que pour des ordres 
assez élevés. 

3) Pour « < 6, l’origine des indentations semble 
difficile à élucider. Tenant compte de l’orthogo- 
nalité des modes et de leur action purement addi- 
tive, nous proposons l'interprétation suivante : 


5€ indentation : (x — 5,3), superposition des 
4e modes électrique et 5€ magnétique. 
4e indentation : (x — 4,5), superposition des 


3e modes électrique et 4° magnétique. 

Les trois précédentes, réduites à des inflexions, 
résultent de la superposition des trois premiers 
modes des deux genres. [l ne semble pas possible de 
préciser davantage. 


b) 22 RÉGION 18 < à << 23 (fig. 8). — Les inden- 
tations de X{(x, n) correspondent encore aux pre- 
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Fic. 8. — Section efficace de diffusion (—) et 127 maxima 
Magnétiques — — —, 


des modes multipolaires | . 
1 Électriques 


miers maxima des modes successifs, mais leurs 
amplitudes sont irrégulières, certaines ont même 
disparu. Il ne s’agit là que d’une apparence, due à 
un intervalle trop large des points calculés. En 
effet, si l’on porte sur la courbe de diffusion les 
points tabulés, on observe que, selon qu'ils coïn- 
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cident ou non les avec maxima d’excitation, l’ondu- 
lation est forte ou faible. Remarquons d’autre part 
qu'une des indentations (x — 21,4) correspond, 
non plus au premier, mais au second maximum 
d’excitation du 22° mode magnétique. 


c) 32 RÉGION 23 << x << 30 (fig. 9). — Les inden- 
tations sont faibles et coïncident cette fois avec les 


1° max du 81° mode électrique 
K fang des modes 
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Fi. 9. — Section efficace de diffusion (—-) et 2° maxima 
: : : Magnétiques — — —, 
des modes multipolaires Électriques ........ 


seconds maxima des différents modes. Celles qui 
sont dues aux premiers maxima d’excitation ont 
disparu, à l’exception du 31€ électrique. Nous 
admettons, comme pour la région précédente, que 
cette anomalie est due à ce que les pointes des 
premiers maxima sont si fines qu’elles échappent à 
l'analyse. 

Pour corroborer notre hypothèse, nous avons 
indiqué, dans le tableau ci-dessous, la plus grande 
valeur numérique d'amplitude d’excitation tabulée 
au voisinage des premiers maxima contenus dans 
cette région. 


Mons AMPLI- Move AMPLI- 
ÉLECTRIQUE Œm  TUDE MAGNÉTIQUE Xm TUDE 
28 23,3 ° 0,1700 28 22,8  0,0803 
29 24,0 0,2748 29 23,6 0,1503 
30 24,8 0,0774 30 24,3  0,0425 
51 25,5: 0,9289 31 25,1  0,0444 
32 26,2 0,0412 32 25,8 0,4309 
33 27,0 0,0794 33 26,5  0,0209 
34 27,7  0,0478 34 27,3  0,0305 
35 28,5 0,0200 35 28,0  0,0389 
36 29,2 0,1463 36 28,8  0,0080 
37 29,9 0,0136 37 29,5  0,0168 
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Positions des premiers maxima des modes multi- 
polaires. — Il serait trop long de déterminer direc- 
tement les positions des premiers maxima des 
modes multipolaires à partir des équations (1) et (2) 
qui sont très complexes. Seul, un examen attentif 
des amplitudes tabulées permet de les localiser 
à 0,1 « près. 

Pour préciser comment varie leur position en 
fonction de x, nous portons, en abscisses (fig. 10), 


D 


Seconds max. des 
modes électriques 


max. des 


modes magnétiques 


Premiers max. des 
modes magnétiques 


Premiers max. des 
modes électriques 


Numero d'ordre des modes multipolaires 


5 10 15 20 25 30 x 


Fic. 10. — Position des premiers et seconds maxima 
des modes multipolaires électriques et magnétiques. 


les valeurs de « se rapportant successivement aux 
deux premiers maxima des modes électriques et 
magnétiques, en ordonnées, les rangs des multi- 
pôles correspondants. À l'exception des tous 
premiers, les points ainsi obtenus se groupent sur 
quatre droites sensiblement parallèles. En d’autres 
termes, les maxima de même origine sont à peu 
près équidistants. 

On peut représenter les valeurs de « corres- 
pondant à ces maxima par les formules approchées 
suivantes : ë 


1e7 max. électrique 


œ 0,764 m + 1,86 
2e max. électrique GX) 

œ 

œ 


,276 m + 4,76 
0,756 m + 1,65 
0,765 m + 4,67. 


1eT max. magnétique 
2e max. magnétique 


Dans tout l'intervalle 10 < x < 30, l’écart entre 
les valeurs déduites des tables et celles données par 
ces formules est inférieur à une unité du premier 
ordre décimal. Il croît légèrement en dehors de cet 
intervalle. Son sens prouve que les courbes pos- 
sèdent une légère concavité tournée vers le bas, que 
l’on peut négliger en première approximation. 


HU 


Détermination quantitative de l’amplitude des 
indentations. — a) LE MAXIMUM D'EXCITATION 
PEUT ÊTRE CONSIDÉRÉ COMME INFINIMENT ÉTROIT. 
— On peut déterminer l'amplitude des indentations 
par le raisonnement suivant : dans l'expression (10) 
qui représente le coefficient de diffusion, tous les 
termes sont positifs et additifs. D’autre part, nous 


avons montré plus haut que les premiers maxima 
des ordres multipolaires supérieurs sont constitués 
par une pointe très aiguë et isolée, on peut donc 
séparer dans la série le terme qui atteint la valeur 1. 
Il n’agit que localement puisqu'il devient négli- 
geable dès que l’on s’écarte de la valeur de « corres- 
pondant à son maximum. En conséquence, l’en- 
semble des autres termes donne la courbe moyenne 
et le terme résonnant l’indentation. 

Ce raisonnement donne une approximation cor- 
recte, dès que la pointe est assez aiguë (grandes 
valeurs de x). La contribution du terme de rang m 
est alors : 


A À (2m +4). (11) 


Le premier maximum du terme 43, 
(a )1=0;929 


coïncide presque avec l’un des points calculés et 
permet une vérification de ce résultat. La hauteur 
de l’indentation qui en résulte est : 


AK calculé — 0,194 
AK mesuré — 0,18. 


Si l’on extrapôle la loi linéaire donnant les posi- 
tions des maxima, on constate que les amplitudes 
diminuent selon une loi en 1/m lorsque m croît 
indéfiniment. 


b) LE MAXIMUM D’EXCITATION A UNE LARGEUR 
NON NÉGLIGEABLE. — Dans ce cas, l’amplitude de 
l’indentation dépend del’acuité du maximum d’exci- 
tation du multipôle correspondant. Nous n’avons 
pas effectué les calculs qui s’y rapportent parce 
que, dans l’intervalle où ils s'imposent (fig. 8) les 
indentations sont assez larges pour être bien 
définies par les tables. 


Application : Forme correcte de la courbe de dif- 
fusion dans l’intervalle 23 < «x < 30. — Nous 
avons noté plus haut que, dans cet intervalle les 
indentations dues aux premiers maxima des modes 
avaient disparu à la tabulation. Nous avons cal- 
culé, selon notre méthode, leurs amplitudes et nous 
les avons portées sur la courbe moyenne de dif- 
fusion (fig. 11). 

La structure secondaire désordonnée qui avait 
étonné différents auteurs (3) s’est régularisée. La 
succession des indentations simples et dédoublées 
de la courbe de Penndorf a disparu. 

La nouvelle structure est formée par la super- 
position de quatre séries d’indentations corres- 
pondant aux deux premiers maxima des modes 
électriques et magnétiques. L’évolution de chaque 
série est la suivante : les ondulations d’abord larges 
et peu intenses s’affinent en même temps que leur 
amplitude croit. Cette dernière passe par un maxi- 
mum puis décroît lentement tandis que la largeur 
de l’indentation diminue constamment. La pre- 
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mière série qui apparaît est celle des 1275 maxima 
magnétiques, ‘puis, celle des 17 maxima élec- 
triques vient s’y intercaler. Pour une valeur 
de « > 22, apparaissent successivement les 2€5 ma- 
xima magnétiques et électriques. Comme l’inter- 


2,5 
F2 
20 
13 2 25 TR | 25 EE) C7 
Fi1c. 11. — Structure détaillée de la courbe de diffusion. 
(1) D’après les tables de Penndorf — courbe très 
altérée par une analyse numérique insuffisante. 
(2) Courbe correcte. — Détails restitués par notre 
méthode. 


valle entre indentations est le même dans toutes 
les séries, on obtient l'aspect de la figure 11. Si l’on 
extrapôle au delà de « — 30, le même processus 
doit continuer avec apparition de nouvelles inden- 
tations correspondant aux 3%$ puis aux 721em€s ma- 
xima des deux ordres. La structure ultime est 
étrangement complexe. : 
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Signalons que sans l'interprétation nouvelle que 
nous avons donnée, il eut été très difficile d'obtenir 
ce résultat. Il eut fallu effectuer les calculs numé- 
riques à intervalles très serrés, le travail eut été 
considérable et surtout l’extrapolation impossible. 


Conclusion. — La structure détaillée de la courbe 
de diffusion de la lumière par les particules sphé- 
riques transparentes est due aux premiers maxima 
des modes multipolaires successifs. L'intensité des 
indentations caractéristiques de cette structure 
peut se déterminer en séparant l’action du terme 
multipolaire qui passe par un maximum de celle de 
l'ensemble de tous les autres termes. 

Nous avons, par une méthode nouvelle, recons- 
titué l’aspect de la structure détaillée de la courbe 
de diffusion correspondant aux valeurs de « infé- 
rieures à 30. Elle est constituée d’une succession 
complexe d’indentations régulièrement disposées 
dont certaines sont très aiguës, leurs positions sont 
données par des lois sensiblement linéaires. Nous 
montrons comment l’extrapolation de ces résultats 
permet de prévoir la structure de la courbe pour les 
grosses sphères. 

Nous tenons à exprimer notre gratitude à M. le 
Pr Grivet pour l'intérêt qu’il a porté à cette étude. 
Nousremercions M. Dews, du Geophysical Research 
Directorate de Cambridge (U. S. A.), qui a rendu ce 
travail possible en nous procurant gracieusement 
un exemplaire des tables de Penndorf. 


Manuscrit reçu le 4 avril 1958. 
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PROPAGATION GUIDÉE DES ONDES ÉLECTROMAGNÉTIQUES 
DANS UN PLASMA ÉLECTRONIQUE GYROÉLECTRIQUE 


Par B. UNAL et Tnéo KAHAN, 


Institut Henri-Poincaré, Paris. 


Résumé. — On établit l'équation aux dérivées partielles du quatrième ordre régissant la pro- 
pagation d’ondes électromagnétiques dans un guide d’onde rempli d’un magnéto-plasma. 


Abstract. — À partial differentiel equation of fourth order is derived which describes the elec- 
tromagnetic wave propagation in a wave guide filled with a magneto-plasma. 


1. Introduction. — Nous nous proposons 
d'étudier la propagation des ondes électroma- 
gnétiques dans un milieu ionisé gyroélectrique 
doué d’une perméabilité magnétique y scalaire et 
d’untenseur diélectrique ((£<)) caractérisant lanature 
anisotropique du plasma gyroélectrique [1]. L’ani- 
sotropie est créée par un champ magnétique 
constant À, qui sera pris parallèle à oz, et les 
éléments du tenseur ((e)) et uw sont fonction de ce 
dernier [2]. D’autre part, (€) et seront en 
général fonctions de la fréquence « de l’onde élec- 
tromagnétique variant suivant la loi exp (— iot). 

Le déplacement électrique D est lié au champ 
électrique E par la relation 


D — ((<)) E, 


es 1x in 
— 1x € re | 
— in —1Ù 3 ! 


B—uH 


D'autre part 


Lorsque À, est homogène, les composantes du 
tenseur ((e)) sont des constantes. 


2. Équations de base. — La propagation de nos 
ondes est régie par les équations de Maxwell : 
— i((e)) k E = v AH, V.E=0 
ik H=VAE, V.B— 0. 


I 


(1) 


Tirons E de la première équation et portons-le 
dans la deuxième 


£ al 
E—=(()-VAH 


0 


et 
uk H—=AA((e) 1! V A H] 
V'H—0 - B) 
où 
((e)}= 1 = ((x)) 
dé 0 de M, K +ik N+in 
td ne AU R —"# M; L + il 
de Ag. de RE LL 11 M 


avec 
1 2 
dx = M; al Es ch 
: 1 : 
ox —= E + ik im — LE K), 
2 1 , 
ax = N + in = (x — ie n), 
æ 7 Y 1 2 
Ey — x» Xy M; PAS n°), 
À 
ay = L + il A 1e;iC), 
œ Y 2 1 2 
GNT Pr NY RUES Ms = (e E2 »°), 
d = EE, s9es — e C2 — €, n° — €, x?. 


Le système (2) constitue les équations de base du 
champ pour déterminer H. 
Pour les mettre sous une forme convenable, 


. écrivons-les dans un repère cartésien p, 4, v; la 


projection sur l’axe p(x, y, z) donne : 
uk H, = dul((e)-V À H}, — D,[((e)) VA H}u 
(p, &, y) en permutation circulaire 
= dulei(V À H)i] — doi (V À Hi]. 
(i— x,y,z;sommationsur lesindicesrépétés sous-entendue). 
En portant 
(V AH); = 0; Hx — d H; 
dans cette expression, il vient pour les compo- 
santes de 
{4 ki H, FA D H;, 
(à == (a, dy nn a du) dj FE (oi, dy me x du) dk (3) 
C’est un opérateur différentiel du second ordre. 


D'une façon explicite le système (3) a pour 


expression 
(D5 — uk) Hy + Où Hy + De H, = 0, 


DH, + (D — uk) H, + DH, = 0, 


I 
© 


D Hs + D Hy + (D — y ki) He 
de Hz + dy Hy GTR d: Hz — 0, 
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Éliminons A, entre la {re, la 2e et la 4 équation ; 
nous obtenons ainsi 


DORE RE 4 PR 


f étant l'opérateur différentiel du 6€ ordre sui- 
vant : 


Li = À [DE Di — De D] de + [DE Dé — Dé Di] D 
+ [Dé Dé — De Di] dy + Le? h0 de — 1 A6 De + De] D 
+ de KGTDS De + O9 de] À De 


En utilisant (3), on obtient pour #; l'expression 
suivante : 


Mu HUM vs + Mid + M 
+ (M, LE M:) 0 — 2% de dy — 21 dy d4 — 2N dz del 


+ le 08 + 008 + € 08] v} D = LD 


avec : 


Vi di EM DE Ve Lio, 
VAR er 


Pour les autres composantes, il vient par per- 
mutation circulaire : 


Er Hy= La H,= 0, avec: Fr = 1, 


Por He = Pin Hz = 0, avec Fm = LH. 


Comme 


DO AN FA 


si 
DS HE NO, d Hi 0, 
alors 1, sera solution de : 
LHE=0;: 


De même, avec les conditions. : 


de My 0, dHy Z0 


et 
dy H, Æ 0, de H, # 0 
LHy = 0, 
LH, = 0. 


Nos équations de base prendront donc la forme 
finale : 


LH 
VE 0: (a) 


I 
© 
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Ce nouveau système contient le même nombre 
d'équations que celui de départ. Cependant les 
deux systèmes ne sont pas équivalents. Le sys- 
tème (1) est du second ordre, tandis que (4) est du 
quatrième ordre. Le dernier admet en plus des 
solutions dépourvues de signification physique. 
Afin d'éliminer ces solutions physiquement inac- 
ceptables on ajoute une des équations du système 
(4) comme condition subsidiaire. 

Prenons la projection sur Oz de (1) ; il vient 


de! lu A + M3 08 HMS, + 2Ù)] He — (A — ik) 02 Hi} 
+ Du fu ke + M0 + Mi + d4)] Au 


— (K + ik) D He} +R = 0. 
où 


R = dAl(L + üÙ) dx — (N + in) dy] (de Hy — dy Hx) 
— 2K da dy(dr He + dy Hy). 
3. Première approximation. — a) Prenons » et € 
suffisamment petits et tels que 
n #0, n=n—0; 
bYEi te: 


Dans ces conditions l’on a : 


M; — M, — Mn d = cg(e — x?), 
M? KA 1 KO 
e ei — x? mé 
MIRE 


L'= {uk +uKUM + M:) vs + 2M%] + MM;N» 
+ (MM + M3) — 41] vè à 
+ (M? Ie) D — 29, (Lo No 1. 
L'= (y? + Ki) (V8 + KE) — 20UN 0e + Ld): 


Posons d, — 1y; la coordonnée z étant cyclique, 
la condition subsidiaire devient : 


Da À eu 6 + MY?) Ho — ik y? Hu} 
+ ouf lu 76 + My?) Hy + ik Y? Ha} 
+ ÉvI(L € ü1) de — (N + in) d9] (De Hy — dy Ha) =: 


L'étude analytique de ces équations aux dérivées 
partielles fera l’objet d’une publication ultérieure: 


Manuscrit reçu le 24 janvier 1958. 
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DISSOCIATION ET IONISATION DE L'AIR PAR UNE ONDE DE CHOC 


Par J. THOUVENIN, 


Agrégé de l’Université, Ingénieur au C. E. A. 


Résumé. — La composition de l’air et la valeur de son énergie interne sont calculées pour des 
températures allant de 3 500 à 11 500 °K et des densités de 4 à 12 fois supérieures à la densité 
normale. D’autre part, on détermine, en fonction des mêmes paramètres, l'augmentation d’énergie 
interne de l’air due au passage d’une onde de choc. En égalant les deux expressions de l’énergie, 
on obtient une relation entre la température T'et le rapport des volumes moléculaires V,/V devant 
et derrière le front de choc. Les valeurs des autres grandeurs physiques, pression, vitesse du front, 
vitesse d'écoulement, taux d’ionisation dans le front, sont alors déterminées si l’on connaît l’un 
des paramètres T ou V,/V. Inversement, la mesure de l’une quelconque de ces grandeurs phy- 
siuqes permet de fixer la valeur de toutes les autres. 

Le calcul montre que pour des chocs intenses (vitesses supérieures à 7 000 m/s), la dissociation 
de l’oxygène est totale, celle de l’azote avancée à plus de 50 %, et la proportion d’électrons libres 
supérieure à 0,1 %. 


Abstract. — The Composition and internalenergyof airarecomputed for a range of temperature 
from 3 500 up to 11 500 °0K and a range of density from 4 up to 12 times normal density. In 
another connection, the increase of internal energy of air by a shock wave traveling through it is 
evaluated in terms of the same parameters. By adjusting both expressions of energy, a relation 
between the temperature 7 and the ratio of molecular volumes V,/V ahead and behind the shock 
front is obtained. The others physical variables, pressure, front velocity, material velocity, 
degree of ionization, can be then computed if either of parameters T or V,/V is known. Conver- 
sely, measurement of any one physical variable makes it possible to get values of all the others. 

Present calculations show the oxygen to be completely dissociated by strong shock waves (velo- 
cities above 7 000 m/s), the nitrogen by a rate of 50 % higher, and the concentration of free 
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electrons to be over 0,1 %. 


1. Introduction . — Les études modernes sur les 
propriétés des gaz traversés par une onde de choc 
intense tiennent compte de l’effet important que 
jouent, dans l’évaluation de l’énergie interne, la 
dissociation (s’il s’agit d’un gaz polyatomique) et 
l’ionisation. L'énergie délivrée au passage du front 
est en effet assez grande pour exciter ces degrés de 
liberté « lointains » : de fait, une onde de choc de 
7 000 m/s traversant l’air ambiant laisse l'oxygène 
‘entièrement dissocié, l’azote à moitié dissocié, et 
une proportion d'électrons voisine de 0,1 %. Les 
phénomènes d’ionisation ont une évidence expéri- 
mentale directe dans la luminosité intense qui 
accompagne le passage de l’onde de choc. 

L'effet de la dissociation et de l’ionisation sur le 
calcul des relations de choc a d’abord fait l’objet de 
travaux de Dôring et Bürckhardt [1], [2] sur l’air. 
Toutefois, il semble à la lecture des documents 
cités, qu'une erreur systématique se soit intro- 
duite dans leur calcul de la fonction de partition 
de vibration, erreur qui se traduit par des valeurs 
légèrement trop grandes de l’énergie interne (d’une 
dizaine de kilojoules par mol.g). En outre, un 
travail expérimental récent de Christian et 
Yarger [3] a montré que l’énergie de dissociation 
de l’azote effective dans une onde de choc était 
9,764 eV.et non 7,373 eV : Bückhardt avait utilisé 
la seconde de ces valeurs, toutes les deux étant 
également admises par les spectroscopistes. 


Signalons aussi la Thèse de Bond [4] qui, tra- 
vaillant sur l’argon et s’affranchissant ainsi des 
difficultés inhérentes aux molécules diatomiques, 
a pu, grâce à une analyse cinétique, déterminer les 
diverses phases de l’ionisation et l’étalement du 
front de choc qui en résulte. 

Nous reprenons ici, pour le! cas de l'air, la 
méthode de calcul proposée par Bürckhardt, en 
l’affectant toutefois des corrections signalées plus 
haut. Les questions s’introduisent dans l’ordre 
sumwant : 

1° Quelle’ est, à une température 7 donnée, et à 
une compression donnée (définie ici par le rapport 
des volumes moléculaires à la température ordi- 
naire 7, et à la température 7”, soit V,/V), la 
compression de 1 mol.g d’air ? 

20 Quelle est son énergie interne, en fonction 
de T'et de V,/V ? 

30 Connaissant celle-ci, quelle est la température 
effectivement atteinte, et quelle est la compression 
obténue derrière une onde de choc de vitesse 
donnée D ? 


2. Constantes d’équilibre des réactions de disso- 
ciation et d’ionisation de L’air. — a) Considérons 
une réaction gazeuse à volume constant, suscep- 
tible d'atteindre un équilibre à la température T: 


ABZZ A +B, (1) 
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On définit le taux de réaction x, à un stade donné 
de l’évolution, par le nombre de particules A (ou B) 
actuellement présentes rapporté au nombre ini- 
tial NV, de particules AB contenues dans le 
volume V : 


Ne 
Ne MOINS 

NAB 

A NE 


sOIt: 


Na = (1 — 2%) No. @l 


La valeur de x à l’équilibre est celle qui rend 
minimum le potentiel thermodynamique qui gou- 
verne la réaction. Puisqu’il s’agit ici d’une réaction 
à volume constant, ce potentiel est l'énergie libre F. 
Évaluons F en fonction de x. 

L'expression statistique de l'énergie libre F, 
de NV, particules d’espèces p, occupant le volume V 
à la température T, est : 


Fo = No No kT[t + Log |, (3) 
No 
ous: 
=) — énergie du niveau fondamental de la particule p. 


Z» = fonction de partition de la particule p rapportée 
à son niveau fondamental. 


On remarque que l’expression de F, ne contient 
aucune constante arbitraire. En effet, la Mécanique 
quantique permet de calculer la valeur de e,, et 
l’on sait que pour une molécule diatomique en 
particulier, cette valeur est différente de zéro. 

Appliquons l'expression (3) aux particules À, B 
et AB qui constituent le système en évolution. En 
remplaçant les N, par leurs valeurs (2) on obtient 
pour : 


= Fa + Fe + Fa 
une fonction de x, dont on annule la dérivée. On 


obtient ainsi l'expression classique donnant le taux 
d'équilibre : 


ETEpTE 
2 T 5 A B ‘AB 
É 1e = VAN VAT: CARTE TONER 
LT N, Las 
En effet: W — ea + ep — esp représente 


l'énergie absorbée dans une réaction fictive où les 
particules se trouveraient toutes dans leur niveau 
fondamental : c’est cette quantité que l’on a 
convenu d'appeler l'énergie de réaction. Dans notre 
notation, cette énergie est positive si elle est 
absorbée dans le sens 


AB A+B—W 
W = Es + Ep — €AB 2 0. (4) 
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La constante d'équilibre à volume constant de la 
réaction (1) s’écrit donc : 


a? He Za ZB eo 


C ) 
No ZaB : (5) 


1— 2 


Cette expression est formellement identique à 
celle donnée par Bürckhardt [1]. Mais en la resti- 
tuant, nous avons précisé le sens et l’origine des 
diverses grandeurs : 


W — énergie absorbée dans la réaction standard 
(au zéro absolu), c’est-à-dire énergie de 
dissociation ou d’ionisation. 


Z, = fonction de partition du constituant p 
rapportée à son niveau fondamental, 
c’est-à-dire : 


Zn = E g» e KT 16) 
n 


où «, représente l'énergie du niveau fondamental 
de la particule p et e, l’énergie du niveau »# de 
poids statistique g,. Or, dans l'expression de Z 
pour les molécules diatomiques (0,, N,) qui com- 
posent l'air, Bürckhardt n’a pas fait intervenir 
(e, — &) mais «,, de sorte que les constantes C des 
réactions de dissociation qu'il a calculées sont trop 
grandes du terme exp (e,/kT). 


b) CALCUL DE LA FONCTION DE PARTITION D'UNE 
MOLÉCULE DIATOMIQUE. — Il est classique de 
séparer de l'énergie totale de la particule l’énergie 
de translation de son centre de gravité, ce mouve- 
ment de translation étant indépendant des mou- 
vements internes. On écrit donc : 


Li LENS Tr (7) 


La fonction de partition de translation d’une 
particule de masse m se déplaçant dans le volume V 
est donnée par l’expression : 


(8) 


2x mkT\°12 
Ztr. => ( L ] 


h? 


La fonction de partition interne Zf est relative 
aux mouvements internes de la molécule : rotation 
autour d’un axe normal à la ligne des centres ; 
vibration des atomes le long de la ligne des centres ; 
excitation des électrons. 

Ces mouvements ne sont pas rigoureusement 
indépendants. Ainsi, l’on admet que le modèle de 
vibration de la molécule est celui d’un oscillateur 
harmonique. Toutefois, si l’ensemble est animé 
d’un mouvement de rotation autour d’un axe 
perpendiculaire, la force centrifuge tend à écarter 
les atomes, et cette force dépend de l’écartement. 
Les oscillations ne peuvent donc être considérées 
comme harmoniques. De même, le moment 
d'inertie de la molécule varie en fonction de l’écar- 
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tement des atomes, et l’ensemble ne peut rigou- 
reusement être traité comme un rotateur rigide. Le 
couplage de ces deux phénomènes est décrit par un 
modèle plus précis de la molécule : celui du rota- 
teur vibrant, 

Ce couplage peut être négligé aux températures 
ordinaires (jusque vers 1 000 °K pour O, et N;). 
Dans l’application aux ondes de choc, où des 
températures de 6 000 à 10 000 CK sont atteintes, 
il est clair que nous devons en tenir compte. Un 
calcul rigoureux de Z° est malheureusement quasi- 
impossible à conduire entièrement, et une méthode 
d’approximation doit être employée. 

Suivant une évolution sommaire, la correction 

due au couplage rotation-vibration porterait sur 
une majoration de Z? de l’ordre de 10 % 
10 000 °K. Nous préférons cependant réserver à 
une prochaine note l’exposé d’une méthode précise 
de correction, et donnons ici les résultats complets 
pour le cas idéal où il n’existe aucun couplage. 
Enfin, les fréquences de rotation et de vibration 
dépendent légèrement de l’état d’excitation des 
électrons. Nous en avons tenu compte. 

Ces précautions prises, les valeurs des niveaux 
d'énergie pour les trois mouvements, supposés indé- 
pendants, de la molécule sont : 


.oJ(oJ + 1) (9) 


Erotation — XCBn 


By = constante rotationnelle (em—t) pour l’état d’exci- 
tation n. 


J — nombre quantique rotationnel = 0, 1,2, 


o — facteur de symétrie — 2 dans le cas des molécules 
diatomiques symétriques. 


Evibration — Àc w,(v + 1/2) (10) 
où : 
©, — fréquence fondamentale de vibration (em—1) pour 
l’état d’excitation n. 


v — nombre quantique vibrationnel = 0,1,2, ... 
£excitation — AC Vi (11) 
Où : 
YA = fréquence émise lors du saut d’un électron de l'orbite 
excitée sur l'orbite de l’état fondamental. 


Les poids statistiques des différents niveaux sont 
respectivement : 


Srotatioo —= 26 J +1 (12) 
£vibration — 1 

Lexcitation — £n — dégénérescence de l’état excité n. 

Le niveau fondamental (J — 0, v —0, n —0) 


a pour énergie : 
Eo = 1/2 ko. 


DISSOCIATION ET IONISATION DE L’AIR 641 


La Mécanique quantique prévoit qu’au zéro 
absolu les molécules diatomiques ont une énergie 
non nulle, et que cette énergie est due à un état ‘de 
vibration d'amplitude finie. 

Compte tenu de l’expression (6), Z' s'écrit : 


à __ int 0 
Zu = ZE gm.e 7 (14) 
J,v.n 
avec : 
Eint. — €o = hCBy.0J(o] + 1) + he w,.0 + hey 


gint. = (26 J + 1) g 


Les sommations sur J et v se séparent. Il n’en 
serait pas ainsi si le modèle plus rigoureux du 
rotateur vibrant avait été adopté, et l’on voit 
combien les calculs seraient alors compliqués. La 
sommation sur J peut être remplacée par une inté- 
grale pour des températures supérieures à 20 °K : 


(Ce) LR SnorE D) KT 
[ (207 + 1)e ‘FT dr. (5) 
La sommation sur v est simple : 
ss EEE ï : 
> AC a (16) 
F 12-SeurAT 


Finalement, l’expression de Zf rapportée à son 
niveau fondamental est : 


AT D gn € 
che Br hco,, 
AS 6 NET 


Zi = (17) 


Cette expression diffère de celle donnée par 
Bürckhardt, et que nous rappelons ci-dessous : 


Zi = Fe LEE _ 


EL) 
a 


rm En 9 re On 
KT 
2 ho, 
2° kT 


par des termes e qui traduisent le fait 
que les énergies ne sont pas rapportées au niveau 
fondamental, mais à l’énergie zéro, ce qui est 
incompatible avec la définition (4) de l’énergie de 
réaction W. 


€) CALGUL DE LA FONCTION DE PARTITION D'UN 
ATOME. — Le seul degré de liberté interne absor- 
bant de l’énergie étant l'excitation des électrons, 
l'expression (6) se réduit à : 


(18) 
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avec 
Roy, 
VA — Y En € xT (19) 
12 
et : m — masse de l’atome A. 


d) EXPRESSION DE LA CONSTANTE D’ÉQUILIBRE. 
— Revenons à la réaction (1). La particule qui 
subit la réaction, AB, peut être une molécule ou un 
atome. En effet, nous serons conduits à envisager 
l’ionisation directe de la molécule, par une réaction 
du type : 

OO 


Dans ce cas, la particule À est aussi une molé- 
cule, la particule B un électron. Rappelons que, 
pour l’électron : 


Zit=9, (20) 


Pour conserver la généralité, nous écrirons donc 
l'expression de la constante d'équilibre : 


[A] [B] pe cie AA 25 QE (4) 


h° NS. 


ZaB 


__ MmA.mB 
H PRIE 


MAB = MA + MB (22) 


et où les Z* sont données par les formules (17), 
(19) ou (20). 

.. Deux paramètres indépendants déterminent la 
valeur de C': ce sont la température T et le 
volume V. Il est commode de les séparer et de 
poser : 


CE ATIOTITS (23) 
OULE 
K(T) = KT\92 Vo ZA ZE (24) 
% h° } No FANS k 


Désormais, nous prenons pour À, le nombre 
d’Avogadro et pour V, le volume moléculaire 


No 


normal (22 415 cm3). Le rapport V,/V représente 
ainsi la Compression, quantité sans dimension. 


3. Composition de l’air à une température 7 
pour une compression V,/V.— a) Les six réactions 
du type (1) que nous considérons sont les suivantes : 


0e -0 OL 0 eV 
OHETOS - HMS ATEN 
DR ET OF Eee Tee" (25) 


Ne N4N — 9,764 eV 
No ENT en eV 
NE Z NT Le 2445 eV. 


Il faudrait sans doute tenir compte de la for- 
mation d’oxydes de Pesotss en particulier NO : 


N+ 0 NO 


Toutefois, les calculs de Bürckhardt ont montré 
que la proportion de NO était négligeable à partir 
de 7 500 ©K, passant par un maximum de l’ordre 
de 7 % vers 4 500 0K. Comme les températures 
d'ondes de choc qui nous intéressent sont supé- 
rieures à 7 500 °K, nous n’avons pas fait intervenir 
la formation de NO dans nos calculs, quitte à 
commettre une légère erreur dans le cas des ondes 
de choc faibles. 

Les valeurs des énergies de Ascer et d’ioni- 
sation sont celles données par Herzberg [| Spectra of 
Diatomic Molecules (1955) — Atomic Spectra 
(1944)], sauf l’énergie de dissociation de l’azote 
que nous avons prise égale à la valeur indiquée 
par Christian et Yarger [3]. 

Les constantes X(T) pour ces 6 réactions ont été 
calculées à partir des données spectroscopiques 
relevées dans Herzberg (op. cit.) pour les molécules 
et dans Bacher et Goudsmit | Atomic Energy States 
(1932)] pour les atomes. 

Les valeurs des six X(7) sont portées dans le 
tableau I pour neuf valeurs de la température 


variant de 1 000 en 1 000 ©K, de 3 500 °K à 


11 500 K. 


TABLEAU I 


CONSTANTES D'ÉQUILIBRE K;(T') 


K; K> K3 K4 K; Ke 
POK * OL O+O0: Om OF Len DE TOP He CNE N EN NS NTE Lo 
3 500 0,02329 0,514 X 10—17 4106 X 10—19 3,774 xX 10% 2,938 x 10—21 9,644 X 1020 
4 500 0,9025 1,077 X 10—12 4,357 X 10—14 4,916 X 10— 4,168 xX 10—18 6,314 X 1018 
5 500 8,988 4,279 X 10—10.. 1,082 X 10—11 4,712 x 10—4 8,422 X 10—13 7,755 X 10-12 
6 500 43,34 2,777 X 40—8 1,151 X 10—? 0,01123 1,709 X 10—19 1,105 X 10— 
7 500 199 60528402 3,640 X1407 0,1169 8:688:X 4107 4,282 RARE 
8 500 322,3 6,500 x 10—5 5,244 X 10—7 0,7152 1,797 X 10—7 7,115 X 107 
9 500 633,7 R 29510 L'H07 ATOS 3,044 20020 6,604 x 40m 
10 500 1089 2,004 X 10—4 2,519 X 10— 9,977 1,430 X 10—5 4,033 X 10—5 
11 500 1694 7,224 X 10—1 1,082 X 10—4 269,0 7,339 X 10—5 1,808 x 10—4 
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Les compressions réalisées derrière le front d’une 
onde de choc pouvant aller jusqu’à 12, nous avons 
utilisé pour le calcul de € = K(T). VV, les sept 
valeurs suivantes : 


PL ST 0 96 10 411 = 12; 


> 


b) Pour un couple donné (7, V,/V), nous dis- 
posons de six équations : 


[Oj® OT ter =. “[0#] Le) 
D eo fe op … 
[Nje ENS] el EN] el 
2 De à IN) 


avec Cr AT) VTIV: 
Trois équations supplémentaires sont fournies 
par la conservation de la masse et celle de la charge : 


D lO) AO 
[O2] NO re 9. dre 1 
Ian (27) 
NAN + ER ed 
[O1 + [NX] + 104]+ [NY] = fe] : 


x. et «, désignent les proportions normales d’oxy- 

gène et d'azote dans l’air : &, —0,212 ; «2 = 0,788 

Rappelons que les quantités entre crochets repré- 

sentent les concentrations moléculaires, c’est-à- 
dire : _ 

p 

IP] N, 


N, étant le nombre de particules d'espèce p dans 
1 mol.g d’air et V, le nombre d’Avogadro. 

Les neuf équations (26) et (27) permettent de 
calculer les neuf concentrations [p]. Le système 
est simple à résoudre si l’on néglige l’ionisation 
devant la dissociation, ce qui est parfaitement légi- 
time dans la gamme de températures de choc 
choisie (le taux d’ionisation sera au maximum 
de 4 % lorsque les dissociation seront complètes). 
Dans ce cas, on calcule d’abord [O,] et [O] au 
moyen de : 


et de même [N,]et [N]. Puis : 
Le] = VGz [02] + Cs [OT C5 [No] + Ce IN] 
et enfin : 


mt [0] an [ 2] 
[0*] Er Ce [e] , [NS] ST EM 

= 0 a ENT 
IQ Cao > [Nt] = Ce fo: 


Introduisons deux nouvelles quantités : 
1%e— [0,1 + [0]'+ [0] + [O*+] + [N:] 
+ IN] + [NX] + [NY] +{e] 
(ele à nombre d’électrons 
1 +< nombre total de particules 


(28) 


(29) 


Lo — 
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£ est l'augmentation relative du nombre de parti- 
cules due à la dissociation et à l’ionisation. x, est 
le taux d’ionisation. Ces onze grandeurs, qui carac- 
térisent la composition de l’air portée à T 0K sous 
une compression V,/V, ont été calculées pour tous 
les couples (7, V,/V) définis à la fin du para- 
graphe 34. 


4. Énergie interne de 1 mol.g d’air en fonction 
de T'et de V,/V. — L'énergie interne de W, parti- 
cules d’espèce p, rapportée au niveau fondamental 
de ces particules est : 


En = No T | 


d Log %) (30) 
nn v' 


d Log T 


Séparons l'énergie de translation. Nous obtenons: 


3 
E, = NÀT É “e r) (31) 
en posant : 
 d Log Zÿ À 
ŸP = 3 Log T° Ge 


a) MorécuLes. — Dans l’expression (17) de Z, 
posons : 


Rev, 


AT 
By ho, 
1—e KT 


(33) 


on obtient aisément : 


PS1 EF Pa Oo + ... 


| 
Yu tr OO Er de (C2 En mL Un 
(34) 

avec : 

ll ke oy l 

| Pn KL ho, 

CAD Er 
| er (35) 
RUACURT fe 


Dans l'expression de Yy, le terme 1 représente 
la contribution de la rotation, qui est totale au- 
dessus de 20 ©K ; le second terme, celle de la 
vibration, qui est comprise entre 0 et 1; le troi- 
sième, celle de l'excitation des électrons ; ce terme 
passe par une suite de maxima pour des valeurs 
de T voisines de hc v,/k, he v:[k, … le plus grand 
maximum (pour 7’ voisin de hc v,/k, soit 11 000 ©K 
pour O:, 70 000 °K pour N,) étant de l’ordre de 0,3. 


b) AToMESs. — En posant : 


he 
Va 


On = gne 
on a : 
PeAGAANs Ga NE AL ee 
nn te (36) 
SR IAO ET 0: io 
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Y, reste toujours compris entre 0 et 1 et passe 
par des maxima pour des températures voisines 


de hcv,lk, hcv;lk 


(premier 


23 000 °K pour O, 30 000 K pour N). 


maximum vers 


c) EXPRESSION DE L'ÉNERGIE TOTALE. — Elle 
sera obtenue en additionnant les valeurs de Æ, 


relatives aux neuf espèces de particules 


pour les 


5) 
électrons, E, —+N, kT | et en tenant compte des 


2 


différences des niveaux fondamentaux des atomes 
et des molécules, c’est-à-dire des énergies de 


réaction. On a ainsi : 


1000 


\ 
au 
3 
.S 
D 
ÿ 
500 
0 7 5 10 K/y 
F1G. 1. — Énergie interne (trait plein) et énergie due au 


choc (tirets) de 1 mol.g d’air en fonction de V,/V. 


TABLEA 


où : 
[OT [NI 


[g] — Sre pr POUR la dissociation (W,, W,) 


[0#],[0+],[0#], [N+] pour l’ionisation (W,, W3, W;, Wa). 


Les valeurs de Æ pour les divers couples 
(T, V,1V) sont portées dans le tableau II. On a 
représenté sur la figure 1 les variations de E en 
fonction de V,/V pour une température 7 donnée 
et sur la figure 2 les variations de Æ en fonction 
de T pour une compression V,/V donnée : la 
pente de cette dernière courbe donne la chaleur 
spécifique molaire de l’air à volume constant. 


5. Augmentation d’énergie interne au passage 
d’une onde de choc. — Nous rappelons brièvement 
les équations du choc. L'indice 0 désignant le 
milieu non perturbé, les lettres sans indice le milieu 
perturbé, p étant la pression, ? la densité, D la 
vitesse du front et Ü la vitesse d’écoulement 
mesurées par un observateur immobile, la conser- 
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K/v. 1 


Kilojoules 


2500 5000 7500 10000 


FiG. 2. — Énergie interne de 1 mol.g d’air 
en fonction de 77. 


UT 


ÉNERGIE INTERNE DE 1 MOL.G D’AIR A T°0K SOUS UNE COMPRESSION VV. 
VALEURS EN KILOJOULES 


TK VIV 


3 900 
4 500 
5 500 
6 500 
7 500 
8 500 
9 500 
10 500 
11 500 


161,41 
231,86 
293,75 
374,46 
498,02 
666,96 
857,77 
1 035,06 


_- 
/ 


150,18 
217,60 
280,15 
349,52 
447,22 
582,83 
748,99 
923,42 
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vation de la masse, de la quantité de mouvement 
et de l’énergie, s'expriment par les trois relations : 


Po D = p(D — U) : 
Po Co D p + po DID U) (38) 


Pos lp gl : lp rm 
Es ++SD=E+l +5 (DU; | 


Les énergies internes introduites dans la troi- 
sième équation sont relatives à 1 gramme de fluide. 
Plus loin, nous les remplacerons par les énergies 
molaires. De ces trois équations, on déduit aisé- 
ment : 


ie DES Lu 
D? P Po Po Fo Po (39) 
Tir 
d P Vo 
y ; 
= D (1 F (40) 


où : M, — masse moléculaire de l’air — 28,86 
et la relation d’Hugoniot : 
1 1 il 
E — = — — — — 
Eo 9 (P + Po) de -) 
1 
= 2m, À Pole M. 


En revenant aux énergies internes de 1 mol.g, 
On à : 


EE, = lp + po) (Vo— V). (41) 


Nous pouvons exprimer la pression p en fonction 
. de T et de V, car nous connaissons l’équation 

d'état de l’air : 
p = nombre de particules par cm x XT 

=: NU + €) KT. (42) 

4 

Nous écrivons ainsi que l’air derrière le front de 
choc est un gaz parfait. Ceci est légitime, car le 
volume spécifique (10 fois plus faible environ que 
celui de l’air normal) est très supérieur au 


covolume. Nous écrivons aussi que l'équilibre entre 
les trois degrés de liberté de translation est atteint : 
l’on sait que c’est chose faite à deux ou trois libres 
parcours moyens derrière le front. D’autre part : 


NAT, 5 


Po Por Eo = 3 NokTo 
de sorte que : 
pi T Vo : 
Do = (1 + €) 1 Y . (43) 


En portant dans (41), on obtient l'expression de 
l’énergie interne de 1 mol.g d’air après le passage 
de l’onde de choc : 


Es (: 7) |: + (1 + de Fi] 


Les conditions choisies pour le milieu aval sont: 
To = 2730K po = 1,014.10$ baryes 
Fi 22%15;cm?. 


6. Grandeurs thermodynamiques dans le front 
d’une onde de choc. — En exprimant que cette 
valeur de Æ doit être égale à celle donnée par la 
relation (37), on obtient une relation entre les deux 
paramètres jusqu'alors indépendants, T'et V,/V. 

Pratiquement, on calcule les valeurs de ÆE 
données par la relation de (44) pour les mêmes 
couples (7, V,/V) que ceux qui ont servi au calcul 
de (37). 

Ces valeurs sont portées dans le tableau TIT et 
représentées sur la figure 1 (courbes en tirets). 
Deux courbes £ = f(V,/V) tracées à température 
constante se coupent en un point, qui représente la 
valeur de V,/V relative à cette température. 

Nous obtenons ainsi neuf couples (7, V,/V) qui 
sont valables derrière des ondes de choc de gran- 
deurs variées. De là nous calculons la pression p 
par (43), ce qui nécessite une interpolation entre 
deux valeurs voisines de (1 + e), puis les vitesses D 
et U par (39) et (40). Toutes ces grandeurs, ainsi 


Fr 
TABLEAU III 


ÉNERGIE DE 1 MOL.G D’AIR DANS LE FRONT D’UNE ONDE DE CHOC 
DE TEMPÉRATURE 7 0K ET DE COMPRESSION V,/V 


To = 273 0K 


VALEURS EN KILOJOULES 


Po = 1 atm. 
HP NN HS 7 d 

3 200 95,05 

4 500 126,55 165,55 

9 900 164,29 215,18 

6 500 202,51 265,91 

7 500 242,10 318,49 

8 500 290,35 380,67 

9 500 351,40 458,51 

10 500 424,79 952,67 


11 500 504,72 657,21 


9,5 10 1 12 
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TABLEAU IV 


EURS DE RANDEURS THERMODYNAMIQUES DAN E FRONT D UNE 
VAL S S G S (e) I S DANS LE FRONT D'UNE ONDE DE CHOC DANS L’AIR 


Conditions initiales : 16 = 27930K  p, = 41 aim VV, —122 415 cm*. 


ThaKR Vo1V DD) D NU) U Ge (J 
m/s m/s ohm—1.cem—t1 
3 500 75153 92,87 2 900 2 495 10? 
4 500 8,272 145,1 3 093 3 159 1820/8405 RSS CE 
5 500 9,089 208,7 4 287 3 815 17 ASE DS OPEL 
6 500 9,424 268,8 4 858 4 342 18015040 BR MOEE 
7 500 9,499 334,5 5 410 4 855 4,995 -X 10 _0,20 
8 500 10,326 423,0 6 066 5 478 LATE 0,71\ 
9500 411,013 541,1 6839 6218 5,143 X 10—1 2,0 
10 500 11,632 685,3 7 678 7 018 LACET OS 4,7 
11 500 12,034 841,6 8 497 DR 2 CALE MIORS 0 
% = pp, = rapport des pressions derrière et devant le front. 
D — vitesse de front. 
U — vitesse d'écoulement. 
45 — taux d’ionisation dans le front. 
o — conductivité électrique du front {mobilité des électrons dans l’air standard = 5 200 em°?.volt=1.sec—1), 


TABLEAU V 


COMPOSITION DE L'AIR DANS LE FRONT D’UNE ONDE DE CHOC. 


Composition normale (p, — 1 atm 


T'oK 
3 200 
4 500 
5 500 
6 500 
7 500 
8 500 
9 500 
10 500 
14 500 


É-12500 


T °K 


ce 
| 7500 


D m 
/s 
0 2500 5000 7500 
Fic. 3. — Température derrière le front 


[O2] 

0,1993 

0,1481 

7,543 x 1072 
2,910 X 10—2 
1,455 X 10? 
5,466 X 10— 
3,035 x 10—5 
1,886 X 10—3 
1,262 x 10—5 


Gaz parfait 


Y=1#4 
De 


[0] 


0,02549 
0,1277 
0,2731 
0,3658 
0,4009 
0,4131 
0,4179 
0,4202 
0,4215 


[N2l 


0,7880 
0,7877 
0,7848 
0,7728 
0,7409 
0,6795 
0,5867 
0,4703 
0,3474 


en fonction de la vitesse de l’onde de choc. 


= 


273 0K): [O;,] — 0,212 
[N] 
2,041 X 10—5 4,6 
6,879 X 10—4 1,404 
6,380 X 10—3 1,989 
3,035 X 10—2 1,217 
9,492 X°40—2 :: 5,782 
0,2171 2,264 
0,4027 7,260 
0,6354 1,931 
0,8812 4,384 


0 


2500 


[N:] = 0,788. 

[e] 1+e 
OS 1,01276 
4 0 1,06420 
x 10—8  1,13976 
x 10—5  1,19810 
x 10—5  1,24765 
<DEs 1,31553 
RIDE 1,41175 
1077 1,53170 
X 10—%  1,66009 


5000 


Gaz parfait Y= 1,4 


7500 10009 


F1G. 4. — Compression derrière le front 
en fonction de la vitesse de l’onde. 
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que le taux d’ionisation x,, sont portées dans le 
tableau IV. Des interpolations en V,/V permet- 
tent de calculer la composition de l’air derrière le 
front (tableau V). 


0 2500 5000 7500 


F1G. 5. — Pression derrière le front 
en fonction de la vitesse de l’onde. 


Fic. 6. — Composition de l’air derrière le front 
en fonction de la vitesse de l’onde de choc. 


Nous avons choisi, pour faciliter les calculs, la 
température comme paramètre. Il est clair qu’une 
grandeur physique beaucoup plus accessible à 
l’expérience est la vitesse du front, puisque le 
passage de l’onde s'accompagne de phénomènes 
lumineux faciles à enregistrer. C’est pourquoi, sur 
les figures 3, 4, 5, 6, 7, nous avons représenté les 
variations des différentes grandeurs en fonction de 
cette vitesse. La figure 8 représente la compres- 
sibilité de l’air par onde de choc, c’est-à-dire ce que 
l’on a convenu d'appeler l'Hugoniot de l’air. 


La vérification expérimentale de ce calcul peut 
se faire en mesurant, en fonction de la vitesse D du 
front d’ondes de choc d’intensités différentes, une 
autre grandeur physique. Deal [5] à ainsi mesuré 


Fic. 7. — Taux d’ionisation 
en fonction de la vitesse de l’onde de choc. 


Gaz parfait 


0 500 
F1. 8. — Hugoniot de l'air. 


simultanément D, vitesse du front, et U, vitesse 
d'écoulement, et obtenu un accord satisfaisant 
avec des calculs semi-théoriques effectués par 
Duff [3]. Ses résultats confirment en particulier 
que la valeur de l'énergie de dissociation de l’azote 
effective dans une onde de choc, c’est-à-dire une 
perturbation rapide, est 9,764 eV, Toutefois, les 
conditions de ses expériences ne lui permettaient 
pas d'atteindre des vitesses de front supérieures 
à 4 500 mys. 

Dans le domaine des ondes de choc intenses 
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(vitesses de l’ordre de 6 000 à 10 000 m/s), nous effec- 
tuons actuellement les mesures d’une autre grandeur 
physique, qui est le taux d’ionisation. Cette gran- 
deur nous est accessible par l’intermédiaire de la 
conductivité électrique, ce qui suppose que l’on 
connaisse la mobilité des électrons dans les condi- 
tions de température et de pression élevées qui 
règnent derrière le front. Si l’on admet que l’énergie 
acquise par un électron entre deux collisions dans 
un champ électrique extérieure est très inférieure à 
l’énergie moyenne d’agitation des particules, on 
peut écrire : 

U = Uo-VTITo-(polp) (45) 
uo désignant la mobilité de l’électron à (po, To) 
et u la mobilité à (p, T). Dans le cas de pressions 
très élevées, les collisions sont assez nombreuses 
et le parcours moyen des électrons assez faible pour 
que la formule (45) soit valable. D’autre part le 
nombre x d’électrons par cm° est relié au taux 
d’ionisation par : 


n = No.le]|V = (NolV).xo. 1 + €) = (pJXT)xo. (46) 


Non 
La conductivité électrique s’exprime alors par: 


No etc to MIE 


La valeur de u, est tirée des résultats expéri- 
mentaux de Crompton, Huxley et Sutton [6], repre- 
nant eux-mêmes les anciennes mesures de Nielsen 
et Bradbury. On trouve ainsi, à 20 % près : 


(47) 


ORALE EN 


o = 5 200 cm?.voltr 1.sect. 
Cette valeur, portée dans (47), donne : 


6 — 2,2.104.æ, V1[L'ohm—!.cm—t. (48) 


Les valeurs de 6 déduites de cette formule sont 
portées dans le tableau IV. Comme elles dépendent 
directement de la valeur x,, elles ne sont exactes 
qu’à 20 % près. Toutefois, compte tenu de cette 
incertitude systématique et des erreurs expéri- 
mentales probables, la croissance quasi-exponen- 
tielle du taux d’ionisation en fonction de la vitesse 
du front doit permettre un contrôle aisé de la 
justesse des calculs présentés. 


Manuscrit reçu le 14 mars 1958. 
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ÉMISSION MOLÉCULAIRE ET RECOMBINAISON ATOMIQUE DANS L'’AZOTE 


Par J. AKRICHE, L. HERMAN, 


Laboratoire des Recherches Physiques à la Sorbonne. 


et H. GRENAT, 


Observatoire de Meudon. 


Résumé. — Les auteurs calculent la répartition de la population sur les différents niveaux de 
vibration d’une molécule diatomique excitée, formée par collision de deux atomes. Pour la tempé- 
rature de 3 000 °K, on obtient un maximum de population pour v’ = 4 ou 5. Pour la température 
de 300 °K, le maximum est déplacé vers v’ = 9 ou 10. 

Ces données sont comparées à celles obtenues par l’étude expérimentale de la phosphorescence 
du premier système positif de l’azote émis dans un mélange d’hélium et d’azote. Pour interpréter 
les observations, on propose un Prop eAus faisant intervenir l’état métastable de l’hélium atomique 


et l'ion moléculaire d’azote. 


Abstraet. — The authors calculate the distribution of the population of the different vibra- 


tional levels of an excited diatomic molecule formed by the collision of two atoms. 


For the tempe- 


- rature of 3 000 ©K à maximum is obtained for v’ = 4 or 5 ; at 300 °K the maximum is for the 


levels »’ — 9 or 10. 


These calculations are compared with the experimental results obtained from the afterglow of 


the first positive system of nitrogen emitted in the presence of helium. 


The proposed explanation 


of these observations involves the metastable state of atomic helium and the molecular ion of 


nitrogen. 


Introduction. — Le déclin de l’ionisation et de la 
luminescence des gaz rares et de quelques gaz 
diatomiques- a été étudié, ces dernières années, à 
l’aide de techniques très différentes. Les résultats 
obtenus dépendent de la grandeur de la concen- 
tration électronique initiale du milieu ionisé et de la 
pression du gaz. 

Suivant les idées récentes, la disparition des 
électrons serait régie par le processus de recom- 
binaison dissociative [1]. 


M} + électron —> M’ + M’. (1) 


Celui-ci correspond, dans le cas des atomes, à la 
recombinaison diélectronique (qui est l’inverse de 
l’autoionisation). La section eflicace est faible pour 
cette réaction qui a lieu généralement pour les 
électrons ayant une énergie cinétique élevée, dont 
le nombre est petit dans un gaz électronique pos- 
Sédant une distribution thermique. Dans le cas de 
molécules diatomiques, on admet qu’il existe un 
nombre élevé de molécules instables à courbe de 
potentiel répulsive et possédant un point d’inter- 
section avec la courbe de potentiel de l’ion. 

- L'énergie des électrons dans lé processus (1) peut 
être faible, ce qui augmente son rendement. 

L'émission des raies atomiques des gaz rares est 
due à l’émission spontanée des atomes suivant la 


réaction, ë 
M" M + hv. (2) 


Le taux de la recombinaison radiative, 
M+ + électron — M’ + Av, (3) 


est très faible par rapport à celui de (1). Le rapport 
des coefficients est de l’ordre de 105 et (3) peut 


généralement être négligé devant (1), dans le cas 
de basses températures du gaz. 

Cette interprétation séduisante n’est pas satis- 
faisante pour les gaz diatomiques, comme Nÿ 
et Hef. En effet, lorsque la pression atteint plu- 
sieurs millimètres de mercure, le spectre émis n’est 
pas atomique, mais moléculaire. Pour interpréter 
cette émission on est conduit à admettre l’existence 
d’une recombinaison moléculaire radiative, 


MX + électron — M,(v') (4) 
M, (v’) > M, + hv. 


Le coefficient de (4) ne peut être très différent de 
celui de (3). Le déclin de l’ionisation serait déter- 
miné par (1) et non par (4). 

Les mesures quantitatives faites par Bryan, Holt 
et Oldenberg [2] sur l'émission du spectre de l’azote 
ne peuvent être interprétées uniquement à l’aide 
des réactions (1) et (4). La comparaison du nombre 
de quanta émis au nombre de recombinaisons du 
processus (4) conduirait à admettre que la dispa- 
rition d’un ion et d’un électron est accompagnée de 
l'émission de 105 quanta de rayonnement, ce qui 
est difficile à expliquer. Dans ce cas, on doit envi- 
sager pour l’émission une recombinaison atomique 
par choc double ou par choc triple. Pour les pres- 
sions de plusieurs mm de Hg, utilisées dans nos 
expériences et dans celles de Bryan, Holt et Olden- 
berg, les chocs triples jouent un rôle important et 
on aurait, 


M'+M + M, + M,(v’) + M (5) 
M, (v') — M,(v) + hv. 
Dans le cas de N, et He, il suffira de retenir les 
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réactions (1) et (5). Pour obtenir le rendement 
lumineux observé il suffit d'admettre que 1073 
ou 10-4 du nombre total de molécules est dissocié 
en atomes. 

I] ne serait pas difficile de vérifier lequel des deux 
processus (4) ou (5) joue le rôle prépondérant. Pour 
cela, il suffit d'appliquer au gaz rayonnant une 
impulsion électrique brève ayant pour effet d’aug- 
menter le champ électrostatique ou la température 
électronique. Les radiations dues aux proces- 
sus (3) et (4) diminueront d’intensité, mais non 
celles du processus (5). 


Distribution de la population des molécules sur 
les niveaux de vibration. — Dans les réactions (4) 
et (5), la probabilité de formation d’une molécule à 
l'état M;(v') avec un nombre quantique de vibra- 
tion v' dépendra du facteur de superposition de 
Franck-Condon, c’est-à-dire de la position relative 
des courbes d’énergie potentielle des molécules 
stables M, et MX ou de M; et de la molécule 
instable (M' + M). La première hypothèse a déjà 
été discutée précédemment [3]. Nous allons exa- 
miner les résultats que l’on obtient dans le cas de 
la deuxième hypothèse. 


Recombinaison atomique. Collisions centrales. — 
On manque de données précises concernant les 
courbes de potentiel répulsives pour les molécules 
d’azote. Nous avons donc utilisé les courbes de 
potentiel de la molécule instable H, (b 32,) et celle 
de la molécule stable H, (a ?Z,). 

Considérons, par exemple, deux atomes d’hydro- 
gène qui entrent en collision. Ils ont, à une grande 
distance l’un de l’autre, une vitesse relative v et le 
système possède à l'infini une énergie totale v?/2, 
où est la masse réduite. Si on considère la courbe 
de potentiel Z de l’état répulsif de la molécule, les 
atomes se rapprochent, en gardant leur énergie, 
suivant une droite horizontale AB (fig. 1). 

Aù point A, leur énergie cinétique sera nulle et 
leur énergie potentieelle égale à uv?/2. La pro- 
babilité pour qu'il y ait émission d’un quantum est 
proportionnelle à la probabilité de transition de AB 
vers un niveau inférieur v”. Pour connaître la popu- 
lation N(v") de ce niveau v”, il faudra faire inter- 
venir, par conséquent, la distribution des vitesses 
des atomes de l’état répulsif. 

Pour une température donnée 7, on admettra 
une distribution maxwWellienne des vitesses, donc 
proportionnelle à v?.exp (— uv?/2kT). Le nombre 
de molécules arrivant au niveau v” sera donc égal 
au nombre de quanta émis, et on peut écrire 


ST az 
N(v") = ff à A 
— cte f vXwr.v?.exp (— v?/2KT) [fa Ver dr]° dy. 


où Ÿ, est la fonction d’onde de l’état répulsif au 
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niveau Z, et W,,, celle de l’état inférieur au niveau v” 
considéré. 

Il faut donc calculer la probabilité de transition 
de AB vers v”. 

4) Condon a mis au point une méthode dite de 
réflexion qui consiste à admettre la fonction d’onde 
de l’état répulsif comme constante au point d’inver- 
sion, Dans ces conditions, la probabilité de tran- 
sition ne dépend plus que de la fonction d’onde du 
niveau inférieur. Par conséquent, on peut écrire 
N(v") = >> = cte.v°.exp (— uv?[2RT) F”2(v’).v8(v").dv. 


\ 


Ur) 


Pré. 4: 


A l’aide de cette méthode, nous avons calculé les 
populations des niveaux 0" = 0, 1, 2, 3, 4, 5, 6, 7, 
de l’état stable a %Z,, en supposant que cet état et 
Pétat répulsif b 32, situé au-dessus de l’état stable, : 
ont la même limite de dissociation. Les fonctions 
d'onde des niveaux v” —0, 1, 2, ... 7 ont été 
calculées par la méthode de distorsion indiquée par 
Pillow [4]. 

Pour calculer la distribution des vitesses, nous 
avons admis une température de 3 000 0K. 

Les populations ÆV(v”) trouvées dans ce cas, 
figurent dans le tableau 1 : 


TABLEAU I 


POPULATION DES NIVEAUX DE VIBRATION 
DE L’ÉTAT à 92, OBTENUE PAR LA MÉTHODE DE RÉFLEXION 


Du 0 1 2 3 A 5 6 7 
Ny = 0,8 5,0 16 29 40 40 35 30 


Ne? 


La distribution dela population est caractérisée 
par un maximum très plat s'étendant entre v" — 4 
et 7, alors que le maximum de la distribution des 
vitesses se situe vers les niveaux v” —7 et 8. 

2) Nous avons ensuite essayé la méthode plus 
exacte qui tient compte de la forme de la fonction 
d'onde de l’état instable. Nous avons utilisé la 
fonction d’onde de W. K. B. 


Jr) A1 J(E — V)V4.sin 2x | [” Vème — V) arf], 
re 


où E est l'énergie totale, 

ro la distance internucléaire au point d’inversion, 
au niveau instable, | 

V, l'énergie potentielle à la distance inter- 

nucléaire r. 

Ainsi, nous avons pu tracer les courbes des fonc- 
tions d’onde pour différentes énergies cinétiques des 
atomes. Pour l’état inférieur a 32, nous avons 

- utilisé les fonctions d’onde de Pillow. Les calculs 
ont été faits pour v” = 4, 5, 6, 7. Nous n’avons pas 
pu considérer les niveaux 0, 1, 2, 3, car le terme de 
distribution des vitesses est très petit dans cette 
région. 

Les populations des niveaux V(v”), calculées de 
cette manière sont données dans le tableau II. 


TABLEAU II 


POPULATION DES NIVEAUX DE VIBRATION 
DE L'ÉTAT @ 2, OBTENUE 
PAR L'EMPLOI DE LA FONCTION D’ONDE W. K. B. 


vu = 4 5 6 7 
Ny =40 43 54 66 


La comparaison des tableaux I et II montre que 

l'emploi de la fonction d’onde de W. K. B. pour 

_ b$E, a pour effet de déplacer le maximum, qui 
reste très plat, avec des valeurs de v” > 7, 


Effet de la rotation moléculaire. — Dans les 
calculs précédents de la population des différents 
niveaux de vibration de l’état stable, on a admis 
implicitement que les collisions sont centrales, 
c’est-à-dire que les deux atomes se rencontrent en 
se déplaçant sur une ligne droite. Cela revient à 
retenir, comme courbe de potentiel, la courbe clas- 
sique de Morse. En réalité, la plupart des collisions 
atomiques possèdent un paramètre de collision p, 
différent de zéro. La courbe de potentiel de Morse 
U{r) correspondant à p —0. Au paramètre p 
correspond un moment angulaire P = pv, p, où p 
est la masse réduite et v,, la vitesse initiale rela- 
tive des deux atomes. 

À la place de l'énergie potentielle U,(r), 
pour p = 0, on a, en général : 

Un(r) = Us(r) + P?[2ur? = énergie potentielle 
effective. P?/2ur? est l'énergie cinétique d’un 
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corps de masse u tournant avec le moment angu- 
laire P sur un cercle de rayon r, 

Nous avons tracé, sur la figure 2, une série de 
courbes de potentiel effectives U,(r), pour difté- 


distance inlernucléaire x (À) 


16 20 24 28 


ü 


Pre 02; 


rentes valeurs du paramètre p. Pour choisir les 
limites de variation de op, considérons l’énergie 
cinétique de rotation, 


Fan = P?[2urt = (4 v8/2) (e?1r*), 


et prenons pour valeur minimum de la distance 
internucléaire la valeur r, qui correspond au point 
de rebroussement. Si l’énergie cinétique 
TARN EENONE) PRE 
on à : 
pale MN 


L'examen des courbes de potentiel effectives 
montre qu’elles présentent une sorte de barrière 
vers 7 — 1,4 à 1,5 À. D’autre part, les fonctions 
d'ondes s’annulent pratiquement pour r > 3 À. 
De sorte que, 


DIR S rIS 0 À 
1,7 À2 < r? <:9,0 A2, 
Ceci nous situe le domaine de variations de p?. 


Pour bien encadrer les valeurs utiles, nous avons 
fait varier 92 entre 0,2 À? et 40 A2, 
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Dans la suite du calcul, nous avons admis 
T = 3 000 ©K, d’où 4T — 2 090 cm1, Comme 
dans le calcul précédent, on écrit que la population 
d’un niveau v" de l’état électronique stable a E, est, 
en utilisant la méthode de Condon, 


N(v”,a3Z,)= cte.v?(v,).v?.exp (—u v?/2KkT).W?(v;). 


Tous calculs faits, on obtient les populations des 
niveaux de vibration du tableau III, ligne 2. 


TABLEAU III 


POPULATION DES NIVEAUX DE VIBRATION 
DE L'ÉTAT @ $Z4 OBTENUE EN TENANT COMPTE DE L'EFFET 
DE LA ROTATION MOLÉCULAIRE (7 = 3 000 0K) 


v' — 0 1 2 3 4 5 6 7 
N{v'} 0,2 - 08 26/2140: 791 NI LT SO 
rotation 

PA CC ONE RS RE CRE ON EE TMC T 
Condon 

N(v") : h LOT LS NE TRE 60 
W. K. B. 


A titre de comparaison, nous avons reproduit 
dans les troisième et quatrième lignes, les nombres 
des tableaux I et Il. On a pris arbitrairement 
N(4) = 40 pour les trois cas. On voit que l’allure 
générale de la distribution ne change pas lorsque 
l’on tient compte de la contribution de la rotation 
à la courbe de potentiel. La seule différence con- 
cerne les premiers niveaux de vibration. 

Les données du tableau III se rapportent à la 
température atomique 7 = 3 000 ©K. Il est facile 
de refaire ces mêmes calculs pour d’autres tempé- 
ratures. Lorsque celle-ci baisse, le maximum de la 
population se déplace vers les plus grands nombres 
quantiques. Pour 7 = 300 ©K, le maximum est 
situé entre v = 9 et v = 10. 


Étude expérimentale 


Pour pouvoir comparer les résultats des calculs 
précédents aux données expérimentales, nous avons 
repris l’étude de l’émission de l’azote. On sait que 
dans des conditions expérimentales bien déter- 
minées, l’azote, préalablement excité par une 
décharge électrique, montre une phosphorescence 
dont la durée peut varier de 104$ à 10 minutes 
(et plus). La composition du spectre de phospho- 
rescence montre que l’on a, en réalité, affaire à 
différentes types de réaction. L'émission de courte 
durée donne un spectre ne présentant pas de sélec- 
tivité d'intensité marquée dans les séquences du 
premier système positif. Par contre, la phospho- 
rescence de longue durée, connue sous le nom de 
« postluminescence de Lewis Rayleigh », montre 
une émission très sélective, les bandes à 2! —6 
et v' = 12 du premier système positif possédant 
une très forte intensité relative. 
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Nous avons utilisé un troisième type de phos- 
phorescence, lequel, tout en étant de longue durée, 
ne montre pas de sélectivité marquée de popu- 
lation des niveaux de vibration de l’état initial. 


Dispositif expérimental. — Le tube à décharge 
que nous avons utilisé était constitué par un ballon 
en pyrex de Î litre muni d’une électrode de tung- 
stène et rempli d’hélium sous la pression de 4 mm 
de Hg et d’azote sous une pression estimée à envi- 
ron 0,1 mm de Hg. L’azote était préparé à partir 
du nitrure de sodium décomposé sous l’influence 
de la chaleur. La dessication de ce mélange était 
assurée par de l’anhydridé phosphorique contenu 
dans une tubulure latérale en communication avec 
le ballon. 

L’excitation du gaz est obtenue à l’aide d’une 
décharge électrique de haute fréquence traversant 
l’électrode en tungstène pendant quelques secondes. 
L’intensité et la durée de la phosphorescence ne 
semblent pas dépendre de la durée de la décharge, 
qui a toujours été supérieure à une seconde. Au 
début de la postluminescence, la luminosité est très 
uniforme dans tout le ballon et l’on ne peut noter 
aucune diminution de l'intensité à la surface inté- 
rieure de la paroi : celle-ci n’exerce aucune action 
de destruction sur les particules actives excitant 
l'émission, 


Dispositif stroboscopique. — La phosphores- 
cence dure plusieurs minutes et il était suffisant, 
pour photographier son spectre, d’utiliser un stro- 
boscope mécanique lent. Un basculeur à mercure 
était actionné par un moteur synchrone dont l’axe 
de sortie était animé d’une rotation d’un tour par 
seconde. Un disque convenablement découpé per- 
mettait d’éliminer la décharge ; il était placé 
devant la fente du spectrographe type A de la 
S. G. O., à optique de verre. Quelques photo- 
graphies dans l’infrarouge proche ont été prises 
avec un spectrographe ouvert à f/1, de 50 mm de 
distance focale. 

Les spectres obtenus montrent la présence de 
bandes positives et négatives de l’azote. 

Nous avons mesuré l'intensité relative des 
diverses bandes, par photographie, suivant la tech- 
nique de Ja photométrie photographique bien 
connue. À partir de ces intensités relatives, /(v’, 0"), 
nous avons déterminé la population M(v') de 
l’état B SII(v'), par la relation : 

N{(v'):=. Cte.J(v!,- 01) /P(v', 07) So 
où v(v', v”) est la fréquence et P(v', v”), la proba- 
bilité relative de transition du niveau supé- 
rieur B #II(v') vers le niveau inférieur À $Z(v"). 
Les valeurs numériques des probabilités de tran- 
sition avaient été calculées par Pillow, Nicholls et 
ses collaborateurs [5]. Les valeurs de (v') figurent 
dans le tableau IV. 
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TABLEAU IV 


PoPpuLAtION V(v') DES NIVEAUX DE VIBRATION DE L’ÉTAT 
INITIAL B SIT pe LA TRANSITION B 8II — À 3X pr LA 
MOLÉGULE N;,, AU DÉBUT DE LA PHOSPHORESCENCE. 

12 

31 


= 5 6 7 8 9- 
ND) = 120 29 53 79 gt 


10 
79 


11 
67 


Ces populations peuvent être comparées à celles 
calculées précédemment dans l'hypothèse d’une 
recombinaison atomique. Les trois distributions 
calculées ont été portées dans la figure‘3 (courbes I, 
IT et IIT). 


AE 


50 


fe] 2 7 6 8 10 12 


Fire, 3. — Population des niveaux de vibration de l’état 
électronique initial B $IT du premier système positif. Les 
points de la courbe IV représentent les données obtenues 
par les mesures photométriques et les courbes I, IT et IIT, 
celles obtenues par le calcul. Les courbes I et II se rap- 
portent aux collisions centrales. Dans la courbe III on 
tient compte de la rotation moléculaire. 


La courbe IV est celle du tableau IV. Étant 
donné que la courbe de potentiel répulsive de N; 
n’est pas connue et qu’elle a été remplacée par 
. celle de la molécule H;, on peut conclure, en raison 
de ces approximations, que l’accord entre le calcul 
(courbe IV) et l'expérience est assez satisfaisant, 
si on utilise la température TZ — 300 °K. Nous 
avons, toutefois, cherché une confirmation de 
l'hypothèse de la recombinaison atomique dans 
l’étude expérimentale de l’intensité émise. 


Coefficient de recombinaison. — En vue de déter- 
miner le coefficient de recombinaison, nous avons 
mesuré l’émission en valeur absolue du gaz phos- 
phorescent. Pour le faire, nous avons utilisé un 
étalon secondaire d’émission constitué par une 
lampe à filament de tungstène ayant la tempé- 
rature de couleur de 2 160 ©K [6]. Comme récep- 
. teur nous avons utilisé une cellule photoélectrique à 
19 étages de Lallemand, alimentée par une tension 
stabilisée. Le courant photoélectrique était admis 
à l'entrée d’un adaptateur d’impédance dont le 
courant de sortie alimentait un enregistreur à 
plume. Pour isoler diverses régions spectrales, on 
disposait de verres colorés placés devant la cellule 
photoélectrique. Nous avons pu ainsi réaliser la 
mesure : 
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1) de l'intensité et du nombre de quanta émis 
par le gaz par em? et par seconde, et 

2) de la variation de l’intensité en fonction du 
temps. 


Mesure de l’intensité absolue. — Nous avons 
placé le ballon et la lampe étalon secondaire à le 
même distance, d'environ 5 mètres, de la cellule 
photoélectrique munie d’un filtre approprié. Pour 
réduire l’intensité de la lampe à filament de tung- 
stène, nous avons utilisé un coin photométrique de 
densité optique connue. On obtient ainsi ; pour la 
région comprise entre 5 950 À et 6500 À, le 
nombre de quanta émis par seconde et par em? 
par le gaz lumineux. Les résultats obtenus figurent 
dans le tableau V. 


TABLEAU V 


NOMBRE DE QUANTA ÎV ÉMIS PAR CM? 
ET PAR SECONDE PAR LE GAZ PHOSPHORESCENT 
DANS LA RÉGION 5 950 À-6 500 À 
(PREMIER SYSTÈME POSITIF) 


30 
0,47 


15 45 60 


1,28 


Temps en secondes : 0 


102 N (car 31): 3,37 0,24 0,12 


Loi du déelin. — Nous avons étudié également 
la variation avec le temps { de l’intensité relative J, 
en lumière rouge, verte, bleue et en lumière totale. 
Grâce à la grande sensibilité de la cellule photo- 
électrique, il était possible de suivre le déclin du 
rayonnement pendant plus de 7 minutes. Les résul- 
tats obtenus sont reproduits sur les figures 4 à 9. 

Nous avons utilisé deux échelles pour les inten- 
sités : l’échelle  logarithmique. (log 7) et 
l’échelle 7-12, Si la réaction était monomoléculaire, 
(log Z)/i serait constant et les points des figures 5, 
6 et 8 se placeraient sur une droite. Si le processus 
était bimoléculaire (1-12)/1 serait constant et Îles 
points des figures 4, 7 et 9 se placeraient sur une 
droite. 

L’examen des figures ci-dessus montre que l’on 
ne peut représenter toutes les observations par un 
seul processus. Les figures 5, 6 et 8 montrent que la 
fin de l’émission correspond à une réaction mono- 
moléculaire. Par contre, le début des courbes des 
figures 4, 7 et 9 semblent se rapporter à une réaction 
bimoléculaire. Ces courbes pourraient également 
être représentées par deux segments de droites. 
Toutefois, l’interprétation semble bien plus difli- 
cile dans ce cas. On serait conduit à admettre 
l’existence de deux processus bimoléculaires, avec 
une vitesse de réaction beaucoup plus lente au 
début et plus rapide à la fin de la phosphorescence. 
Mais, il serait difficile de trouver deux réactions 
bimoléculaïres ayant un taux de réaction croissant 
avec le temps. 
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(0) 1 2 3 4 5 6 


FiG. 4. — Variation de linverse de la racine carrée de l'intensité lumineuse en fonction du temps. La courbe La 
été obtenue à l’aide d’un écran rouge transparent dans le domaine spectral 6100 À-6 700 À et la courbe IL en lumière 
totale, sans interposition de filtre entre la source et la cellule photoélectrique. La courbe [ semble être composée 
de deux tronçons de droites et la courbe IT de trois tronçons. 


La 


De ANT UT NS ST EN NE TE 

Fic. 5.— Les mesures utilisées ici sont les mêmes que  Fic. 6. — Les points expérimentaux ont été obtenus en 
celles de la figure 2. La courbe I a été obtenue sans filtre lumière verte (4 870 À < x < 6100 À). La partie recti- 
et la courbe II, avec un filtre rouge. La seule différence ligne à la fin de la phosphorescence est bien plus marquée 
vient de l’échelle des abscisses qui est logarithmique. que dans le cas de la figure 5, 


Les points se placent dans leur ensemble sur une courbe 
et non sur une droite. Toutefois, ces courbes ont des 
droites comme asymptotes. 
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0,75 


0,5 


025 


[o] TE 2 3 - 4 


Fig. 7. — Les points expérimentaux de la figure 6 ont été 
retracés en portant en ordonnée 7Z—1/2?, On constate que 
la région correspondant aux temps compris entre 1 el 
3 minutes ne peut être représentée par une droite. 


0 1 e 3 4 


Fi. 8. — Variation de log Z en fonction du temps. Les 
points ne ont été obtenus en lumière bleue 
(4 870 À > x > 4100 À). Ici, la région comprise entre 
2 et 4 minutes est convenablement représentée par une 
droite, mais non le début de la phosphorescence. 


Coefficient de recombinaison. -— Considérons 
d’abord le début de la phosphorescence où le pro- 
cessus semble effectivement être du second ordre. 


S + 


(0) 1 2 3 4 


F1G. 9. — Les données de la figure 8 sont retracées en 
portant en ordonnées 1—1/2, à la place de log.Z. La pre- 
mière partie de l’émission est bien représentée par une 
droite, mais non la seconde. 


Les données quantitatives précédentes per- 
mettent la détermination du coefficient de recom- 
binaison du processus bimoléculaire : 


N,+N —>+N,(B8"Il). (6) 


Il n’est pas nécessaire de préciser la nature des 
particules N;. Nous admettrons seulement que les 
concentrations de N; et N; sont les mêmes el 
égales à N. 

Pour déterminer, à partir des mesures d’inten- 
sité, le coefficient de recombinaison dans le pro- 
cessus (6), nous admettrons que la disparition de 
deux particules N est accompagnée de l’émission 
d’un quantum de rayonnement : 1 — J/hv. On a : 


21 = 2Jfhv = aN, = — ANfdt 
HUE = LUE + (2e) + HAE — 10), 


où Z, et t, définissent les conditions initiales. Ainsi, 
la pente des courbes de la figure 4 permettent de 
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déterminer les coefficient &«. Les valeurs obtenues 
avec le filtre rouge sont : 


œ, = 0,77.10—12 cm8.s—1, pour 0 < £ < 2,2 minutes 


mi — 272401 cms. pOur 1 > 2,2 minutes. 

Discussion. — Les valeurs des coefficients de 
recombinaison obtenues ci-dessus sont assez sur- 
prenantes. Dans leur étude de la phosphorescence 
de courte durée, Bryan, Holt et Oldenberg [2] 
trouvent &« — 2,5.10-8 em.s-1, alors que, pour la 
recombinaison entre les ions et les électrons, ils 


obtiennent « — 2.10-6 em?.s-—1, 


Recombinaison atomique. -—- Suivant le méca- 
nisme exposé dans la partie |, pour la recombi- 
naison avec émission de radiation, par collision de 
deux atomes, on devrait s’attendre à un coefficient 
de l’ordre de 10-1$.cm%.s-! qui est beaucoup plus 
petit que les coefficients trouvés dans le présent 
travail, 

Le coefficient de recombinaison par choc triple 
est de l’ordre de 10-%2.cmf.s-!, Pour la concen- 
tration totale d'environ 1017.cm-% utilisée ici ; le 
coefficient apparent pour une collision double 
serait de 10—15.cm$.s—1 qui est encore beaucoup 
trop petit. : 

L'étude de l’influence de la pression du gaz sur la 
valeur de « fournirait d’intéressantes indications 
sur l’ordre de la réaction. De telles mesures ont été 
faites par Willey, Rayleigh et Kash [7] pour l'azote 
pur. Ces auteurs ont montré que l’émission est liée 
aux collisions triples pour le premier système posi- 
tif, et aux collisions doubles pour le second sys- 
tème positif de N, et les bandes violettes de CN. 

Nous n’avons pas encore étudié l'effet de 
l’accroissement de la pression de l’hélium. 


Recombinaison électronique. — Considérons 
maintenant les processus dans lesquels. inter- 
viennent les électrons. La recombinaison disso- 
ciative, habituellement admise pour expliquer les 
valeurs élevées trouvées par la technique des impul- 
sions de haute fréquence, ne peut être utilisée ici. 
La recombinaison diélectronique donne des coef- 
ficients beaucoup trop petits. Il reste à considérer 
la recombinaison radiative entre les ions et les 
électrons. Nos valeurs (4 — 10-12 em$.s-1) sont en 
accord avec celles que l’on trouve par le calcul 
théorique pour les atomes, et l'hydrogène en parti- 
culier, pour la température ambiante. 

Malgré ces résultats surprenants, considérons de 
plus près cette hypothèse. Lors du passage de la 
décharge électrique, on observe un spectre d’émis- 
sion composé, en dehors de bandes positives, de 
raies He I et de bandes négatives de Nf. En phos- 
phorescence, les raies atomiques sont invisibles, 
alors que les bandes négatives sont intenses. Comme 
on l'avait souligné précédemment, à plusieurs 
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reprises [8], la grande intensité des bandes 
négatives de l’azote excité en présence d’hélium, 
vient de l’action de l’état métastable He(35S,) sur 
les molécules normales N,(1X). Celles-ci se 
trouvent ionisées et excitées au niveau B'?2X,,. En 
résumé, en présence d’atomes métastables d’hélium, 
très nombreux en phosphorescence, les molécules 
d’azote sont ionisées. Il n’est donc pas surprenant 
d'observer un spectre de recombinaison radiative. 

Pour rendre compte des données expérimentales 
ci-dessus, nous sommes conduits à admettre le pro- 
cessus Suivant : es S 


He(%S,) + N,(X ?5,) 


N'(Br 2) = NF IX/ ES) ET) 


NY (X'25,) + électron a N, + hv,, 


recombinaison radiative (8a) 
N; > N, + 4 
recombinaison radiative (8b) 
N, - —> N° + N prédissociation (8c) 
N + N° -> N,(B SII, v’) 
recombinaison atomique. (84) 


Dans cette suite de réactions, nous aurions pu 
identifier N; à N;; dans ce cas, on aurait une 
recombinaison dissociative. Mais, suivant les idées 
récentes, une telle réaction aurait une section efli- 
cace environ 106 fois plus grande que celle trouvée 
dans nos expériences. Nous devons donc admettre 
que le processus primaire est une recombinaison 
radiative, mais que le processus conduisant à 
l'émission est lié à la collision de deux atomes. De 
cette façon, on expliquerait à la fois la distribution 
de la population des niveaux v' et la valeur du 
coefficient de recombinaison obtenue expérimen- 
talement. | 

Rappelons que dans l’azote pur, ou en présence 
de petites quantités d'oxygène, le spectre d’émis- 
sion et le mécanisme d’excitation semblent être 
très différents [9]. Le fait d'observer des modifi- 
cations sensibles dans la répartition spectrale 
lorsque l’azote est dilué dans un gaz rare, s’expli- 
querait par la superposition de l’effet de deux ou 
plusieurs réactions. 


Fin de la phosphorescence. — La réaction est 
bimoléculaire au début de la luminescence en raison 
de la concentration élevée en ions N+(X' ?E). La 
contribution du processus (7) serait peu impor- 
tante pendant les premières minutes. Par contre, 
à la fin de la phosphorescence, le taux de dispa- 
rition de ces ions tend à être égal au taux de for- - 
mation dans (7) et le processus devient mono- 
moléculaire. Le taux de la réaction serait alors 
déterminé essentiellement par la durée de vie de 
l’état métastable He(%S,) qui est très grande et par 
le taux de diffusion vers les parois. Le taux de 


NT 


diffusion trouvé ici n’est pas très différent de celui 
déterminé par Anderson et ses collaborateurs, si 


_ Pon tient compte de la différence des dimensions 


des tubes à décharge et des poids moléculaires de 
l’azote et de l’hélium. 


Ces auteurs donnent pour l’azote pur et la t > 


pérature de 30 0C, la durée de vie appa- 
rente rt, — 19 s. Si l’on tient compte des plus 
grandes dimensions de notre ballon et du fait que 
l'atome d’hélium est 7 fois plus léger que la molé- 
cule N,, on obtient la valeur +, = 585, alors que 
la valeur moyenne expérimentale déduite des 
figures 5, 6 et 8est +; — 66. 


Diffusion. — Un caleul approché de la vitesse de 
diffusion des particules actives montre qu’elles ne 
sont pas détruites à chaque collision sur la paroi, 
mais peuvent subir des milliers de collisions avant 
de retourner au niveau fondamental. 

La durée de vie + constitue la limite inférieure 
pour la durée de vie de la particule métastable 
diffusant vers la paroi. Elle peut fournir des rensei- 
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gnements intéressants sur la nature de l'atome ou 
de la molécule qui est à l’origine de la phospho- 
rescence. 


Conelusion. — Les résultats expérimentaux obte- 
nus dans le présent travail semblent montrer que le 
coefficient de recombinaison entre les ions molé- 
culaires et les électrons sont de l’ordre de 
10-12 .cm$.s-1. Cette valeur est très différente de 
celles obtenues jusqu'ici par des procédés modernes 
d’impulsions en haute fréquence. Nous devons 
attendre que des mesures directes de concentration 
électronique viennent confirmer notre interpré- 
tation. On serait alors conduit à réviser les idées 
actuelles sur la recombinaison moléculaire disso- 
ciative et à rechercher à nouveau l’origine des 
valeurs élevées des coefficients de recombinaison 
obtenues par la méthode d’impulsion (et non par 
les anciennes mesures à l’aide de sondes). 


Manuscrit reçu le 18 mars 1958. 
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Résumé. — Afin de mettre en évidence leffet Delbrück, nous avons poursuivi notre étude [1] 
de la diffusion élastiqué des rayons y de 2,62 MeV du Th C” par les atomes de plomb, de bismuth 
et d'uranium [2]. Nous avons pu déterminer les sections efficaces de-ce phénomène avec une pré- 


cision de 15 à 20 %. 


Abstract. — Elastic scattering measurements were performed for y-rays in order tofind a 


proof of Delbruck scattering [1]. 


and scatterers of lead, bismuth and uranium [2]. 


been found with an accuracy 15-20 %. 


I. Introduction. — 1. LA DIFFUSION ÉLASTIQUE 
PAR L'ATOME. — La diffusion d’un photon implique 
un changement de sa quantité de mouvement 
vectorielle, le respect des lois de conservation étant 
assuré par les effets subis par le système quantique 
diffuseur, qui prend le recul. On peut classer les 
processus de diffusion selon l’état final dans lequel 
est laissé le système diffuseur 

La diffusion élastique par l’atome correspond à 
un état final où l’atome est resté intact et s’est 
contenté de prendre une certaine quantité de mou- 
vement. Les lois de conservation montrent alors 
que le photon perd très peu d’énergie dans la dif- 
fusion (moins de 20 eV pour des atomes lourds et 
des photons de 2 MeV). D’autre part, si on consi- 
dère ce type de diffusion pour un ensemble 
d’atomes, l’état final est un état où un seul atome 
a pris le recul. La diffusion élastique de l’atome est 
une diffusion inélastique pour le réseau cristallin 
dont il fait partie. Les probabilités de diffusion se 
composent comme des probabilités classiques très 
petites, la probabilité de diffusion du photon est la 
somme des probabilités de diffusion par chacun des 
atomes. ) 

Toutefois, si on ne peut pas négliger, devant 
l'énergie cédée à l’atome, son énergie de liaison dans 
le réseau, la liaison transforme le phénomène. Il 
peut arriver que l’ensemble du réseau subisse l'effet 
de recul. L'état où un atome seulement est mis en 
mouvement n’est alors qu’un état intermédiaire. 
L’amplitude de probabilité résultante est la somme 
des amplitudes de probabilité de diffusion par 
chacun des atomes, il y a interférence entre les 
atomes. Dans l’état final, le réseau est conservé, 
c’est une diffusion élastique pour le réseau. La perte 
d'énergie est encore plus petite que dans le cas 
précédent. 

Pour des rayons y de plus d’un MeV, la diffusion 
élastique par un réseau ne peut jouer de rôle appré- 
ciable que pour des angles de diffusion inférieurs 


Measurements are described for 2,62 MeV +y-rays from Th C” 


The cross sections of elastic scattering have 


au 1/100 de radian [3], où l’énergie communiquée à 
l’atome est du même ordre que l’énergie de liaison 
dans le réseau. Dans les expériences que nous 
décrivons ici, on peut négliger ce phénomène. 

La diffusion inélastique par l’atome correspond à 
un état final où un seul atome a subi le recul, mais 
où il a été modifié dans sa structure. Le photon a 
perdu une partie plus ou moins grande de son 
énergie et laisse l’atome ionisé, ou accompagné 
d’une paire d'électrons positif et négatif, ou bien 
seulement excité. La probabilité de diffusion est la 
somme des probabilités de diffusion relatives à 
chacun des états finals possibles. 

Pour des rayons y de l’ordre du MeV, la liaison 
des électrons dans l’atome est petite devant l’éner- 
gie communiquée à un électron du cortège par 
interaction directe, dès que l’angle de diffusion 
dépasse quelques degrés. La plupart des diffusions 
par l’atome sont inélastiques et correspondent à un 
état final où l’électron a pris beaucoup plus d’éner- 
gie que l’ion résiduel. C’est presque une diffusion 
élastique sur un électron libre. C’est l'effet Compton. 
La majorité des photons diffusés ont une énergie  v 
voisine de l’énergie donnée par la relation de 
Compton : 


(L.1) 


où À vw est l'énergie des rayons y incidents, 
mc? est l’énergie de masse de l’électron, 
0 est l’angle de diffusion. 


Néanmoins, la liaison de l’électron dans l’état 
initial introduit une dispersion de la diffusion iné- 
lastique et le spectre d’énergie des rayons y dif- 
fusés sous un angle déterminé s’étend de façon 
continue jusqu’à une énergie voisine de À w. Ce 
spectre d'énergie ne peut pas être tout à fait négligé 
devant la diffusion élastique, au voisinage de 
k vo [4]. 
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2. Les composantes de la diffusion élastique. — 
La probabilité de transition de l’état initial (sys- 
tème diffuseur et photon incident) vers l’état final 
(système diffuseur et photon diffusé) se calcule à 
partir des amplitudes de probabilité relatives à 
l’absorption et à l’émission. On considère des états 
intermédiaires, où le photon incident est absorbé 
ou bien le photon diffusé déjà émis ; l'amplitude de 
probabilité résultante est la somme des amplitudes 
de probabilité relatives à la transition de l’état 
initial vers l’état final en passant par chacun des 
états intermédiaires. Pour les états intermédiaires, 
ce ne sont pas les probabilités mais les amplitudes 
de probabilités qui s’additionnent. Il y a inter- 
férence entre les composantes d’une même dit- 
fusion. 

On peut classer les différentes composantes de la 
diffusion élastique selon le type d’état intermédiaire 
qu’elles font considérer. 


a) L'EFFET THOMSON Du NoYAU. — Le noyau 
seul subit une interaction pour passer de l’état 
initial à l’état intermédiaire et il ne prend que de 
l'énergie cinétique. On peut décrire ce processus 
comme un effet Compton sur le noyau, les électrons 
du cortège étant entraînés avec lui dans l’état 
final. Quand le photon a une énergie de l’ordre du 
MeV, très petite devant l'énergie au repos, Mc?, du 
noyau, l'amplitude de diffusion pour cet effet se 
réduit à : 


2 
ap = — 0,78.10—16 em À HE Los 0) 
pour la conservation de la polarition du photon 
, Z? 
ar = — 0,78.10—16 cm Pa, + cos 0) 


pour le changement de la polarisation du photon. 


où Z et À sont le nombre et la masse atomique, 
0 l’angle de diffusion. 


b) L'EFFET DE RÉSONANCE NUCLÉAIRE. — Dans 
l’état intermédiaire, le noyau se trouve dans un 
état excité. Ce processus ne doit jouer de rôle 
appréciable que si l'énergie du rayon y est très 
voisine de l’énergie d’un niveau d’excitation du 
noyau, circonstance tout à fait exceptionnelle. En 
effet, Levinger [5] a montré que les très hauts 
niveaux d’excitation des noyaux, de plus 
de 15 MeV, ne peuvent avoir une grande influence 
sur la diffusion des photons dont l'énergie est de 
l’ordre du MeV. 


c) L'EFFET DELBRÜCK OU EFFET DE POLARI- 
SATION-DU VIDE. — L'état intermédiaire est cons- 
titué par l’atome avec une paire d'électrons positif 
et négatif créée dans le champ électrique du noyau. 
Le calcul de l'amplitude de cet effet se réduit au 
calcul de la diffusion cohérente d’un photon par 
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création de paires virtuelles, dans le champ élec- 
trique du noyau assimilé à un champ coulombien 
extérieur. En outre, la sommation sur les états 
intermédiaires, s'étend à tous les états possibles 
constitués par une paire d'électrons dans le champ, 
y compris ceux où l’électron négatil se trouve dans 
un des niveaux d'énergie déjà occupé par un des 
électrons du cortège. Les termes calculés, ap et ay 
relatifs à la conservation et au changement de la 
polarisation du photon, correspondent à la dif- 
fusion par un noyau dépouillé de ses électrons. 
Pour un angle de diffusion 0 = 0, ap est cal- 
culé [6] et ap est nul pour des raisons de symétrie. 


d) L'EFFET RAYLEIGH. — L'état intermédiaire 
diffère de l’état initial de l’atome en ce qu’un 
électron du cortège est passé sur un niveau excité 
ou non lié. Le calcul des amplitudes a et ag corres- 
pondant à la conservation et au changement de la 
polarisation, se réduit au calcul de la diffusion 
oee des électrons liés dans un champ 
coulombien extérieur. On étend la sommation sur 
les états intermédiaires aux états d'énergie néga- 
tive des électrons, pourtant interdits par le prin- 
cipe d’exclusion. Ceci provient du fait que la pré- 
sence des électrons du cortège interdit la création 
d’une paire avec l’électron négatif sur un niveau 
déjà occupé ; les termes calculés en tenant compte 
des énergies négatives compensent les termes cal- 
culés en trop dans l'effet de polarisation du vide. 
Pour les amplitudes de diffusion a et ar de l’effet 
Rayleigh, on dispose de formules valables pour des 
photons d’énergie inférieure à 0,5 MeV [7] et de 
calculs numériques jusqu’à l’énergie de 1,3 MeV [8]. 

Si l’on peut effectivement négliger la résonance 
nucléaire, l'amplitude a de la diffusion élastique est 
la somme des amplitudes de diffusion par les effets 
Thomson, Delbruck et Rayleigh, puisque ces effets 
ne diffèrent que par leurs états intermédiaires. 


& = 4r + ap + ar \ 
(conservation de la polarisation). | 
" , , , P (1,3) 
a = ar + ap + am | 
(changement de la polarisation). 


Les amplitudes a et a’ correspondent à des états 
finals différents puisque le photon diffusé n’a pas la 
même polarisation. Néanmoins, au cours de nos 
expériences, nous mesurons la somme des nombres 
de’ photons relatifs aux deux polarisations. La 
section efficace mesurée sera par conséquent : 


Gi [2/22 af?, (I.4) 

3. Principe de l’expérience. — Le principe de 
l'expérience découle de ce qui précède. La compo- 
sante due à la polarisation du vide se manifeste par 
son influence sur la diffusion élastique. Elle est 


mise en évidence par une différence significative que 


660 


l’on pourrait constater entre la valeur mesurée de la 
section efficace de diffusion élastique et la valeur 
théorique de celle-ci, calculée en tenant compte des 
effets Thomson et Rayleigh seulement. A cause des 
phénomènes d’interférence, cette différence peut 
être négative ou positive, selon les signes CEA na 
tudes de diffusion des composantes. 

L'expérience consiste à mesurer les sections elfi- 
caces de diffusion élastique par l’atome, afin qu’on 
puisse comparer aux prévisions théoriques les 
valeurs obtenues expérimentalement. 

En disposant, à la vue d’une source et d’un détec- 
teur à scmtillation, un diffuseur de densité p, de 
poids atomique À, de coefficient d'absorption y 
pour les photons de l’énergie considérée, on enre- 
gistre un taux de comptage dû à la diffusion 
élastique par l’atome : 


n=N£ ft) JE 


À 
r? R? 


e—utæ+ x) AW (L.5) 


où N est le nombre d’Avogadro ; d V un élément de 
volume du diffuseur ; 0 l’angle de diffusion d’un 
photon émis par là source, diffusé par d V et traver- 
sant le scintillateur ; c(0) la section efficace 
en cm?/st, correspondant à cet angle ; J l'intensité 
de la source en nombre de photons émis, par sté- 
radiant, dans la direction de dV ; S la « section 
efficace » du système détecteur, vu de l’élément de 
volume ; r et À les distances de d V au centre-de la 
source et au centre du scintillateur respectivement ; 
x la distance que les photons parcourent à l’inté- 
rieur du diffuseur, entre la source et l’élément dV 
et X celle qu'ils ont à parcourir dans le diffuseur 
avant d'atteindre le détecteur (fig. 1). 


Détecteur 


AE RE be 


L'intégrale de (1.5) est étendue à tous les points 
du diffuseur. 

Certes, les événements enregistrés sont dus à des 
photons émis de divers points de la source, diffusés 
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en différents points du diffuseur et détectés en dif- 
férents points du scintillateur. Ils correspondent à 
des angles de diffusion 0 qui ont des valeurs 
diverses, comprises entre une valeur minimum Owin 
et une valeur maximum 0,3. Néanmoins, il existe 
un angle 0,, compris entre Onin 64 Omax, tel que 


op f JS 


—p(t+ X) 
À) rR° 2 


n = 6(0,) N< (1.6) 


Par la mesure de tous les termes qui inter- 
viennent dans cette formule, on détermine la 
section efficace c(0,), correspondant à un angle de 
diffusion 0, que l’on sait compris entre les deux 
valeurs extrêmes Oin 0t Omax, Ces deux valeurs étant 
déterminées par la géométrie du dispositif expéri- 
mental. 

Nous décrivons ici nos mesures, pour différents 
angles 0,, de la section efficace de diffusion élas- 
tique, dans le plomb, l'uranium et le bismuth, des 
photons de 2,62 MeV émis par le Th C”. 


# 
* * 


IT. Le dispositif. — 1. LA GÉOMÉTRIE. — Pour 
ces expériences, nous disposions d’une source de 
quelque 150 me Msth (Ra 228), qui contenait envi- 
ron 50 me de Th C” (T1 208). À 2,62 MeV, les 
sections efficaces de diffusion élastique sont si 
petites que, pour obtenir des taux de comptage 
suffisants, il a fallu faire, du dispositif, une étude 
approfondie [91], dont nous donnons la conclusion 
ici. 

La meilleure géométrie est celle qui, pour un 
domaine d'incertitude angulaire AO = O5ax — Omin 
déterminé, rend maximum le taux de comptage » 
(formule 1.6). Il faut que le diffuseur épouse aussi 
bien que possible la forme d’une portion d’un lieu 
géométrique, celui des points qui diffusent, vers le 
centre D du détecteur, les photons issus du centre O 
de la source sous un angle déterminé. On peut 
trouver des expressions approchées pour A0 et 
pour ñ [9] et. en résolvant un problème d’extrema 
liés, on s aperçoit que le diffuseur doit se réduire à 
peu près à un tôle roulée en forme de cylindre 
d’axe OD; placé à égale distance des centres O de la 
source et D du détecteur. L’épaisseur de la tôle 
doit être du même ordre de grandeur que le libre 
parcours moyen {/x des photons ; la hauteur du 
cylindre, égale au tiers environ de la distance OD. 
C'est la géométrie que nous avons adoptée, la 


figure 2 en donne une application pour l’étude de la 


diffusion sous 550. 

D'autre part, nous avons étudié le problème de 
l’élimination des rayonnements parasites de 
2,62 MeV, que le détecteur ne peut distinguer du 
rayonnement diffusé élastiquement et dont il faut 
se protéger au moyen d’écrans absorbants. Nous 
placions entre la source et le détecteur, un bouchon 


| 


| 


ï 
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(fig. 2) de hauteur telle que le rayonnement direct 
ne produisait plus, dans un canal de 60 keV autour 
de 2,62 MeV, qu'un taux de comptage inférieur 
au 1/10 du taux de comptage produit par les rayons 
cosmiques laissant la même énergie dans le scintil- 
lateur, Nous avons étudié la disposition et la forme 


Château protecteur 
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du faisceau 
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des pièces, de manière à réduire au minimum les 
diffusions élastiques parasites par d’autres objets 
que le diffuseur. Nous avons vérifié finalement que 
la présence de la source n’augmentait pas de façon 
appréciable le bruit de fond mesuré en l’absence de 
diffuseur. 

Pour les angles de diffusion de 25° à 709, où la 
précision sur l’angle de diffusion est très impor- 
tante, nous employions un bouchon de plomb 
de 35 cm de hauteur et de 14 em de diamètre et il 
y avait de 50 cm à 1 m entre la source et le détec- 
teur. Pour 1000 et 1250, les sections efficaces, très 
petites, nous ont obligés à réduire toutes les dis- 
Lances en jeu dans les expériences et il a fallu 
employer un bouchon d'uranium de 21 em de 
hauteur et 9 cm de diamètre, ce qui nous per- 
mettait de mettre seulement 30 cm entre la source 
et le détecteur. Avec le bouchon d'uranium nous ne 
constations pas d'augmentation sensible du bruit 
- de fond dans la région de 2,62 MeV ; ceci tient au 
fait que le rayonnement de l'uranium se compose 
entièrement de photons de moins de 2,62 MeV. 


Nous avons eu enfin à nous protéger contre un. 


rayonnement très gênant de 2,62 MeV, aisément 
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reconnaissable dans le bruit de fond mesuré en 
l’absence de source. Nous avons pu démontrer que 
son origine résidait entièrement dans la teneur 
en Th C” des murs du laboratoire (teneur corres- 
pondant à 12 g du Th par tonne de béton). Afin de 
réduire son influence, 1l a fallu disposer autour du 
détecteur, dans toutes les directions que n’em- 
pruntait pas le rayonnement diffusé intéressant, 
des écrans de plomb différents pour chaque expé- 
rience, et nous avons réussi de cette manière à 
réduire l’effet de la radioactivité naturelle dans un 
rapport 1/2 à 1/3. 


2. LE DÉTECTEUR. — Pendant la mesure de n 
(qui permet de déterminer la section efficace par la 
formule 1.6), le détecteur reçoit des photons para- 
sites provenant de la diffusion inélastique par 
l’atome, de l’annihilation des positrons formés dans 
le diffuseur, annihilation en vol ou au repos, de la 
diffusion élastique ou nôn, des photons de toutes 
énergies émis par la source de Msth et du rayon- 
nement de freinage des électrons émis et ralentis 
dans le diffuseur. Tous ces processus produisent des 
photons très nombreux d’énergie nettement infé- 
rieure à 2,62 MeV et d’autres, moins nombreux, 
dont le spectre décroît quand l’énergie croît mais 
ne devient pas tout à fait négligeable devant la 
diffusion élastique dans la région de 2,62 MeV 


(fig. 3). 
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Notre détecteur à scintillation, était un cristal 
de Nal(T1) de 4 em sur 5 em monté sur un photo- 
multiplicateur EMI 6260. On analysait les impul- 
sions, mises en forme à 410$ s, dans une voie lente 
aboutissant à un sélecteur à cinq canaux, réglé 
comme le montre la figure 3. On a représenté éga- 
lement sur la figure 3 le pic photoélectrique des 
photons de 2,62 MeV dont la largeur à mi-hauteur 
était de 5 %. 

Afin de réduire le comptage dû aux empilements 
de petites impulsions, il a fallu ajouter une voie de 
sélection rapide, où les impulsions étaient mises en 
forrne à 10-7 s, aboutissant à un discriminateur 
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dont le seuil correspondait à une amplitude 
de 2,2 MeV et commandant l’ouverture d’une porte 
sur la voie lente. Le schéma de principe du dispo- 
sitif d'analyse est représenté sur la figure 4. 
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Ce système de double sélection présentait les 
caractéristiques de la sélection rapide du point de 
vue de la résolution en temps et, en divisant par 4 
ou 5 le nombre des empilements, permettait d’éli- 
miner l'influence des très nombreux parasites 
d'énergie comprise entre 0,5 et 2 MeV. Il présentait 
néanmoins les caractéristiques de la voie lente du 
point de vue résolution en énergie et permettait 
l'élimination des parasites d’énergie proche de 
2,62 MeV ; en effet, l'introduction de la sélection 
rapide ne modifiait ni la hauteur ni la forme du 
pic photoélectrique des rayons y du Th C". 

Mais, pendant les expériences nous avons trouvé 
que le plus grand nombre des impulsions était 
d'amplitude très inférieure à 0,5 MeV: Trop petites, 
elles pouvaient rarement, dans la voie lente, 
s’empiler jusqu’à 2,6 MeV ; par contre, il arrivait 
qu’elles se superposent à des impulsions intéres- 
santes et en faussent l’analyse par le sélecteur. Ce 
phénomène agissait en élargissant le pic photo- 
électrique des rayons y du Th C”. Nous placions 
au-dessus du cristal scintillateur des filtres de 
plomb, de 1 à 3 cm d’épaisseur, qui réduisaient le 
rayonnement de moins de 0,5 MeV dans un rapport 
1/100 à 3/100 mais laissaient passer 1/4 à 2/3 du 
rayonnement diffusé élastiquement. Nous avons 
vérifié que le spectre d’une source de Th C” n’était 
plus déformé notablement par la présence du rayon- 
nement diffusé. 


III. La méthode. — 1. DÉTERMINATION DE LA 
SECTION EFFICACE. — Pour appliquer la for- 
mule (1.6) au cas de la géométrie décrite ci-dessus, 
nous avons étudié l’anisotropie de la source et du 
détecteur, puis développé en série la fonction à 
intégrer dans (1.6). On fait une erreur toujours 
inférieure à 1,3 % en supposant, sur une généra- 
trice du cylindre diffuseur, cette fonction cons- 
tante et égale à sa valeur pour le point situé à 
mi-hauteur. En outre, l'intensité J durayonnement 
de la source possédait la même symétrie cylindrique 
que le dispositif, tout eomme la fonction 

1 
r2 R? 

évidemment. Par contre nous n’avons pas pu consi- 
dérer que la grandeur #, liée au détecteur, füt symé- 
trique, probablement à cause d’un manque d’uni- 
formité dans l’épaisseur du plomb qui filtrait le 
rayonnement diffusé. Il a fallu mesurer la gran- 
deur S en suspendant une source auxiliaire, d’une 
fraction de millicurie de Th C”, à l’intérieur du 
cylindre diffuseur de telle manière que son centre 
d'activité se trouve sur le plan de symétrie hori- 
zontal du diffuseur. On mesurait le taux de comp- 
tage n' produit par la source auxiliaire et on cal- 
culait la moyenne n’ relative à différentes positions 
de cette source sur la circonférence du diffuseur. 

Le calcul complet [9] montre alors que la section 
efficace 6(0,) peut s'exprimer par : 


eu +X 


c(0,) == f8— (IIL.1) 
n 
où 
_ A 
FE 


FSURES HO 
£ = 7, 1Rk sin 0, sin Sr 


J'et J’ sont les nombres de photons de 2,62 MeV 
émis par unité d’angle solide, respectivement par la 
source de l’expérience et par la source auxiliaire. 
À, pet u sont le poids atomique, la densité et le 
coefficient linéaire d'absorption du corps diffuseur ; 
N, le nombre d’'Avogadro. ZL est la distance du 
centre de la source au centre du scintillateur ; k, la 
hauteur du diffuseur et Z tg 0,/2 le diamètre 
moyen qu’on lui attribue, compte tenu du fait que 
les diffusions enregistrées ont lieu, en moyenne, à 


PA 
une certaine distance Su LE È = 1) de sa 


surface intérieure. Que 


e est l’épaisseur du diffuseur ; k = 1 — e 
un coefficient toujours voisin de 1. 


k', encore plus voisin de 1, corrige une erreur 


faite lors de la mesure de S décrite ci-dessus, où la 
source est suspendue à l’intérieur et non à la place 
du diffuseur, 
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2. L'ANGLE EFFICACE DE DIFFUSION 0,5. — Dans 
nos expériences, les angles de diffusion extrêmes 
Gnnin Et Onax définissaient un domaine d’incertitude 
assez restreint pour l’angle 6, qui correspondait 
à la mesure. On peut exprimer la fonction (0) par 
son développement limité : 


6(9) = 50) + (0 — 05) 5 (Ou). (111.2) 
Les formules (1.5) et (1.6) et (111.2) permettent 


alors de définir 0, comme une moyenne des angles 
de diffusion 


fo se ute+ x) dV—0. (IL) 

En toute rigueur, cette moyenne devrait 
s'étendre à tous les angles définis par chaque point 
de la source, chaque point du diffuseur et chaque 
point du détecteur. Mais, en comparant les formes 
de la source, du détecteur et du diffuseur (fig. 2), on 
s'aperçoit facilement que l’on peut considérer la 
source et le détecteur concentrés en leurs centres O 
et D. En faisant les mêmes approximations qu’au 
paragraphe précédent, on trouve une expression 
approchée 


80 = 3 0 + 3 0% 
(III.4) 
où lg 2e =" à f1 er 


La valeur de 6 déterminée par (111.4) est la 
valeur nominale de l’angle de diffusion, celle que 
l’on fait correspondre à la mesure de section efti- 
cace. Dans le domaine On — Omar elle est assez 
voisine de la valeur exacte, puisqu'on a tenu 
compte du développement en série de (0) jusqu’au 
terme de degré 1. Si l’on introduit un terme du 
2e ordre dans (111.2), on s'aperçoit que 0, est en 
réalité supérieure à la valeur déterminée par (111.4) 
si d?c/d0? a le même signe que do/d0’ qu’elle lui 
est inférieure dans le cas contraire. 


3. LES MESURES. — Pour mesurer le rap- 
port J’JJ, nous avons exposé, à la vue du détec- 
teur dépouillé de ses filtres, les deux sources uti- 
lisées dans l'expérience, après les avoir placées 
dans les mêmes récipients et orientées de même 
façon que dans nos mesures. Il a fallu attacher une 
valeur diflérente de J’/J à chaque mesure de 
section eflicace. Pour la détermination de ce rap- 
port, on considérait seulement les taux de comp- 
tage au sommet de pic photoélectrique des rayons 
gamma de 2,62 MeV, cela afin d'éviter des erreurs 
dues à l’émission de photons d’énergie inférieure. 
Pour faciliter les mesures, nous avons comparé 
chaque source à une source d’activité intermédiaire. 

Nous avons mesuré J’/J à trois reprises, à cause 
de l’évolution du RdTh dans le Msth des sources 


qui étaient d’âges différents et inconnus. Par inter- 
polation, nous pouvions obtenir le rapport J'/J 
relatif à chaque expérience. L'erreur sur ce rapport 
était toujours inférieure à 5 %. 

Dans chacun des 5 canaux du sélecteur, corres- 
pondant à une « section efficace » S déterminée de 
l'appareil détecteur, les rapports n/n' devaient 
avoir la même valeur, puisque ces taux de comp- 
tage étaient dus à des photons de même énergie. 

La convergence des mesures de ces rapports dans 
les différents canaux, nous permettait de contrôler 
la stabilité du dispositif électronique et l’élimi- 
nation des parasites. Dans la formule (111.1), nous 
introduisons la moyenne pondérée des valeurs 
de nn’ trouvées dans les canaux 2, 3 et 4 de notre 
sélecteur (fig. 3), avec un poids 1 pour les 28 et 
4e canaux et un poids de 2 pour le 3. 

Nous n’avons jamais pu réduire à une valeur 
négligeable le bruit de fond dû à la radioactivité 
naturelle et aux rayons cosmiques. Il a fallu donc, 
en escamotant la source, mesurer ce bruit de fond, 
n», et le soustraire du taux de comptage 
ñ, =n + n, mesuré en sa présence. On déter- 
minait n. et r, au moyen d’une dizaine de mesures 
alternées, de deux heures environ chacune. Ainsi 
l'erreur statistique sur An était négigeable. A 
chaque interruption de comptage nous effectuions 
un étalonnage au moyen de la source auxiliaire. 


IV. Observations et résultats. — 1. ALLURE 
GÉNÉRALE DES SPEGTRES. — Le rayonnement dif- 
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Fi. 5. — Spectres sans double sélection pour 559, 


fusé produisait, dans le spectrographe à scintil- 
lation, un grand nombre d’impulsions dont les 
amplitudes correspondaient à des énergies infé- 


| 
| 
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rieures à 2,6 MeV. À titre d'exemple, nous repré- 
sentons, en trait plein, sur la figure 5, le spectre des 
impulsions produites par les photons diffusés par le 
plomb ; en trait pointillé, le spectre obtenu dans 
les mêmes conditions avec un diffuseur d’atomes 
légers, constitué par une cuve de fer remplie 
d’eau [1]. Les spectres des impulsions d’amplitude 
inférieure à 2,6 MeV s'expliquent entièrement par 
lPaction des parasites de moins de 2,6 MeV dont 
nous avons dressé la liste au paragraphe 2 du 
chapitre IT. 

Dans la partie inférieure des spectres, où la fré- 
quence des impulsions était de l’ordre de 10% s-1 
dans un canal de 60 keV, les taux de comptage 
relatifs aux atomes lourds et aux atomes légers 
étaient approximativement les mêmes ; ils se diffé- 
renciaient nettement ensuite, dans la partie supé- 
rieure, mais là le nombre de impulsions était beau- 
coup plus petit. On peut donc penser que le phéno- 
mène d’empilements était du même ordre de gran- 
deur dans les deux cas. 

Autour de 2,62 MeV, les spectres relatifs à la 
diffusion par les atomes légers ne présentaient pas 
de pic, une fois déduite la contribution du rayon- 
nement cosmique et de la radioactivité naturelle. 
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o 
x 
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2,6 ; spy 
Amplitude des impulsions en MeV 
Fra. 6. — Spectres obtenus pour l’angle de diffusion de 45°, 


La figure 6 montre, en trait plein, un de ces 
spectres pour l’angle de diffusion de 450, en trait 
pointillé, le bruit de fond, que l’on avait soustrait 
auparavant. Dans les spectres relatifs aux atomes 
légers, 1l est vraisemblable que le phénomène 
d’empilement était prépondérant. 

Sauf pour le plomb sous 250, les taux de comp- 
tage dus aux photons diffusés par les atomes 
lourds, étaient bien plus élevés que les taux de 
comptage relatifs aux atomes légers, et les spectres 
présentaient un pic à l’abscisse de 2,62 MeV, même 
après soustraction du bruit de fond. Mais le pic 
était moins prononcé que dans le spectre étalon. 
Sur la figure 7, nous avons représenté des spectres 
relevés lors de l’étude de la diffusion sous 459 par 
les atomes d'uranium. 


JOURNAL DE PHYSIQUE 


No 2? 


Pour le plomb sous 250, on ne trouve pas de pic 
dans le spectre des impulsions produites par les 
photons diffusés et les taux de comptage ne sont 
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Fi. 8. — Spectres obtenus pour l’angle de diffusion de 259, 


que 2 à 3 fois supérieurs aux taux de comptage 
correspondant au diffuseur d’atomes légers (fig. 8). 
Il est vraisemblable qu'ici le système de double 
sélection décrit au chapitre 11, ne suffit pas à rendre 
négligeable le phénomène d’empilement. 


2. ANALYSE DES SPECTRES ET PRÉCISION DES 
MESURES. — Dans le cas des atomes légers, 
l'absence de pie montre que la diffusion élastique 
participe aux comptages de façon négligeable seu- 
lement. Dans le cas des atomes lourds, seule la 
diffusion élastique de l’atome peut expliquer la 
présence d’un pic à l’abscisse de 2,62 MeV. On en 
déduit que la section efficace de diffusion élastique 
croît plus vite que le coefficient massique d’absorp- 
tion u/6 qui intervient dans la formule (111.1). 

Les spectres relatifs à la diffusion par les atomes 
lourds ne sont pas tout à fait proportionnels aux 
spectres relevés pendant les étalonnages avec la 


danses 
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source de Th C”. Ceci prouve l’existence de phéno- 
mènes perturbateurs. Le taux de comptage dans le 
premier canal du sélecteur est supérieur au taux de 
comptage dans le second (fig. 7), ce qui établit la 
présence de parasites, rendant inutilisable la 
mesure du rapport n/n' dans le premier canal. 
Cependant l’asymétrie du pic permet d’évaluer 
pour les trois canaux médians, la contribution de la 
diffusion inélastique, dont le spectre d’impulsions 
doit nécessairement décroître rapidement en fonc- 
tion de l'amplitude. D’autre part, le taux de comp- 
tage produit par les atomes légers représente une 
borne supérieure pour la contribution des empi- 
lements, ce que l’on peut évaluer également par la 
variation du taux de comptage en fonction de 
l'épaisseur de l’écran de plomb disposé sur le détec- 
teur (voir chap. IT fin du paragraphe 2). 

Nous avons pu établir que la contribution des 
parasites était faible et ne pouvait entièrement 
rendre compte des défauts de proportionnalité dans 
les 2e, 3e et 4€ canaux. D’autre part, nous avons 
constaté des phénomènes liés au fonctionnement 
du dispositif électronique, qui pouvaient expliquer 
un certain élargissement du pic dans le spectre de 
la diffusion [9]. Nous avons pu évaluer directement 
l'importance de ces phénomènes et nous avons 
trouvé qu’ils ne pouvaient être entièrement res- 
ponsables des défauts de proportionnalité. Il est 
probable que les défauts de proportionnalité résul- 
tent de l’action simultanée des parasites et de ces 
phénomènes liés au dispositif électronique. 

Tous les phénomènes perturbateurs ont ten- 
dance à effacer le pic et à rendre, dans les 22 et 
4e canaux, le rapport n/n' plus grand que dans 
le 3 canal. À partir des écarts entre les valeurs de 
ce rapport, mesurées dans les trois canaux médians, 
on déterminait une limite supérieure de l'erreur 
dont ces phénomènes pouvaient affecter nos 
mesures. Néanmoins, les parasites avaient ten- 
dance à augmenter les taux de comptage alors que 
les phénomènes liés au fonctionnement du dispo- 
sitif électronique, au contraire, les diminuaient 
dans les trois canaux médians. Ceux-là faisaient 
surestimer les sections efficaces tandis que ceux-ci 
les faisaient sous-estimer. Puisque nous subissions 
à la fois l’effet de ces deux causes d’erreur, il est 
vraisemblable que les erreurs attribués à nos 
mesures sont nettement supérieures aux erreurs 
réelles dont elles sont entachées. 

Les erreurs dues aux fluctuations statistiques 
étaient négligeables, l'erreur sur J'/J était peu 
importante et l'erreur attribuée aux sections effi- 
caces tient essentiellement aux défauts de propor- 
tionnalité. 

Lors de l’étude de la diffusion sous 250 par le 
plomb, les défauts de proportionnalité étaient tels 
que l’on ne pouvait pas déterminer la participation 
de la diffusion élastique. Nous n’avons pu donner 
à la section efficace qu’une borne supérieure. 


4. LES RÉSULTATS. — Pour le plomb, le bismuth 
et l’uranium nous avons trouvé les sections effi- 
caces indiquées dans le tableau ci-dessous. 6 est la 
section efficace de diffusion élastique en em? par 
stéradiant, pour les rayons gamma de 2,62 MeV 
0, est l’angle efficace de diffusion, déterminé selon 
la méthode décrite au paragraphe 2 du chapitre IIT. 
Omin Et Omax Sont les angles extrêmes sous lesquels 
on enregistre des diffusions pendant la mesure de 6. 


TABLEAU 
ERREUR 
NATURE DU 
2e NS DER 6 EN CM? REDATINE MO AO in LOS 
SUR 6 
RON 260 220 — 300 
Sel0rs2e . 20% 39 300 — 430 
1,6 10—28", 20 % 460 400 570 
Plomb 7 1102 90 580 520 680 
Do 0e 200 710 , 620 — 820 
LRO 200% 41019 930 1130 
HS 0) 20995 1269 1160 — 1350 
ALTAQEÉE 207 &60 4092— 570 
: HRElO 22 2007 739  640— 840 
PouiR 3,8 10-%. 200, 14030 94041140 
5,3 10—2 410.9, 1260 1160 — 1350 
10e 30 % 2500940990 
4,8 107 ?8 AXDAGE 340 (290 420 
2 1028 20% 450 400 540 
Uranium 159 10527 282007 569 500 — 640 
S6402 20 719 620 — 820 
HOMO? AD 1019 939 — 1130 
LA A Nm) 20202 1250 1169 — 1350 


x Résultats pour Pb 


Rayleigh+Thomson pour Pb 
À d'après Bethe 
\ 


d'après Brown et Mayers 


Section efficace (cm? par stéradian) 


————— 
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Fic. 9. — Sections efficaces de diffusion élastique des 
rayons y de 2,62 MeV dans Pb et Bi. 
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Ces valeurs sont portées sur les figures 9 et 10. 


Rayleigh + Thomson 
D d'apres Bethe 


Section efficace (cm? par stéradian ) 
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FiG. 10. — Sections efficaces de diffusion élastique 
des rayons y de 2,62 MeV dans l’uranium. 


Les expériences effectuées avec les diffuseurs 
d’atomes légers permettent d’affirmer que la sec- 
tion efficace de diffusion élastique croît plus rapi- 
dement que le nombre atomique Z. 


V. Discussion. — 1. VALIDITÉ DES MESURES. — 
Dans chacune des expériences correspondant à une 
détermination de section efficace, il y avait un pic 
dans le spectre des impulsions produites par les 
photons diffusés et les écarts entre les différentes 
valeurs du rapport r/n'ne dépassaient jamais 15 % 
pour les trois canaux médians du sélecteur. Ceci ne 
peut s'expliquer que par l’effet de la diffusion élas- 
tique de l’atome et les mesures de celle-ci ne peu- 
vent être entachées d’erreurs supérieures à celles 
que nous leur avons attribuées dans le tableau. 
Par conséquent, contrairement à l'hypothèse 
avancée par Brown [8] pour expliquer, à 1,33 MeV, 
les résultats obtenus par d’autres expérimen- 
tateurs, il est exclu qu’un parasite tel que la dif- 
fusion inélastique des électrons liés ait faussé nos 
déterminations de sections efficaces dans une pro- 
portion plus grande que celle que nous avons 
évaluée. 

Avec notre dispositif électronique, les taux de 
comptage produits par la diffusion sur les atomes 
légers étaient faibles et permettaient d’estimer que 
le nombre des empilements étaient négligeables 
dans les spectres relatifs aux atomes lourds. Nous 
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n’avons ainsi pas eu à recourir au procédé de sous- 
traction de ces deux spectres, pour éliminer les 
effets des empilements, procédé peu légitime dans 
les expériences décrites ici où les spectres des 
impulsions de petite amplitude ne sont pas abso- 
lument identiques (fig. 5) dans le cas des atomes 
lourds et dans le cas des atomes légers. 


2. LA RÉSONANCE NUCLÉAIRE. — Parmi les 
composantes de la diffusion, élastique (chapitre I, 
paragraphe 2), la diffusion par résonance nucléaire 
a des propriétés caractéristiques : Elle prend de 
grandes proportions lorsque l’énergie du rayon 
gamma est, dans le système du centre de masse, 
voisine d’un niveau d’excitation du noyau diflu- 
seur, et diminue très rapidement dès que l’énergie 
du rayon gamma s’éoigne de l’énergie du niveau. 
Elle se réduit alors à la contribution, négligeable 
dans notre cas [5], des niveaux très éloignés, qui se 
manifestent dans les réactions photonucléaires à 
partir de 15 MeV. 

Dans une source de Th C”, l'émission du photon 
de 2,62 MeV se produit au moment du passage à 
l’état fondamental d’un noyau de Pb 208. 2,62 MeV 
est donc un des niveaux excités du Pb 208, dont 
l'abondance est de 52 % dans le plomb naturel. La 
différence d’énergie entre le rayon gamma et le 
niveau excité correspond juste à l’énergie emportée 
par le recul du noyau émetteur (près de 20 eV). 
De telles circonstances n’existent pas pour le 
Bi 209, qui constitue entièrement le bismuth 
naturel. On peut donc s’attendre à observer, dans 
la diffusion par le plomb, un effet de résonance 
nucléaire beaucoup plus grand que dans la dif- 
fusion par le bismuth. 

Les valeurs obtenues pour les sections efficaces 
du bismuth (Z — 83) et de l’uranium (Z = 92) 
donnent une idée de la dépendance en fonction du 
nombre quantique Z des composantes Thomson, 
Rayleigh et Delbrück, et permettent d’évaluer leur 
variation, du plomb (Z — 82) au bismuth (Z = 83) 
Cette variation ne serait jamais supérieure à 10 %,. 
Or, les différences entre les sections efficaces que 
nous obtenons pour le plomb et le bismuth ne sont 
pas significatives. Il s’ensuit que la résonance 
nucléaire joue un rôle négligeable dans la diffusion 
de nos rayons gamma par le Pb 208. 

Ainsi, pour une différence d'énergie de 20 eV 
entre le rayon gamma incident et le niveau d’exei- 
tation, la section efficace de diffusion par résonance 
nucléaire est plus petite que 2.10—*$ cm?/st, même 
si les déphasages ont les valeurs les plus défa- 
vorables à une mise en évidence de cet effet. Ces 
expériences confirment les prévisions selon les- 
quelles il serait impossible de mettre en évidence la 
diffusion de résonance nucléaire, pour une diffé- 
rence d'énergie de plus que quelques eV, prévisions 
confirmées d’ailleurs par les expériences de 
Moon [10]. 
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3. LES CALCULS DE L'EFFET RAYLEIGH. — 
Puisque la diffusion par résonance nucléaire est 
hors de cause, l'effet Delbrück serait effectivement 
mis en évidence par une différence significative 
constatée entre les sections efficaces mesurées et les 
sections efficaces calculées à partir des seuls effets 
Thomson du noyau (formule 1.2) et Rayleigh des 
électrons-liés (voir chapitre [, paragraphe 2). Or, 
pour l'effet Rayleigh, nous disposons d’une formule 
établie par Franz [7] pour des photons d’énergie 
petite devant l’énergie mc? de l’électron au repos. 


ar — — © f{1 + cos 0) \ 


conservation de la polarisation 
7 oi (V1) 


’ 


GR —f4 — cos 0) 
changement de polarisation. | 


où f est le facteur de structure dépendant du chan- 
gement de quantité de mouvement du photon, 
facteur égal à Z quand ce changement est égal à 
Zéro. 

En introduisant dans (V.1) le facteur de struc- 
ture, calculé par Bethe [7], des deux électrons X de 
l’atome, nous obtenons, pour la section efficace due 
aux effets Rayleigh et Thomson seulement, les 
courbes tracées en pointillé sur la figure 9 et en 
traits pleins sur la figure 10. Mais la validité de 
la formule (V.1) n’est pas démontrée dans notre 
cas, où l’énergie des photons est cinq fois supé- 
rieuree à l’énergie d’un électron au repos. Elle s’est 
déjà révélée fausse pour une énergie de 1,33 MeV 
des rayons gamma ; il en est très certainement de 
même à 2,02 MeV. 
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À partir des calculs tout à fait exacts de 
Brown [8] pour l'énergie de 1,30 MeV, on peut faire 
une extrapolation jusqu’à 2,62 MeV, en supposant 
que l’amplitude dépend, au facteur 1 Æ cos 0 près, 
du changement de quantité de mouvement seu- 
lement. On obtient alors la courbe en trait plein 
sur la figure 9. Il y aurait une différence signifi- 
cative qui pourrait s'interpréter par l'effet 
Delbrück, mais la légitimité de cette extrapolation 
est encore sujette à caution. 


4. CONCLUSION. — Pour mettre en évidence 
l’effet Delbiück dans la diffusion élastique, il faut 
des mesures de sections efficaces, qui sont décrites 
ici pour l’énergie de 2,62 MeV, et des calculs théo- 
riques, qui sont à continuer. 

On peut constater néanmoins, que les sections 
efficaces à 2,62 MeV sont de 10 à 100 fois inférieures 
à la valeur de la section efficace de diffusion 
Delbiück sous 00, calculée par Rohrlich et 
Gluckstern pour le plomb à cette énergie 
(0,9 1026 em?/st). Bien que l'effet Delbiück diminue 
avec l’angle de diffusion, il intervient proba- 
blement, s’il existe, de manière significative dans la 
diffusion élastique à 2,62 MeV. On peut s’attendre 
à ce que nos résultats puissent servir de test à 
lPexistence de la diffusion des photons par la polari- 
sation du vide. 
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DIFFUSION ÉLASTIQUE PAR LE PLOMB ET L'ÉTAIN DES RAYONS y 
DE 1,33 MeV ET DE 1,17 MeV DU Co, ET DE 1,12 MeV DU ‘Zn. 


Par MM. E. HARA, J. BANAIGS, P. EBERHARD, L. GOLDZAHL 


Laboratoire de Physique Atomique et Moléculaire du Collège de France 


et J. 


MEY, 


Centre d'Études Nucléaires de Saclay. 


Résumé, — On a mesuré la diffusion élastique des rayons gamma de 1,12; 1,17et 1,33 MeV, et 
comparé les résultats à de nouveaux calculs théoriques. On ne décèle aucune contribution de 
l’effet Delbrück, mais on met en évidence un effet de diffusion inélastique jusqu'ici négligé. 


Abstract. — Total elastic scattering is measured for 1.12, 1.17 and 1.33 MeV y-rays and com- 
pared with new theoretical calculations. No effect of Delbrück scattering is found, but an inelastic 
scattering appears, which we had neglected untill now. 


1. Introduction.— Dans un article précédent [1], 
ont été décrites des mesures de diffusion élastique 
des gamma de 1,33 MeV et de 1,17 MeV du Co. 

Les sections efficaces, pour 1,33 MeV et pour des 
angles de diffusion compris entre 35° et 900, furent 
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trouvées bien inférieures à celles calculées à partir 
du facteur de structure de Bethe, pour la somme 
des effets Rayleigh et Thomson, extrapolation 
hasardeuse d’une formule établie dans le cas des 
rayons X, De cette différence on pouvait conclure 
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Fic. 3. — Schéma de fonctionnement. 


en admettant la validité même très approximative 
de ce calcul à existence d’un nouveau processus de 
diffusion élastique qui avait été identifié à l'effet 
Delbruck. 
A 1,17 MeV, les sections efficaces, obtenues par 
une méthode indirecte, furent trouvées légèrement 
inférieures seulement aux prévisions du calcul, ce 
qui appuyait fortement cette interprétation. Etant 
. donnée pour celle-ci l’importance des mesures à 
1,17 MeV, les auteurs de ce premier article déci- 
dèrent de les reprendre et d’effectuer des mesures 
avec les gamma de 1,12 MeV du 6°Zn. En effec- 
tuant ces mesures, nous avons alors constaté que 
les résultats à 1,12 MeV, obtenus par la méthode 
directe, infirment la valeur des sections efficaces 


_ déduites des expériences à 1,17 MeV, ce qui montre. 


que la méthode indirecte n’est pas valable dans ce 
cas. 

D'ailleurs, de récents calculs de la somme des 
eflets Rayleigh et Thomson par Brown et Mayers, 
fondés sur des bases tout à fait satisfaisantes 
montrent un écart considérable avec ceux que nous 
avions utilisés pour interpréter les résultats à 
1,33 MeV, et renversent notre conclusion, comme 
on pourra le voir dans la discussion des résultats, 

Nous avons également repris les mesures à 
1,33 MeV, avec l’étain afin d’en améliorer la pré- 


cision et d'effectuer une mesure à 209, avec le 
plomb pour compléter les mesures en mesurant la 
section efficace à 200. 


2. Dispositif. — a) DiSsPOsiTIF MÉGANIQUE. — 
Les expériences pour des angles supérieurs ou 
égaux à 350 ont été effectuées à l’aide du dispositif 
mécanique utilisé dans l’article cité [1]. Ce dispo- 
sitif ne permettant pas les mésures pour des angles 
inférieurs à 359, nous en avons réalisé un nouveau 
qui permet les mesures entre 200 et 509 (fig. 1). 

Pour atteindre ce résultat, on a remplacé le 
bouchon de fermeture conique de la source par un 
chariot sur rails, supportant un écran de plomb, 
et qui s’escamote totalement lors de la mesure. De 
plus la manœuvre du dispositif est plus facile et 
plus rapide. 


b) DISPOSITIF ÉLECTRONIQUE. — Pour toutes ces 
expériences nous avons utilisé un nouveau dispo- 
sitif électronique permettant, grâce à une bonne 
résolution en temps, d'éviter les erreurs dues aux 
empilements des petites impulsions. Le schéma de 
fonctionnement de ce dispositif est indiqué sur les 
figures 2, 3a, 3b. 

Le photomultiplicateur est muni de deux voies 
de sortie, La première, à constantes de temps assez 
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longues, 1 us environ, à une mauvaise résolution en 
temps mais une bonne résolution en énergie. Elle 
est dirigée vers les sélecteurs et les échelles de comp- 
tage. En l’absence d’impulsions cette voie est main- 
tenue bloquée. Le déblocage de cette première voie 
est commandé par la seconde. Celle-ci, à constantes 
de temps plus courtes, n’a pas une bonne résolution 
en énergie, mais ayant une bonne résolution en 
temps, permet de mieux séparer les impulsions. 
Dans ces conditions, si une impulsion arrive seule à 
l'entrée du système (cas de la figure 3b), la voie 
lente est débloquée pendant 4 us environ, et l’im- 
pulsion est comptée. Par contre, s’il arrive deux 
impulsions en un court intervalle de temps, plus 
court que 3 us environ, la voie principale reste 
bloquée et les deux impulsions ne sont pas comptées 
(cas de la figure 3a). Les taux de comptage étant 
dans notre cas très faibles (quelques coups par 
minute) la perte de comptage ainsi occasionnée est 
négligeable. 

D'autre part, l'emploi d’un photomultiplicateur 
Du Mont 6292 et d’un cristal de INa de 4 cm de 
diamètre et de 5 cm de hauteur a permis d’obtenir 
une résolution de 5,7 % contre 7 % dans les expé- 
riences précédentes. 

Enfin les mesures ont été effectuées avec un 
nouveau sélecteur à 5 canaux, contrôlé par un 
sélecteur à 1 canal. 
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3. Mesures et résultats. -— Comme dans les expé- 
riences précédentes, nous avons tenu compte des 
empilements résiduels en soustrayant du spectre 
de la cible étudiée celui d’une cuve à eau équi- 
valente. Grâce au nouveau dispositif électronique 
Pimportance de ce dernier spectre est bien moins 
moindre qu'auparavant. 

Comme il est essentiel dans ces expériences de ne 
mesurer que la diffusion élastique, nous n’avons 
utilisé que la partie du spectre approximativement 
proportionnelle au spectre directement obtenu avec 
une source auxiliaire, le défaut de proportionnalité 
fournissant la plus grande partie de l'erreur attri- 
buée à la mesure ; en particulier la contribution 


de la diffusion inélastique y est incluse. Avec le Co - 


on ne retrouve pas dans le spectre diffusé le pic 
de 1,33 MeV et la seule partie du spectre propor- 
tionnelle au spectre de la source auxiliaire est celle 
située sur le sommet du pic et au delà. Autour 
de 1,17 MeV, le pic que l’on retrouve n’est propor- 
tionnel au spectre direct en aucune partie, ce qui 
enlève beaucoup de sûreté à la méthode indirecte. 
Par contre, avec le $°Zn, nous retrouvons dans le 
spectre diffusé le pic de 1,12 MeV, pour des angles 
supérieurs à 500; pour les angles inférieurs, on 
trouve une partie proportionnelle, comme à 
1,33 MeV, qui permet la mesure de la section 
efficace. 


TABLEAU 1 


SECTIONS EFFICACES DE DIFFUSION ÉLASTIQUE 6 MESURÉES DE 1,12 A 1,33 MeV 


ANGLE DE DIFFUSION 200 390 500 700 90° 4150 

os en mby/ster 43 + 30% 3,81 + 20 % 0,580 + 20 % 0,439 + 20 9, 0,145 + 20 9% 0,104 120007 
pour Pb à 1,33 MeV (179-230) (309-400) (479-530) (659-750) (859-950) (1089-1220) 

sen mb/ster k,5 + 50% : 0,41,€ 509% 0:049/E 30 % ‘0,011 + 20 % 0,0064 + 20 %  0,0095 2 2007 
pour Sn à 1,33 MeV (179-230) (350-410) (479-530) (689-740) (869-960) (1189-1220) 

o en mby/ster 81.4 30,9%65:.8,2:4,209 1,246 20 0/2 00/36: +20 10,30 2020007 000 ORAN 
pour Pb à 1,12 MeV (179-230) (299-440) (479-530) (659-750) (849-960) (1030-1280) 


Les résultats pour le plomb et l’étain à 1,33 MeV 
et le plomb à 1,12 MeV, sont indiqués dans le 
tableau I. 

Pour le plomb à 1,33 MeV et les angles compris 
entre 350 à 1159, ce sont les résultats déjà indiqués 
par Goldzahl et Eberhard [1]. Les valeurs sont 
reportées sur les figures 6, 7, 8. 

Les résultats à 1,17 MeV pour le plomb, sont 
indiqués dans le tableau ci-dessous où nous repor- 


TABLEAU 2 


RAPPORT DES SECTIONS EFFICACES 
A 1,17 ET 1,33 MeV Pour Pb 


ANGLE DE DIFFUSION 509 700 4159 
o 4 
Ar 6 + 2 6- + 2 2,1 + 0,3 
S1:33 


tons le rapport des sections eflicaces à 1,17 et 
1,33 MeV (tableau IT). 


4. Calcul de l’effet Rayleigh et Thomson. — 
Tous les résultats expérimentaux sont comparés 
aux résultats théoriques récemment calculés à 
Birmingham par Brown et Mayers [2]. Ceux-ci ont 
effectué des calculs exacts de la diffusion Rayleigh 
par les électrons X du mercure, pour des énergies 
de 0,32-0,64-1,28 et 2,56 mce?. Ces calculs ont une 
base théorique solide, en particulier, à la différence 
des calculs jusque là existants [3], ils utilisent pour 
les électrons dans l’état intermédiaire les fonctions 
d'ondes relativistes dans le champ du noyau. 

Brown et Mayers donnent pour 2,56 me?, 
soit 1,305 MeV, les parties réelles et imaginaires des 
amplitudes de diffusion pour les deux cas, avec et 
sans changement de polarisation. En tenant compte 


née dite nt ne RÉ) fes th LÉ EL LL D tn ne à à ne ne 
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de la diffusion Thomson, la section efficace est 
donnée par la formule : 


do 


HU + 7 + F9 [RE] 


2 


+R ne + 4 


F, et F; sont les parties réelles et imaginaires des 
facteurs de diffusion sans changement de polari- 
sation, F, et Fi correspondant à la diffusion avec 

IL m 
changement de polarisation, et 7° = 2? Ti étant le 
facteur de diffusion Thomson. Les valeurs indi- 


quées par Brown et Mayers sont 7" (5) 
n (L— cos 0 
et F 
L'étude des résultats théoriques conduit Brown 
et Mayers à constater deux propriétés importantes. 
Tout d’abord alors qu’il est nécessaire de calculer 
exactement les amplitudes de diffusion sans chan- 
gement de polarisation, c’est-à-dire les facteurs 7”, 
on s’aperçoit que les amplitudes avec changement 
de polarisation sont données avec une bonne pré- 
cision par le facteur de structure calculé par 
Bethe [3], que l’on peut mettre sous la forme para- 
métrique suivante : 


F(X) = SP ra o 


- KES 
Y Sin 2 


Nez 


E, énergie du + incident en MeV 
0, angle de diffusion 


et il suffit de remplacer (7% + iFi) par F, réel, 
dans la formule précédente. 

D'autre part, en comparant les résultats théo- 
riques pour 0,65 et 1,305 MeV, Brown et Mayers 
ont constaté que dans ce domaine d’énergie tout 
au moins, les facteurs Æ” ne dépendent que de 
Esin 0/2, comme les facteurs F et T (qui est une 
constante). Donc en traçant les courbes F(X) 
et F'(X) à partir des résultats à 1,305 MeV, on 
pourra calculer les sections efficaces dans ce 
domaine d’énergie. 

Nous avons tracé ces courbes pour le plomb, 
pour lequel nous avons déduit les amplitudes de 
diffusion à partir des résultats de Brown et Mayers 
dans le mercure en admettant comme eux une 
variation en Z5 des sections efficaces, pour un faible 
accroissement de Z. La précision obtenue ainsi est 
largement suffisante. 

Nous avons tracé sur les figures 4 et 5 la 
courbe F(X) du facteur de structure pour le plomb, 
ainsi que les courbes 


GUX) = r8(F, + T)2 + r8 Fi? 
Par = : 


| 
Ze 
= 
+- 
= 
Te 
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X=E sin 6/2 E en MeV 
(0) 0,5 1,0 X 
Fic. 4. — Facteur de structure pour Pb d’après Bethe, 


107%, 


cm? par stéradian 


0 0,5 1,0 X 
Hrcw5: 


de telle sorte que la section efficace est finalement 
donnée par la formule 


do 


do 1 + cos 0 
adQ 


9 


A 


1 — cos 7 


G(x) ( | + (— 


qui est sûrement valable pour le plomb et pour les 
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énergies comprises entre 0,65 et 1,305 MeV. Nous 
avons pu ainsi obtenir les courbes théoriques pour 
1,33 et 1,12 MeV (fig. 6 et 8). 
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F1c. 6. — Sections efficaces de diffusion élastique 
des rayons y de 1,33 MeV dans Pb. 
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F1G. 7. — Sections efficaces de diffusion élastique 
des rayons y de 1,33 MeV dans l’étain. 


Mais d’après Brown [4], le domaine de validité 
de la formule précédente peut être étendu à des 
énergies supérieures. (Cela revient à admettre que, 
même pour des énergies supérieures, les facteurs 
dépendent uniquement de Esin 60/2.) On peut ainsi 
calculer les sections efficaces de diffusion 
à 2,62 MeV, qui ont été utilisées dans une autre 


Nog7 


série d'expériences [5]. Mais la validité de la courbe 
obtenue est celle d’une extrapolation, alors que 
pour 1,12 MeV elle a la valeur d’une interpolation. 

Enfin, Brown et Mayers ont effectué des caleuls 
pour l’étain à 2,56 mc? en utilisant une formule 
semi-empirique établie à partir du facteur de struc- 
ture corrigé, et vérifiée tant que Z «x < 0,6. Nous 
avons porté ces résultats sur la figure 7. 

Tous ces calculs ne tiennent compte que de la 


- diffusion due aux électrons X. La contribution des 


électrons Z a été calculée approximativement par 
Brown et Mayers pour 0,650 MeV, en utilisant le 
facteur de structure et la formule semi-empirique 
pour les amplitudes avec et sans changement de 
polarisation, respectivement. À 1,305 MeV ils n’ont 
effectué la correction que sur l'amplitude avec 
changement de polarisation, et tiennent pour négli: 
geable la correction dans l’autre cas. D’après ces 
calculs, la contribution à la section efficace des 
électrons L ne dépasserait pas 20 %, aussi étant 
donné le caractère peu sûr de cette correction à 
1,305 MeV, nous avons tracé les courbes sans en 
tenir compte. 
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Fi, 8. — Sections efficaces de diffusion élastique 
des rayons y de 1,12 MeV dans Pb. 


5. Discussion et conclusion. —— Nous voyons sur 
les figures 6, 7, 8 qu’à la différence des autres expé- 
rimentateurs, dont les points sont en général situés 
bien au-dessus de la courbe théorique (1), nos 
mesures à 1,12 et 1,33 MeV coïncident avec les 


() G. E. Browx [2], loc. cut, 


dichéiet.sn.s" 
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prévisions de Brown et Mayers dans la limite des 
erreurs expérimentales, soit environ 20 %,. 

À 1,17 MeV, on retrouve sensiblement les résul- 
tats déjà obtenus [1] ce qui rend improbable une 
erreur de manipulation. Toutefois il y a un complet 
désaccord avec les prévisions de Brown et Mayers, 
ce qui montre que la méthode indirecte décrite 
dans [1] n’est pas valable. On ne peut expliquer 
ceci, ni par les empilements résiduels dus à la 
diffusion Compton, dont nous avons tenu compte 
expérimentalement, ni par le rayonnement de frei- 
nage des électrons projetés à l’intérieur de la cible 
car cette méthode indirecte tenait compte de cet 
effet. On peut interpréter ce désaccord par l’exis- 
tence d’une diffusion inélastique de la raie 
de 1,33 MeV, effet négligé dans notre mesure 
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à 1,17 MeV. Cependant son intervention possible 
dans nos expériences fut envisagée par Peierls, et 
d’après une étude théorique récente de Randles [6], 
elle ne serait pas négligeable. 

Étant donné la sûreté des calculs de Brown et 
Mayers, si l'effet Delbrück contribue à la diffusion 
pour des énergies inférieures ou égales à 1,33 MeV 
ce ne peut être que pour une faible part. Pour le 
mettre en évidence il faudrait améliorer à la fois la 
précision des expériences et celle du calcul, notam- 
ment en ce qui concerne les électrons L. 

Nous remercions M. E. Alexandre pour l’aide 
efficace qu’il a apportée dans la réalisation du dis- 
positif et dans les manipulations. 


Manuscrit reçu le 127 février 1958. 
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FORMATION DE L'IMAGE D'UN POINT 
DANS UN SPECTROGRAPHE MAGNÉTIQUE UTILISANT UN CHAMP EN r-! 


Par C. BASTARD et J. LAFOUCRIÈRE, 
Institut de Physique Nucléaire de l'Université de Lyon. 


Résumé. — On étudie les trajectoires du type trochoïdal de corpuscules chargés dans un champ 
magnétique inhomogène à rH constant en vue de la construction d’un spectrographe 6 à grande 


luminosité. 


L'étude des trajectoires gauches montre que l’image d’un point est représenté pes un segment 
de droite perpendiculaire au plan de symétrie du champ. Ë 


Abstract. — À study of the paths of trochoïdal type described by charged particles in an inho- 
mogeneous magnetic field of constant rÆ, with a view Lo the construction of a high-luminosity B-spec- 


tr ograph. 


The study of the twisted paths shows that the image of a point is represented by a straight 
segment perpendicular to the plane of symmetry of the field. 


L'intérêt des études de spectrographie 8 n’est 
plus à démontrer en physique nucléaire. Les deux 
types de spectrographes, les plus utilisés, spectro- 
graphes à lentilles et spectrographes plats, pos- 
sèdent des qualités différentes. Les premiers ont 
une bonne luminosité, les seconds sont- surtout 
intéressants par leur pouvoir séparateur. L'appareil 
étudié ici possède l’avantage de la luminosité des 
lentilles et de la définition des prismes. Il utilise un 
champ magnétique de révolution, possédant un 
plan de symétrie, tel que le produit r A soit cons- 
tant, r étant la distance à l’axe de révolution et H 
la valeur du champ dans le plan de symétrie. Pré- 
cisons, dès maintenant, la différence essentielle 
avec le spectrographe dit à «double focalisation» 
dû à Siegbahn. Dans ce spectrographe les trajec- 
toires des particules sont les orbites stationnaires ; 
dans notre cas, les orbites stationnaires ne sont pas 
stables au point de vue radial et ne peuvent donc 
pas être utilisées. 


I. Trajectoires planes. — Le corpuscule de masse 
au repos M9, de charge e, placé dans un champ 
magnétique à symétrie de révolution défini par le 
potentiel vecteur À, — À, —0 4, — A(rz) décrit, 
comme il est bien connu, des trajectoires gauches 
données par le système 


mr” = mr og"? +eo ; (r A) 


| mr? © + erA = C ‘(1) 
| my = e s+ (rA). 
Soit en introduisant l'intégrale première 
v? — constante : 
Rae AR — v8 
mr? ®" + erA = C (2) 
UE = ri (rA) 


où m est la masse relativiste m1 — (v?/c2)F 22, 
Si z — 0 le corpuscule reste dans le plan de 
symétrie du champ magnétique et le mouvement 


décrit par les deux équations 


TE Tr QE 72 vi à 
mr? po + erÀA = C (3) 


s’efféctue suivant des trajectoires trochoïdales [1] 
ou des trajectoires présentant des branches infinies. 
Seules les premières nous intéressent, nous verrons 
par la suite la condition pour obtenir ces trajec- 
toires. 

La forme particulière du champ magnétique 
rend d’une part le problème complètement inté- 
grable et donne lieu, d’autre part à une focalisation 
radiale parfaite [2]. 

Le potentiel vecteur décrivant un tel champ 
magnétique est solution de l'équation de Laplace : 
D (A) 5 D (PA) + À (r4) = 0 

dr? 22° 


La fonction rA(rz) définie par rA(r0) — Xr dans 
le plan de symétrie vérifiant l’équation (4) est 


z? 1/2 
rA(rz) = Kr (1 +5) À 


Les lignes de force du champ sont done des 
cercles [3]. L’équation des trajectoires planes des 
corpuscules chargées est dans ce cas : 


Pis pe 
r (C — eKr) dr 


ue ro Cl(m?v — e? K?) r? + 2eKCr — COUES 


Introduisons les quantités suivantes : 


Ée , GC C." ER 
eK + mv, enr. TRE LE TER 


ri r = 
à eK — ms eK UTA 


Les trajectoires trochoïdales sont obtenues pour 


te Es dé RE ie 
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k > 1, Péquation de ces trajectoires s’écrivant : 


k(rm — r) dr 
RSS De [re 7) (re 


intégrale qui s’intègre facilement. 
La précession à la valeur simple 


Ô = 2 (1 > 


mx) 


_ cette quantité ne dépendant que de # = K}Jp H, 


avec pA — mivle la focalisation a lieu quel que 
soit l’angle de tir «. D’autre part, le pouvoir sépa- 
rateur qui peut être défini par l'expression 
do/dk —2r/(k?— 1) *? peut être pris aussi grand 
que l’on veut en choisissant k voisin de 1. 


| 

SK S | 
RCE 
NS 

SJ 
FiG:21. 


L'aspect des trajectoires est le suivant (fig. 1): 
Toutes ces trajectoires sont homothétiques par 
rapport au centre de symétrie du champ dans le 

SRE : 
rapport : : = — = . 7 “et «' étant les angles 
* de tir de deux corpuscules. 


II. Trajectoires gauches. — Le système [2] ne 
se ramenant pas à des quadratures, nous avons eu 
recours à une méthode d’approximation en nous 
bornant à des trajectoires voisines du plan z — 0. 

La troisième équation du système (2) s’écrit en 
remplaçant ®' par sa valeur 

# e C— erA DA 

ST r dz 

C — erA dA 1 
— — en série de 
dz 
Taylor au voisinage du plan de symétrie en négli- 
geant les termes du second ordre 
C—erAdA  fC—erADA 
r Dz a ( L ae 


Nous développons alors 


il. 


+:| 
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On obtient finalement : 


ze (- — %) (5) 
de m? r DANCE CEE 


soit plus simplement : 


r r 
PT ee PA 
BI NR AUS EE F3 


équation du second ordre en {, r étant une fonction 
de t obtenue en intégrant le système (3). 
Passant à la variable r on obtient : 


d2z 


k°? dz 
“e 7 (r AT Fm) (ru 


r) —=\ 0; 


Cette équation est l’équation différentielle des 
méridiennes de révolution sur lesquelles se trouve 
située chaque trajectoire gauche. C’est une équa- 
tion linéaire du second ordre qui n’est pas d’un 
type connu et ne se ramène pas à la forme hyper- 
géométrique. Si on exprime les conditions initiales 
par les coordonnées r, 0, ©, de la source, 0, angle 
du plan de symétrie avec la vitesse initiale, w, angle 
de la tangente au cercle de rayôn r, avec la pro- 
jection de la vitesse initiale sur le plan de symétrie, 
la forme des coefficients de l’équation (6) montre 
que les méridiennes correspondant aux mêmes r, 
et 0, mais à des w, différents, sont les mêmes. Nous 
pouvons donc prendre un %, particulier uw, — 900. 

Pour construire les méridiennes nous avons 
résolu numériquement l’équation (5) 


Nous avons utilisé une formule polynominale due 
à Milne [4], 


LA 
2% 
Zn +1 22h  En—1 + L2n 


t étant l'intervalle de temps entre deux points 
consécutifs, le 17 point étant obtenu en supposant 
le mouvement uniforme, ce qui est légitime puisque 
= Débré 0. 

La précision donnée par cette formule a été 
jugée satisfaisante en l’appliquant à la construction 
de la trochoïde plane qui était calculable eu moyen 
de transcendantes élémentaires de l’algèbre. 

Le calcul montre que pour des angles 0, allant 
jusqu’à 100 la méthode précédente était applicable. 
L'aspect des méridiennes est alors la suivante. 
Nous avons pris 


Fo — 20cm pH — 10%, k — 2; (fig. 2). 


Nous observons une concentration au voisinage 
des points z —0 pour € = 13,5 em il y a donc 
deux focalisations : une focalisation radiale et une 
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focalisation axiale, mais les angles de focalisation 
ne sont pas égaux et la seconde se produit pour 


TR ATIE 


HiGeu2 


0 5° : 10° 
Fic. 3. 
Pour obtenir la forme de l’image d’un point 
source, il reste à étudier la variation de la préces- 
sion à en fonction de l’angle 0,. Pour cela nous 


avons résolu numériquement le système (1) en utili- 
sant comme précédemment les deux formules 


#” 
Ê À 2 
Zn+1 = 23n — 2n—1 + l° 


» 
c 22 
Tntai = 27n — ni + L° 2. 


Nous avons trouvé que pour les trajectoires 
voisines du plan de symétrie, à reste constant 
lorsque 0, varie, l’image d’un point source S situé 
sans le plan z — 0 est de ce fait un petit segment 
de droite perpendiculaire au plan z — 0 et situé sur 
le cylindre r — r, à la distance angulaire dde S. 

En reprenant l'exemple considéré plus haut 
ro — 20 em, pH — 10, À — 2 les points d'impact 
des trajectoires gauches se situent sur ce segment à 
des distances z données en fonction de 6, par la 
courbe (3). 

Le domaine de validité de l’équation (5) s’étend 
ainsi jusqu'à 0, æ 100. L’allure des trajectoires 
gauches est alors la suivante : 


Fic. 4, 


Nous ne saurions terminer sans adresser nos 
remerciements à M. le Pr J. Thibaud, Directeur de 
PInstitut de Physique Nucléaire, pour l'intérêt 
qu’il a bien voulu prendre à notre travall. 


Manuscrit reçu le 7 mai 1958. 
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LETTRES A LA RÉDACTION 


INFLUENCE DE L’INTRODUCTION 
D'IMPURETÉS (Cu, Zn, Li) 
SUR LES PROPRIÉTÉS DIÉLECTRIQUES 
DE L’OXYDE DE ZINC 


Par Mlle. Marie-Louise BLANCHARD 
et Mille Maryvonne MARTIN, 


Faculté des Sciences de Rennes. 


Le comportement diélectrique de l’oxyde de zime 
dans le spectre hertzien est très complexe [1]et s’avère 
particulièrement sensible aux traces d’impuretés [2]. 
La position des bandes d’ absorption dipolaire Debye 
et l'énergie d'activation qui s’en déduit sont fortement 


“influencées par la concentration du cation introduit. 


Nous avons repris dans de meilleures conditions 
expérimentales que précédemment [2], [3] l’étude 
d’oxydes différemment dopés : le traitement sous des 
atmosphères variées (azote, oxygène) a permis de pré- 
ciser les résultats que nous allons résumer ici. 


Préparation des échantillons, — L'introduction 
d’impuretés s'effectue par impréenation de la poudre 
d’oxyde de zinc (Vieille-Montagne) [1] par des solutions 
de dilution variable, séchage à létuve à 130 °C 
(2 heures) puis traitement vers 900 °C en atmosphère 
d'azote (4 heures). Les échantillons sont refroidis 
rapidement (trempe) après leur traitement. 


I. ZnO.popé AU cuIvRE. — Les courbes €” — f(T) 
sont représentées pour v — 10 KHz sur la figure 1 


£" 


50 100 150 200 T°K 
FiG. 1. — ZnO : Cu. 
v = 40 kHz ; concentrations : 1. — 107; 124— 41072; 
3. — 10—4: h,— A0 ; — 5, 5 10—$, 


pour les diverses concentrations utilisées. Les maxima 
se déplacent vers les basses fréquences lorsque la 


concentration diminue. De la position des bandes aux 
différentes fréquences, on peut déduire les énergies 
d'activation correspondant aux différentes teneurs en 
impuretés. 


ÉNERGIE D’ACTIVATION 


CONCENTRATION SUN 
107$ 0,07 
40 0,21 
ND 0,23 
10— 0,24 
94072 0,45 


) 


La comparaison des valeurs de U suggère l’existence 
de tros domaines : 


1) Concentrations voisines de 106 : U — 0,07 eV. 

Il est difficile de préciser les résultats pour des con- 
centrations inférieures à cette limite : Quels que soient: 
les soins apportés à la préparation, on introduit néces- 
sairement une concentration en impuretés de cet ordre 
de grandeur. 


2) Concentrations comprises entre 105 et 103 cor- 
respondant à des énergies moyennes de l’ordre de 
0,23 eV. 


3) Concentrations supérieures à 10% 
de l’ordre de 0,45 eV. 


Nature du dipôle. 

Un échantillon ZnO (104 Cu) a été traité à 850 0C en 
atmosphère d'oxygène : aucune bande d'absorption 
rapparaît. Ceci prouve que la présence de lacunes 
d’anions est indispensable à la formation du dipôle, 
celui-ci étant vraisemblablement constitué par l’asso- 
ciation cuivre-lacune. 

Lorsque l’on dépasse une concentration en impureté 
supérieure à 10—%, on passe à un mécanisme à énergie 
plus élevée. Il semble que la quantité de vacances 
d’anions disponibles ne soit plus suffisante pour créer 
des associations du type précédent. On observe en 
effet, l'apparition de ce mécanisme pour des concen- 
trations bien inférieures en impureté, si l’on effectue le 
traitement thermique dans des conditions moins favo- 
rables à la formation de vacances d'oxygène, par 
exemple en opérant dans l'air. 

Pour ces concentrations élevées en impuretés, la 
compensation pourrait être réalisée non plus par 
vacances d’anions, mais par piégeage de charges sur 
les atomes avoisinant l’impureté [4]. 


: forte énergie, 


II. ZnO popé AU zINc. — Les courbes obtenues 
pour y — 10 kHz sont représentées figure 2 Pour 
chaque concentration, il apparaît deux domaines 
d'absorption : L’un aux basses températures (< 700K) 
correspondant à une faible énergie, l’autre aux tempé- 
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ratures moyennes (entre 100 et 200 °K) correspondant 
à une énergie plus élevée. 


F1G. 2. — ZnO : Zn. 


v = 10 kHz ; concentrations : 4. — 106; 2:— 10; 
D l0 mie l0RES 


On remarque que le niveau de faible énergie 
s’affaiblit au profit de celui à haute énergie lorsque la 
concentration croit. 

Le tableau ci-dessous donne les énergies d’activation 
des deux régions en fonctions des concentrations. 


ÉNERGIE 
CONCENTRATION D’ACTIVATION U 

EN eV 
407$ environ 0,02 
; 0,18 
d'OS environ 0,09 
10e environ 0,10 
0,24 
1072 0,27 

III. 7nO popé AU LITHIUM. — Le dopage par ions 


lithium se manifeste par l'apparition de bandes 
d'absorption Debye situées au-dessous de 100 0K 


200 MES 


50 100 150 


FrG. 3. — ZnO ; Li. 


v — 10 kHz ; concentrations : 4. — 10—$ ; 
2, — 104; 3.:+ 10—3. 


(fig. 3) dont l’énergie d’activation reste très faible par 
rapport aux cas précédents. On peut toutefois soup- 
çonner la présence de plusieurs mécanismes d’absorp- 
tion analogues à ceux précédemment observés dans le 
cas du cuivre : un premier niveau d’énergie inférieur 
à 0,01 eV et un second de l’ordre de 0,030 eV. Pour les 
fortes concentrations (fig. 3) le phénomène s’affaiblit 
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considérablement et la bande L'ebye est noyée dans une 
conductivité correspondant vraisemblablement à un 
nouveau mode de compensation de l’impureté. 


Lettre reçue le 3 juin 1958. 
BIBLIOGRAPHIE 


[1] BzancxARD (M. L.), C. R. Acad. Sc., 1957, 244, 767. 

[2] HAGÈNE (B.), CorneTeau (H.) et BLancuarD (M. L.), 
J. Physique Rad., 1956, 17, 450. 

[3] FREyMANN (R.) et coll., Colloque luminescence, Paris, 
mai 1958. J. Physique Rad., 1956, 17, 806. 

[4] Krocer (F. A.) et Vinx (H. J.), Physica, 1954, 20, 950. 


NS 


MISE EN ÉVIDENCE 
DE LA DIFFUSION INÉLASTIQUE DES RAYONS y 
PAR LES ÉLECTRONS LIÉS DE L’ATOME 


Par L. GoLpzanL et J. BANaAïrGs, 


Laboratoire de Physique Atomique et Moléculaire 
du Collège de France. 


I. Introduction. — Au cours d’expériences effec- 


tuées dans ce laboratoire sur la diffusion élastique des 
rayons y du %Co et du Zn [1], la comparaison des 


spectres diffusés a suggéré l’existence d’un processus 


non encore étudié de diffusion des rayons y. Il s’agit 
de la diffusion inélastique des rayons y par les électrons 
liés de l’atome, diffusion dans laquelle l’énergie et la 
quantité de mouvement du photon incident se répar- 
tissent entre l’atome, l’électron diffuseur et le photon 
diffusé. Par suite de cette répartition en trois corps, on 
a un spectre continu pour l'énergie des rayons y dif- 
fusés sous un angle donné. L'énergie maximale de ce 
spectre est pratiquement celle des rayons y diffusés 
élastiquement. 

Randles [2] a effectué un calcul approximatif du 
rapport de la section efficace de cet effet à celle de la 
diffusion élastique, dans une région voisine de celle-ci, 
pour des rayons y de l’ordre de 1 MeV, et pour de 
grands angles de diffusion. D’après ce calcul, ce rapport 
ne serait pas négligeable. 

Nous avons entrepris de mettre en évidence ce phé- 
nomène pour les rayons y de 1,12 MeV du Zn, diffusés 
sous 509 par du plomb, jusqu’à 0,9 MeV environ soit 
dans une bande de 20 % de l'étendue totale du spectre. 


2. Dispositif. — On a utilisé un dispositif de révo- 
lution du même type que ceux employés dans les 
expériences sur la diffusion élastique [1]. 

Une source de 5Zn à 0,5 C environ irradiait la cible 
étudiée et le rayonnement diffusé était mesuré à l’aide 
d’un détecteur à scintillations constitué d’un cristal 
de Nal(T1) de 4 X 5 cm (des mesures effectuées avec 
un cristal de 10 X 7,5 em ont donné de moins bons 
résultats), d’un photomultiplicateur Du Mont 6292, et 
d’un sélecteur R.I.D.L. à 100 canaux réglé sur la partie 
duspectre étudiée. Le détecteurétait protégé parun filtre 
de plomb de 2,5 cm d’épaisseur, afin de réduire le taux 
des empilements. Le dispositif électronique utilisé 
dans ce but au cours des expériences précédentes n’a 
pas été employé ici, car il risque de déformer le spectre 
aux basses énergies et nous verrons que par la méthode 
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No: 


des filtres, nous avons suffisamment réduit le taux des 
empilements. 

On a utilisé 4 cibles différentes : 

— Une cible cylindrique de plomb de 27 cm de 
diamètre, 18 em de hauteur et 1/2 mm d'épaisseur, 
qui nous a servi pour les mesures. 

— Une cible formée de feuilles d’or séparées par une 

- épaisseur calculée de papier, qui a servi pour des véri- 
fications (voir paragraphe suivant). 

— Une cible formée d’un anneau de cuivre de 8 mm 
d'épaisseur et de 24 cm de diamètre qui nous a servi à 
mesurer l’efficacité du cristal. 

— Une cible d'aluminium de même forme que la 
première pour l’évaluation des empilements. 


3. Mesures. — Le spectre diffusé est la résultante 
de nombreux processus dont il faut tirer la partie 
correspondant au phénomène étudié. 

Chaque atome de la cible donne lieu à 3 processus 
de diffusion : 

— La diffusion élastique (effets Thomson, Rayleigh, 
Delbruck, résonance nucléaire. L’atome diffuseur 
prend tout le recul). Elle peut être éliminée du spectre 
diffusé par soustraction d’un spectre fournit par irra- 
diation directe du détecteur, et normalisé à la hauteur 
du pic de diffusion élastique. 

— La diffusion Compton (l’électron diffuseur prend 
tout le recul). Cette diffusion se fait à 50° avec une 
grande perte d'énergie et peut donc être facilement 
.discriminée, mais elle est la cause essentielle des empi- 
lements. 

— La diffusion inélastique étudiée. 

D'autre part, l'épaisseur même de la cible étudiée 
donne lieu à d’autres processus de diffusion : 

— Des diffusions multiples. — Diffusion sous un 
faible angle par effet Compton dans la cible du rayon- 
nement inélastique lui-même, ou du rayonnement déjà 
diffusé élastiquement. Ces effets ont été éliminés ici 
par l’adoption de cibles suflisamment minces pour 
rendre la probabilité de double diffusion très faible. 


— Le rayonnement de freinage des électrons projetés. 
— Si l’on se place, comme dans notre cas, près du pic 
de diffusion élastique, il ne peut s’agir que des élec- 
trons produits par effet photoélectrique dans la cible. 
Ce rayonnement parasite est alors dû à la succession 
de deux processus : L’effet photoélectrique et le rayon- 
nement de freinage des électrons projetés. 

Lors des expériences de diffusion élastique 
à 1,33 MeV [3], la comparaison des spectres diffusés 
par une cible de plomb pur, et une cible constituée par 
une solution de nitrate de plomb, a montré que le 
rayonnement de freinage ne jouait aucun rôle à 500 
jusqu’à des énergies inférieures de 20 % à celle de la 
diffusion élastique. Cette partie du spectre correspond 
à celle étudiée dans notre expérience. 

Nous avons vérifié ceci par une autre méthode : 
Appelons a l’épaisseur d’une cible telle qu’un électron 
projeté la traversant y perde 20 % de son énergie (ce 
qui la réduit à une valeur inférieure à celles du domaine 
étudié). 

Dans ces conditions, pour une cible d'épaisseur e 
plus grande que a, l'intensité du rayonnement para- 
site Z,”est proportionnelle à a.e, alors que pour une 
cible d'épaisseur inférieure à a, cette intensité est pro- 
portionnelle à e2. Comme dans les deux cas l’intensité 
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du rayonnement étudié /, est proportionnelle à e, le 
rapport /,/1, est constant dans le premier cas, et 
diminue comme l’épaisseur de la cible dans le second 
cas. 

Donc, si la comparaison des deux spectres norma- 
lisés sur le pic de diffusion élastique et obtenus respec- 
tivement avec une cible d'épaisseur plus grande que 4, 
et une cible d'épaisseur beaucoup plus faible ne montre 
pas de différences, c’est que le rayonnement parasite a 
une importance négligeable par rapport au rayon- 
nement diffusé (dans la zone étudiée), même pour une 
cible épaisse. 

Pour l’or et des électrons de 1 MeV environ, a est de 
l’ordre de 0,1 mm. Nous avons donc utilisé une cible 
en or de 0,02 mm d’épaisseur. Mais le taux de comptage 
était alors très faible. Nous avons donc été amenés à 
prendre une cible constituée de plusieurs feuilles d’or 
séparées par des feuilles de papier d'épaisseur telle 
qu’un électron de 1 MeV y perde 20 % de son énergie, 
sans produire de rayonnement de freinage appréciable, 
par suite des faibles numéros atomiques des consti- 
tuants du papier. Le rapport /,/1e est alors le même 
que pour une feuille d’or unique. 

Cette méthode a confirmé que dans la zone étudiée, 
le rayonnement de freinage ne jouait aucun rôle. 


DÉFORMATION DU SPECTRE DUE AU FILTRE DE PLOMB 
ET AU DÉTECTEUR. — Le filtre de plomb de 2,5 cma 
suffi pour éliminer complètement les empilements, ce 
que nous avons vérifié par la méthode des cibles équi- 
valentes déjà décrite [3]. Mais ce filtre provoque une 
dégradation du spectre diffusé et modifie l'efficacité du 
détecteur. 

La méthode de soustraction de la diffusion élastique 
tient compte de cet eflet, puisque le spectre direct 
normalisé est obtenu avec le même filtre. Par contre, 
il n’a pas été tenu compte de la dégradation du spectre 
de la diffusion inélastique. Mais la forme décroissante 
du spectre de cette diffusion et le fait qu’une raie 
de 1 MeVenviron est peu déformée par la présence du 
filtre, nous conduisent à penser que la dégradation 
produite est de faible importance. 


VARIATION DE L'EFFICACITÉ DU DÉTECTEUR MUNI DE 
SON FILTRE. — Nous avons mesuré la variation d’efti- 
cacité du détecteur (avec son filtre) en utilisant les 
variations d’énergie en fonction de l’angle du rayon- 
nement diffusé par l’effet Compton produit sur une 
cible par les rayons y de 2,76 MeV du ?#Na. On a tenu 
compte de la variation d'intensité de ce rayonnement à 
l’aide de la formule de Klein-Nishina. 

Afin de rester dans les mêmes conditions expéri- 
mentales, nous avons utilisé le même dispositif avec 
une cible et des collimateurs permettant d’obtenir 
une résolution angulaire inférieure à 3°. 

On trouve que l'efficacité entre 1,12 MeV et 0,9 MeV 
est pratiquement constante. Cela est dû en partie à ce 
que le filtre compense la variation d'efficacité du cristal. 


INVERSION DU SPECTRE. — Le spectre finalement 
obtenu pour la diffusion inélastique n’est pas le spectre 
réel car le détecteur introduit une transformation due 
à l’aspect statistique de son fonctionnement. 

Soit g(Æ) le spectre continu réel de diffusion inélas- 
tique ; on a après le détecteur, en admettant une dis- 
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tribution gaussienne pour les électrons de la photo- 
cathode, un spectre mesuré de la forme 
(EE) 


fP atBet) ——e #6 
== oi e 2k LU 
DONNE oÈ 


où £E est l’écart quadratique moyen pour la courbe de 
résolution correspondant à l’énergie Æ du photon, 
et o(E) l'efficacité du système cristal plus filtre. - 

On ne connaît pas de solution mathématique exacte 
de ce type d’équation intégrale (équation de Fredholm 
de 1re espèce), mais si l’on fait les approximations 
suivantes : 

— kJE est petit (bonne résolution). 

— Dans le domaine correspondant aux valeurs ds la 
distribution de. Gauss contribuant notablement au 
résultat, on peut utiliser un développement 
de g(Æ).p(E) limité au second ordre, on obtient la 
formule : 


Dans nos expériences b(E) varie de façon négligeable, 
et l’on a : 


où E/2 vaut 0,02. 


4. Résultats. — Nous avons tracé sur la figure : 

— Le spectre diffusé (déduction faite du bruit de 
fond). 

— Le spectre direct normalisé. 


— Leur différence, représentant le spectre de dif- 


fusion inélastique donné par le détecteur. 

— Le spectre réel de cette diffusion. 

On voit donc que pour les rayons y de 1,12 MeV 
diffusés par le plomb sous 50° (avec une résolution 
angulaire de 80), il y a approximativement autant 
de photons diffusés inélastiquement entre 1,12 et 
0,95 MeV que de photons diffusés élastiquement. 

Nous remercions M. Philippe Eberhard, pour de 
fructueuses discussions, M. Nguyen-Huu, pour son 
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Comptage 


x Spectre mesuré 
\ Ds 


Spectre de diffusion : 
inélastique mesuré 0 


eu CS 


Spectre réel de diffusion 
inélastique 


03 5 


Fic. 1. — Diffusion inélastique 
des rayons y de 1,12 MeV par le plomb. 


aide dans la réalisation de la partie électronique, et NW 
M. Eugène Alexandre pour sa précieuse collaboration ! 
dans le dessin et la réalisation du dispositif méca- 
nique. 

Lettre reçue le 10 juin 1958. 
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